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Abstract 


In this thesis, the physical properties of the alkaline earth oxide, MgO, doped with 
transition metals and co-doped with a rare earth element (Eu) and a transition metal are studied. 
The approach utilized 1s the density functional theory, with the use of the full-potential linearized 
augmented plane waves (FP-LAPW) method. This method is implemented in the Wien2k solid 
state simulation code. The calculations are performed using the GGA and GGA+mBJ 
approximations for the exchange and correlation potential. This methodological approach allows 
us to obtain accurate results for the structural, electronic, magnetic and optical properties of doped 


and co-doped M£0O. 


In the first section, a thorough study is conducted on MgO oxide using the main 
approximations incorporated in the code (LDA, GGA). This approach aims to compare and validate 
the results obtained against experimental data. Then, doping M£O with 12.5% of Ni, Pd, or Ptin 
the NaCI structure was studied using the GGA-PBEsol approximation. The effect of Z and d 
electrons on the properties of the doped compounds was analyzed. The X0.125Mg0.8750 (x=Ni, Pd, 
Pt) compounds were found to be magnetic semiconductors with a narrower band gap than M£gO. 
They also had a magnetic moment of 2 18 and a high absorption in the low UV regions, making 


them promising for optical applications. 


In a second phase of the study, the effect of doping MgO with Cu, Ag, or Au at 
concentrations of 3.125%, 6.25%, 12.5%, and 25% was studied. The results showed that the nature 
of the doped compounds changed with concentration. At 25%, the compounds were non-magnetic 
conductors. At low concentrations, they were ferromagnetic semi-metals with a magnetic moment 


of 1 u8 and a strong absorption in the UV. 


The last part of the study, MgO was co-doped with Eu and Cu, Ag, or Au at a 
concentration of 12.5% for each element. The co-doped compounds were found to have a narrower 
band gap and a preserved semi-metallic character compared to the pure and doped oxides. The 


optical properties of the co-doped compounds were also improved. 


This process, doping and co-doping of MgO with Eu and transition metals, opens the way 
to a wide range of potential applications in the fields of optoelectronics, spintronics and magnetic 


data storage. 


Keywords: 
Me£0, spintronics, doping, co-doping, transition metals, magnetic properties, optical 


properties, Wien2k code, rare earth co-doping. 


Résumé 


Dans le cadre de cette thèse, les propriétés physiques de l’oxyde alcalino-terreux, MgO, 
mono-dopé avec des métaux de transition et co-dopé avec un élément des terre rare (Eu) et un métal 
de transition sont étudiées. L'approche employée est la théorie de la fonctionnelle de la densité, 
avec l'utilisation de la méthode des ondes planes augmentées linéarisées à plein potentiel (FP- 
LAPW). Cette méthode est implémentée dans le code de simulation des solides Wien2k. Les 
calculs sont réalisés à l’aide, des approximations GGA et GGA+mB)J, pour le potentiel d’échange 
et de corrélations. Grâce à cette démarche méthodologique, des résultats précis des propriétés 


structurales, électroniques, magnétiques et optiques de M£O dopé et co-dopés sont obtenus. 


Dans la première section, une étude approfondie est menée sur l’oxyde MgO en utilisant 
les approximations principales intégrées dans le code (LDA, GGA). Cette démarche vise à 
comparer et valider les résultats obtenus par rapport aux données expérimentales. Par la suite, 
l’approximation GGA-PBEsol est adoptée pour étudier l’effet du dopage de MgO par des métaux 
de transition (Ni, Pd et Pt) à 12.5% de concentration dans la structure NaCT. L'influence du numéro 
atomique (Z) et du nombre d'électrons (d), sur les propriétés des composés dopés est analysée. 
X0.125M g0.8750 (x=N1, Pd, Pt) présentent un comportement semi-conducteur magnetique dans les 
deux directions de spin (up et down) avec un gap énergétique plus étroit que la matrice hôte et un 
moment magnétique de 2 ug indiquant un comportement ferromagnétique. Ils montrent également 
une absorption élevée dans les régions basses UV par rapport au composé pur, ce qui les rend 


prometteurs pour des applications optiques. 


Dans une seconde phase de l'étude, l'impact de la variation de la concentration (3.125%, 
6.25%, 12.5%, et 25%), en corrélation avec le changement du numéro atomique (Z) et de la 
configuration électronique (d) des métaux de transition TM;M£g:1:0 (Cu, Ag, Au) est examiné. 
L'étude montre une transition de nature avec la concentration. À 25%, les composés dopés sont 
conducteurs non magnétiques. A faibles concentrations, ils deviennent semi-métalliques 


ferromagnétique (1 LB) avec une forte absorption dans l’UV. 


La dernière partie de l'étude se concentre sur le co-dopage de MgO avec les mêmes métaux 
de transition (Cu, Ag, Au), en introduisant un élément des terres rare (Eu) à une concentration fixe 


(12.5 % pour Eu et 12.5 % pour TM). Ce co-dopage a conduit à un rétrécissement plus important 


de la bande interdite et à la préservation du caractère semi-métallique de tous les composés. De 


plus, les propriétés optiques sont améliorées par rapport à l'oxyde pur et dopé. 


Ce procédé ouvre la voie à un large éventail d'applications dans les domaines de 


l'optoélectronique, de la spintronique et du stockage de données magnétiques. 
Mots clés : 


Me£0, spintronique, dopage, co-dopage, métaux de transitions, propriétés magnétiques, propriétés 


optiques, code Wien2k, co-dopage aux terres rares. 
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Introduction générale 


Depuis le début du XIXe siècle, les semi-conducteurs ont attiré beaucoup d’attention en 
raison de leurs applications diverses, s'insérant dans la plupart des industries électroniques et 
démontrant leur importance. Généralement, les études sur les semi-conducteurs reposent sur la 
manipulation du transport de charges. D’un autre côté, la spintronique utilise le spin de l’électron 
au lieu de la charge. La combinaison des propriétés de charge et de spin des électrons offre aux 
semi-conducteurs des avantages indéniables, tels qu'une faible consommation d'énergie réduite et 
une vitesse de traitement des données exceptionnelle, positionnant ces matériaux au cœur des 


innovations technologiques contemporaines. 


L’une des avancées technologiques les plus performantes de tous les temps est la 
technologie de stockage de données. Étant donné que la polarisation magnétique demeure stable 
dans les circuits magnétiques, le spin électronique est utilisé dans le stockage de données 
magnétiques pour enregistrer les informations dans des matériaux magnétiques plutôt que d'utiliser 
des charges électriques ou des flux de courant. Cette approche garantit la rétention des données 
même en cas de panne de courant. En raison de ces avantages significatifs, l’électronique basée sur 
le spin émerge comme l’un des domaines d’étude les plus cruciaux pour le développement de 


nouveaux dispositifs fonctionnels à semi-conducteurs. 


La matière joue un rôle essentiel et incontournable dans notre vie quotidienne, ainsi que 
dans le développement de la technologie dans laquelle nous vivons. Les laboratoires scientifiques 
rivalisent pour produire ces matériaux et comprendre leur structure cristalline. Les chercheurs ont 
amélioré les propriétés de ces matériaux et ont contribué à l'émergence de nouvelles familles aux 
applications variées, générant ainsi une richesse industrielle considérable. Ces matériaux 
acquièrent des propriétés spécifiques en fonction de divers paramètres. Ce niveau de précision dans 
la manipulation des propriétés matérielles ouvre la voie à des avancées technologiques 


significatives. 


Une nouvelle génération de matériaux spintronique, appelés semi-conducteurs 
magnétiques dilués, qui suscite un grand intérêt de la communauté scientifique pour ses 
applications potentielles. Les composés binaires II-VI substitués par des métaux de transition 
possèdent des propriétés intéressantes et uniques, qui peuvent offrir des applications 


technologiques potentielles dans les dispositifs spintronique. L'utilisation de métaux de transition 
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pour substituer certains éléments dans des oxydes semi-conducteurs, tels que ZnO, In203, CeO et 
TiO», offre des propriétés magnétiques intéressantes [1-4]. Ces oxydes semi-conducteurs dopés 
sont considérés comme des candidats prometteurs pour les applications spintronique. L’une des 
caractéristiques intéressante observée dans certains oxydes magnétiques dilués est la demi- 
métallicité ferromagnétique. L'idée de ce concept est proposée pour la première fois par de Groot 
et coll. [5, 6]. La demi-métallicité se réfère à la coexistence d'états électroniques métalliques et non 
métalliques à la limite de Fermi, la surface d'énergie où l'on trouve les niveaux d'énergie des 
électrons à zéro kelvin. Elle peut être exploitée pour contrôler le flux de spin des électrons, ce qui 
est crucial dans le domaine de la spintronique où l'orientation du spin est utilisée pour transmettre 


et stocker l’information. 


Dans le contexte des DMS, des études ultérieures ont révélé la présence de 
ferromagnétisme semi-métallique dans différents types de composés qu'ils soient dopés ou non 
dopés. Cela inclut des matériaux tels que les semi-conducteur magnétique dilué eux-mêmes, ainsi 
que les composés de type heusler, les nitrates, les oxyde de type pérovskite et les oxydes de métaux 
alcalino-terreux. Il est intéressant de noter que les oxydes de métaux alcalino-terreux dopés aux 
métaux de transition, aux non-métaux et aux métaux des terres rares peuvent présenter des 
propriétés du ferromagnétisme semi-métalliques [7-10]. Le ferromagnétisme semi-métallique est 


une caractéristique qui combine des propriétés métalliques et ferromagnétiques. 


Parmi les oxydes semi-conducteurs mentionnés se trouve le semi-conducteur MgO, sa 
transparence aux longueurs d'onde visibles, associée à ses propriétés semi-conductrices, en fait un 
matériau attrayant pour des applications telles que les écrans plats, les diodes électroluminescentes 
(LED) et d'autres dispositifs optoélectroniques. En effet, M£gO pur a une large bande interdite, 
cependant, lorsqu'il est dopé avec un élément de transition (TM), il conduit à une combinaison des 
propriétés de magnétisme, d'optique et de transport, ce qui en fait un seul matériau polyvalent 


pouvant être utilisé dans la fabrication d'appareils multifonctions [11]. 


Les progrès des recherches informatiques ont considérablement transformé l’aspect de la 
physique des matériaux en réduisant les contraintes de temps et de coûts pour les expérimentateurs 
scientifiques. L'utilisation de simulations informatiques a permis aux chercheurs de modéliser 
virtuellement les propriétés des matériaux à différentes échelles, offrant des prédictions précieuses 


avant même la réalisation d'expériences physiques et de guider les expérimentateurs vers des voies 
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plus prometteuses. Cette approche prédictive a ouvert de nouvelles perspectives en matière de 
conception de matériaux, d'optimisation des processus de fabrication, et de compréhension des 


comportements complexes à l'échelle atomique. 


Les algorithmes de recherche informatique ont également joué un rôle crucial dans la 
découverte de nouveaux matériaux aux propriétés spécifiques, offrant ainsi des solutions plus 
rapides et économiques. L'interaction synergique entre la modélisation informatique et les 
expérimentations a non seulement guidé les chercheurs vers des voies prometteuses, mais a 
également ouvert la porte à d'innombrables possibilités excitantes dans le domaine de la physique 


des matériaux, accélérant ainsi le rythme des découvertes et des innovations. 


Ces dernières années, et dans le cadre de la fonctionnalité, plusieurs modèles théoriques 
sont suggérés. La théorie de la fonctionnelle de la densité DFT (Density Functional Theory) [12] a 
émergé comme un modèle théorique crucial. Dans le présent travail, l'approche FP-LAPW (Full- 
Potential Linearized Augmented Plane Wave), basée sur la DFT est utilisée pour étudier les 
propriétés structurales, électroniques, optiques, et magnétique des composés dopés M£g;A:x0 (A 
= Ni, Pd, Pt, Cu, Ag, Au) et co-dopé EuMs,;B:1,0 (B = Cu, Ag, Au). Cette méthode offre une 
approche puissante pour comprendre le comportement des matériaux à l'échelle atomique, en 
utilisant la densité électronique comme variable clé. En se concentrant sur les propriétés 
structurales, électroniques et magnétiques, cette approche permet d'obtenir des informations 
détaillées sur la réponse des matériaux aux dopages. Les résultats de cette étude peuvent contribuer 
à éclairer la compréhension des propriétés fondamentales des composés dopés, ouvrant ainsi la 
voie à des applications potentielles dans divers domaines de la science des matériaux et de la 


technologie. 
Dans la présente thèse, notre approche s'organise autour des chapitres suivants : 


Chapitre 1 : Les formalismes de la méthode de la théorie de la fonctionnelle de la densité, ainsi 
que la méthode (FP-LAPW) sont exposés et expliqués. Cette base théorique solide établit le 


fondement nécessaire pour les investigations ultérieures. 


Chapitre 2 : Donne une revue de la littérature, mettant en lumière les travaux antérieurs sur les 


composés à base de M£gO, les domaines d’application.. et définissant les objectifs de cette thèse. 
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Chapitre 3, 4 et 5 : Présentent les résultats obtenus pour les propriétés structurales, électroniques, 
optiques, et magnétiques de MgO dopé et co-dopé. Ces résultats sont discutés et comparés aux 


résultats disponibles dans la littérature scientifique. 


Enfin, dans la dernière section dédiée à la conclusion, sont synthétisés les principaux résultats 
obtenus au cours de cette étude, soulignant les perspectives intéressantes ouvertes par ce travail, 
notamment dans l'étude de multiples propriétés en vue d'une utilisation potentielle de ces matériaux 


dans des applications technologiques avancées. 
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CHAPITRE I 


Cadre théorique et méthodes de calculs 


Dans ce chapitre, les principes physiques fondamentaux de ce domaine sont explorés. Les approximations appliquées pour 
simplifier et résoudre l'équation de Schrôdinger, ainsi que les théorèmes de Hohenberg-Kohn et Born-Oppenheimer sont examinés 


en détail. De plus, les méthodes des ondes planes augmentées linéarisées (FP-LAPW) ainsi que le code utilisé sont introduits. 


"Toute la science exacte est dominée par l'idée d'approximation " 


Bernard Russel 
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Partie A : Contexte théorique 


L.1. Introduction 


Suite aux succès récents dans la description et la prédiction des propriétés des matériaux, 
les simulations atomistiques en général et les calculs de structures électroniques en particulier ont 
acquis une importance considérable dans les domaines de la physique et de la chimie au cours de 
la dernière décennie, en grande partie grâce à l'émergence des technologies informatiques à haute 
performance [1]. En supposant que nous ayons connaissance des types d’atomes composant un 
matériau donné, une approche computationnelle nous permet de répondre à deux questions 


fondamentales : 


> Quelle est la structure atomique du matériau ? 
> Quelles sont ses propriétés électroniques ? 


De plus, ce serait bien d'avoir la réponse à une autre question : 


> Comment peut-on altérer la liaison entre les atomes ou la composition chimique d'un matériau 
pour engendrer de nouveaux matériaux dotés de propriétés spécifiquement prédéfinies ? 

Un certain nombre de méthodes sont développées pour trouver des réponses à ces 
questions. Ces méthodes de calculs de la structure et des propriétés des matériaux peuvent être 
conditionnellement divisées en deux classes : celles qui n'utilisent aucune quantité dérivée 
empiriquement ou expérimentalement et celles qui le font [2]. Les premières sont souvent appelées 
ab-initio ou méthodes de premiers principes tel que la théorie de la fonctionnelle de la densité 
(DFT), tandis que les secondes sont appelées empiriques ou semi-empiriques. Les méthodes ab- 
initio s'avèrent extrêmement bénéfiques pour prédire les propriétés de nouveaux matériaux et 
prévoir les tendances sur un large éventail. Toutefois, si de tels paramètres ne sont pas disponibles 
ou si un système présente des phénomènes inhabituels encore mal compris, il est souvent nécessaire 
de recourir à des calculs ab-initio. Ils exigent davantage de ressources informatiques, limitant ainsi 
le traitement à des unités cellulaires plus restreintes par rapport aux calculs basés sur les champs 
de force. La caractéristique principale des méthodes du premier principe (ab-initio) est qu’elles 


peuvent être mises en œuvre sans aucune information expérimentale préalables sur le système. 


La mécanique quantique est impliquée dans la physique de solide pour étudier les 


systèmes moléculaires et comprendre leurs propriétés physiques et chimiques. Erwin Schrôdinger 
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a développé dans les années vingt, une équation mathématique fondamentale pour décrire le 
comportement des électrons et des atomes. Cette équation constitue un pilier de la mécanique 
quantique et a permis de faire d'énormes avancées dans la compréhension de la structure atomique 
et moléculaire. L'équation de Schrôdinger est une équation aux dérivées partielles qui décrit 
l'évolution de la fonction d'onde d'un système quantique, telle que l'électron d'un atome ou d'une 
molécule. Cependant, résoudre l'équation de Schrôdinger pour des systèmes atomiques ou 
moléculaires plus complexes contenant de nombreux électrons interagissant entre eux n'est pas 


aussi simple que pour un électron unique. 

Extraire une solution de l'équation de Schrüdinger qui se rapproche du système physique 
réel est le but de la physique quantique [3]. 
L2. Théorie 


L'équation de Schrôdinger est développée par Erwin Schrôüdinger en 1926. Elle est au 
cœur des études quantiques sur les cristaux et les systèmes quantiques en général [4], c’est l'une 
des équations les plus importantes en physique. Elle décrit l'évolution de la fonction d'onde d’un 
système quantique dans le temps et permet de comprendre le comportement des systèmes 
quantiques, y compris les cristaux, Ainsi elle fournit une interprétation précise des phénomènes qui 
ne peuvent pas être expliqués de manière classique. Elle est décrite par l'équation suivante : 

HY=EY (1.1) 
Avec : 
H : Hamiltonien. 
Y : La fonction d'onde du système. 
E : Son énergie propre. 


L'Hamiltonien total du système est composé de l'énergie cinétique de toutes les particules, de 
l'énergie des interactions entre elles et dans certains cas de l'énergie interagissant avec le milieu 


extérieur en l'absence de champ extérieur, l'Hamiltonien s'écrit sous la forme suivante : 
H = Te + Ty + Ve + Ve_N + VN-N (L 2) 


Avec : 
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L'énergie cinétique des électrons : 


D eN = Ÿ (y L3 


L'énergie cinétique des noyaux : 


h2 
ps si rm or, nant 2 I. 4 
TN :2M; Es 2M; (Va ) ŒE 


Les interactions coulombiennes répulsives entre les électrons : 


— 1 >> e? LS 
€ Ans . 2 (É—5) ŒS) 


ij#i 


Les interactions coulombiennes attractives entre les électrons et les noyaux : 


ÿ —1 Y e? z,e (6) 

eN — = — L 
ATEo Æia(r, — R,) 

Les interactions coulombiennes répulsives entre les noyaux : 


LL e2 Z,Zy 
NN — Ane, 2 ne =. (L 7) 
TE a;B (R; — R;) 
Avec : i etj pour les électrons, a et b spécifique aux noyaux, M et m sont les masses des noyaux et 
des électrons, respectivement. Z,;e charge des noyaux, R, = R, la distance entre deux noyaux a 


et b, T; — R, la distance entre l’électron et le noyaux a, T; — T'; est la distance entre les deux 


électrons. 
Au repos, l'équation de Schrôdinger est indépendante du temps 

HY(rR)=EY(r;R) (L. 8) 
Où Y représente la fonction d'onde et se rapporte à toutes les coordonnées des particules, qu'il 


s'agisse d'électrons ou de noyaux, E est la valeur propre correspondante. 


La difficulté réside dans le fait que pour les systèmes contenant plusieurs électrons, les 
interactions électrostatiques entre les électrons deviennent significatives, ce qui rend l'équation de 


Schrôüdinger intraitable de manière exacte. Dans le but de rendre sa solution plus simple et plus 
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facile, plusieurs approximations et approches sont développées telles que l'approximation de Born- 


Oppenheimer, l'approximation de Hartree et l'approximation de Hartree-Fock. 


L3. Approximation de Born-Oppenheimer 


L'approximation dans laquelle on considère que les noyaux sont stationnaires par rapport 
aux électrons est connue sous le nom d'approximation du noyau stationnaire ou d'approximation 
du noyau immobile. Cette approximation est couramment utilisée en chimie quantique pour 
simplifier les calculs liés à la structure électronique des atomes et des molécules. Elle est basée sur 
l'idée que les électrons se déplacent beaucoup plus rapidement que les noyaux, ce qui permet de 
traiter les mouvements des électrons de manière indépendante de ceux des noyaux dans de 
nombreux cas. En raison du déplacement beaucoup plus rapide des électrons par rapport au noyaux 
ce qui permet de traiter les mouvements des électrons de manière indépendante de ceux des noyaux, 
cette approximation suppose que les noyau sont stationnaires par rapport aux électrons [5], c'est-à- 
dire que l’énergie cinétique est nulle Ty = 0, et l'énergie d'interaction entre les noyaux est 


constante Vyx , alors l'Hamiltonien électronique peut être défini de la manière suivante : 
EST + Ver Vr (L.9) 


L'équation de Schrüdinger devient alors : 


CEE VE + 2 Du Dia 1 Ve = Ee Ve (I. 10) 


1 ri—Ral 


Pour résoudre l'équation de Schrôüdinger, les méthodes de Hartree-Fock sont couramment 
utilisées en chimie quantique pour les atomes et les molécules, tandis que la théorie de la 
fonctionnelle de la densité (DFT) est une méthode polyvalente utilisée pour une variété de 


systèmes, y compris les solides. 
L.4. Approximation de Hartree et Hartree-Fock 


L'approximation de Hartree est basée sur le principe d'exclusion de Pauli : qui stipule qu’il 
n'est pas permis à deux électrons d'avoir le même ensemble de nombres quantiques. L'autre système 
de particules est résolu en supposant que les électrons sont indépendants les uns des autres c'est-à- 


dire que les interactions électron-électron sont négligées pour simplifier le problème, 
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L'Hamiltonien du système est alors traité comme celui d'un ensemble à N particules non 


interagissantes et la fonction d'onde des électrons est exprimée de la manière suivante [6,7]: 
y CA D, T3 sx IN) — ÿ, (Tr) (D) P:(r3) … PN() (L. 11) 
Dans cette équation, Vy(r y) est une fonction d'onde de N électrons. 


Cette approximation a contribué à transformer le problème de l'étude de la particule N en une étude 
d'un électron indépendant, où l'Hamiltonien est formulé comme une somme d'Hamiltonien partiels, 


dont chacun décrit l'état d'un électron indépendant comme suit : 
H — (Ts Æ Vext a Va)@) (L 12) 


Où T, est l'énergie cinétique, V,,+4 est le potentiel externe et V,, est le potentiel Hartree pour les 


électrons n'interagissant pas, et cela peut s'écrire comme : 


N 2 
_ _1 IE GE [ra drdr 
VW = —— ) ER ee | 


_ L 13 
8TE£) ( ) 


L5. Théorie de la fonctionnelle de la densité 


La DFT repose sur les travaux pionniers de Kohn et de Hohenberg. Ces théorèmes ont 
permis de réduire l'équation de Schrüdinger à 3N coordonnées spatiales aboutissant à une 
représentation en termes de la densité électronique, simplifiant ainsi le traitement des systèmes 
électroniques. Le premier théorème dit que la densité de l’état fondamental (n,) émerge de manière 
unique à partir d'une fonction énergie réduite avec succès. Tandis que le deuxième théorème établit 
que la théorie des fonctions d’ondes électroniques et détermine ainsi toutes les propriétés de l’état 
fondamental d’un système électronique et peut être entièrement déterminée par la densité 
électronique de l’état fondamental, révélant ainsi toutes les propriétés essentielles du système 
électronique [8]. La DFT représente une avancée significative qui remplace la fonction d'onde 
électronique par la densité électronique. Cette théorie est plus puissante et plus moderne que 
l'approximation de Hartree Fock. Elle permet de résoudre le problème à plusieurs corps en le 
reformulant en termes de problème à un seul corps, simplifiant ainsi considérablement les calculs 


pour de nombreux systèmes complexes. 
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L.5.1. Théorèmes de Hohenberg et Kohn 


En 1964, Hohenberg et Kohn [9, 10] ont établi les bases de la théorie de la fonctionnelle 
de la densité. Leurs travaux ont instauré les fondements de cette approche puissante qui a 
révolutionné la compréhension des systèmes électroniques, en se basant sur des principes 


fondamentaux. 


Dans la formulation de la DET. Les noyaux des molécules ou des clusters, traités comme 
d’habitude dans les calculs de structure électronique à N corps, sont considères comme fixes 
(approximation de Born-Oppenheimer), générant un potentiel externe statique V dans lequel les 


électrons se déplacent. L’état électronique stationnaire est alors représenté par une fonction d’onde 


W (Ta ,T2,T3. Ty) satisfaisant l'équation de Schrüdinger à plusieurs électrons : 


fw = [f+9+ 0] = [EY-* V2 + XNV(R) + ENU (ff) | V=ET 
(I. 14) 


Où Hest l’hamiltonien moléculaire électronique, N est le nombre d’électrons, T est l’énergie 
cinétique des N électrons provenant du champ externe, Ü est l’énergie d’nteraction electron 
electron pour le système électronique. Les opérateurs T et Ü sont dits operateurs universels car ils 
sont identiques pour tout systèmes, tandis qu’ils sont dépendants du système, c’est-à-dire non 
universels. La différenciation entre les problèmes séparés à une seule particule et le problème 


beaucoup plus complexe à plusieurs particules provient du terme d’interaction Ü. 


Les équations de Schrüdinger pour les systèmes à plusieurs corps peuvent être résolues 
par de nombreuses techniques complexes basées sur l’extension de la fonction d’onde dans les 
déterminants de Slater. Bien que la méthode Hartree-Fock soit la plus simple, les techniques post- 
Hartree Fock offrent des approches plus sophistiquées pour améliorer la précision de la description 
de la structure électronique. Ces méthodes, incluant également la théorie de la fonctionnelle de la 


densité et la théorie de perturbation à plusieurs corps [11]. 


Cependant, ces techniques rencontrent un défi majeur lié à l’énorme coût computationnel, 
qui rend pratiquement impossible leur application efficace à des systèmes plus grands et plus 
complexes. La DET offre un substitut attrayant, étant beaucoup plus adaptable car elle fournit un 
moyen de cartographier systématiquement le problème a plusieurs corps, avec Ü, vers un problème 


à un seul corps sans Ü. Cette approche permet ainsi de simplifier le traitement des systèmes 
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complexes tout en préservant une précision adéquate La variable clé est la densité des 


particules n(r), qui, normalisée, est donné par : 
n(F) = N [ dé rl J d3 PV … J d° IN * (PP, EN )Y (PP, EN ) (I. 15) 


Cette relation peut être inversée, c’est-à-dire pour une densité d’état fondamental donnée ns (F), 
il est principalement possible d’élaborer la fonction d’onde équivalente à l’état fondamental 


Y (ra T2 T3 28 Ty). 


Autrement dit, W, est une seule fonctionnelle de no, Vo = Ÿ [no] et par conséquent, la valeur 


attendue de l’état fondamental d’une observable O est également une fonctionnelle de no 
O [nol = < Y [no] [0] Y [no] > (I. 16) 
En particulier, l’énergie de l’état fondamental est une fonctionnelle de ns. 
Es = E no] 
= < Y [no] [T + V + U| Y [no] > 
= < [no] [T| Y [no] > + < Y [no] |P] Y [no] > + < Y [no] [0] Y [no] > (1.17) 
Où < Y[no] LA Y [no] > représente la contribution de l’extérieure 


Le potentiel peut être écrit clairement en termes de densité de l’état fondamental no 
V [no] = Jv () no À) dér (I. 18) 


Plus couramment, la contribution du potentiel externe < Y LA Y >, peut s’écrire clairement en 


termes de densité n : 


V [no] = Î V (F) no) dr (1.19) 


Comme mentionné ci-dessus, T [n| et U [n] sont appelés fonctionnelles universels, alors que V [n] 
est appelée fonctionnelle non universelle, car elle dépend du système étudié. Ayant un système 


défini, c’est-à-dire ayant spécifié Ÿ, il faut alors réduire les fonctionnelles : 


E(n)=T{[n]+U[n]+fVPnPdr (I. 20) 
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En considérant n (r)et en supposant des conditions fiables pour T [n] et U [n]|, une minimisation 
réussie de la fonctionnelle énergétique conduira à la détermination de la densité de l’état 


fondamental n, et donc à toutes les autres valeurs des observables de l’état fondamental. 


La technique lagrangienne des multiplicateurs indéterminés peut être appliquée pour 


résoudre les problèmes de variation de minimisation de la fonctionnelle énergétique E (n) 
Es [n] = (#[n]|fs + [W,[n]) (21) 


OùT, indique l’énergie cinétique sans interaction et Ÿ, est un potentiel efficace externe dans lequel 


les particules se déplacent. Clairement, si Ÿ, est n, (7) = n (Fr), sélectionné pour être : 
V.= V+Ü+(T-T.) (EL 22) 


Par conséquent, il devient possible de résoudre les équations de Kohn-Sham de ce système non 


interactif. 
Év2 + M0: M = € 0) (1.23) 


Les équations de Kohn-Sham conduisent à la détermination des orbitales @; qui reproduisent la 


densité n (Tr) du systéme a plusieurs corps d'origine. 


N 
n = n, M = D 16: À (1.24) 
i 
Le potentiel effectif d’une particule unique peut s’écrire plus en détail sous la forme : 
—> — => e2n,(r) 3 où ne 
V,) =V E)+ SJ F4 F + Vxc n,(r)] (I. 25) 


Où le deuxième terme représente ce que l’on appelle l’expression de Hartree décrivant la répulsion 
coulombienne électron-électron, tandis que la dernière expression Vyc est appélée le potentiel 
d’échange corrélation. Ici, Vyc englobe toute les interactions entre plusieurs particules, puisque 
l'expression Hartree et Vyc dépendent de n (Tr) qui à son tour dépend des orbitales @; et de 
l'opérateur Ÿ, correspondant. Pour résoudre de manière cohérente le problème de la résolution de 
l'équation de Kohn-Sham, il est recommandé de commencer par une première estimation de n (Tr), 


puis de calculer l’équivalent et de résoudre les équations de Kohn —-Sham pour les @;. Ce processus 
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itératif permet d'obtenir une description auto-cohérente de la densité électronique et des orbitales 


pour le système considéré. 


L5.2. Equations de Kohn et Sham 


Bien que les théorèmes de Hohenberg-Kohn soient puissants, ils ne peuvent pas offrir une 
méthode appropriée pour calculer la densité de l'état fondamental. Kohn-Sham [12] ont introduit 


une approche pour effectuer des calculs DFT. 


Ils se sont concentrés sur la transformation du système en interaction avec un potentiel 
réel vers un système fictif sans interaction, où les électrons se déplacent dans un potentiel de Kohn- 
Sham effectif Ve, où Vxs. Cette méthode demeure exacte car elle donne la même densité d’état 
fondamental que le système réel en supposant que les électrons non interactifs ont la même densité 


que les électrons réels. 


Dans cette approximation la fonctionnelle F [n(r)] est composée en trois termes : les deux 
premiers étant connus et représentent l'essentiel de l'énergie, tandis que le troisième terme demeure 


une inconnue, F [n(r)] est donné par : 


Fn()]= Tin(n)] + Î Vext 7) n () dr + Egin(n)] + Excin(n)] (1.26) 


Où T,[n(r)] est l'énergie cinétique du gaz électronique sans interaction de densité n(r). 


E[n(r)| est l'énergie classique de Hartree des électrons donnée par : 


1 fn(r)n(r) 


Ex —= 2 Es dr dr (L 27) 


et Exc[n(r)] est l'énergie d'échange-corrélation. Le potentiel effectif Ve est maintenant donné 
par : 


0Exc [n] 


Verr QD = Vext G) + EG + 


(L. 28) 


Nous savons que les électrons sont des fermions, ce qui indique que la fonction d'onde de Kohn- 
Sham est obtenue à partir du déterminant de l'orbitale d'une seule particule, conduisant ainsi à la 


solution de plus basse énergie de l'équation de Schrüdinger à N un électron : 
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h2 
Ve iCr) + Vers (r)| Dir) = e; Di(r) (29) 


Cette équation aux valeurs propres est une représentation particulière des équations de Kohn-Sham. 
La densité électronique n (r) peut être représentée comme la somme des carrés d’un ensemble de 


fonctions d’onde orthonormées ; (r) : 


N 
n (= D I P (1.30) 


Après avoir résolu l'équation de Schrüdinger pour ce système, la densité résultante correspond à la 
densité de l'état fondamental, dénotée n (Tr) ), comme définie dans l'équation (I. 30). Si n (r) est 


calculé alors : 
Ves = Y[n(r)] (1. 31) 


Si nous avons Y &$, nous pouvons alors obtenir absolument la densité [n(r)] de l'état fondamental. 
Cela signifie que la fonction d’onde W,4 est une fonction unique de la densité de l’état 


fondamental. 


Ces trois équations forment les équations orbitales de Kohn-Sham dans leur forme 
standard. Ce système est ensuite résolu de manière itérative, jusqu’à ce que l’auto-cohérence soit 
approchée. Il est essentiel de souligner que les valeurs propres &; n’ont aucune signification 
physique, seulement la somme totale, qui correspond à l’énergie de l’ensemble du système E à 


travers l’équation : 


N 


ÔE 
E = 7 €; — Valn] + Exc [n] - EE n(r)dr (I 32) 


i 
La figure L1 est une représentation schématique de la boucle auto-cohérente pour la 
solution des équations de Kohn-Sham. D'une maniérée générale, il est nécessaire de répéter 
simultanément deux de ces cycles, chacun correspondant à l'un des deux spins, le potentiel de 
chaque spin dépendant de la densité des deux spins. En pratique, différentes approches de la théorie 
de Kohn- Sham peuvent être appliquées en fonction du contexte spécifique des calculs à l’état 


solide. 
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Créer des fichiers de base en Générer les densités atomiques et 
utilisant les entrées initiales analyser le traitement des orbitales 


Produire une densité pour SCF par 
superposition des n atomiques 


Calculer les valeurs et les fonctions Calculer les états du noyau et la 
propres densité électronique du noyau 
[Ts + VIF x = Ex Tr HF ut = Ent Tnt 


Calculer les densités de valence à partir 
de vecteurs propres 


Nyat = Fr nn Mélange Les densités 
wat — K°K “| d'entrée (anciennes) et de 
sortie 


Figure L.1 : Schéma du cycle auto-cohérent utilisé pour la résolution des équations de Kohn- 


Sham. 


L.6. Différents types de fonctionnelles 


Comme précédemment indiqué, tous les termes énergétiques et leurs potentiels associés 
peuvent être évalués, à l’exception de l’échange et de la corrélation, qui constituent le principal 


défi de la DFT. Afin de l’évaluer, une approximation doit être utilisée. 


L.6.1. Approximation de la densité locale 


L’approche la plus simple des fonctionnelles d’échange corrélation est l’approximation 
de la densité locale (LDA). Cette méthode, proposée par Kohn et Sham, s’est avérée 
remarquablement précise utile et difficile à améliorer. L'approximation LDA est largement utilisée 
pour décrire avec précision diverses interactions d’échange-corrélation à courte portée par exemple 


dans des systèmes de liaison covalente. Toutefois, 1l est important de noter que cette approximation 
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ne peut pas fournir une estimation précise de l’interaction d’échange corrélation à longue distance, 


comme celle caractérisant les forces de Van der Waals (VdW) [13]. 


L'approximation LDA est locale, ce qui signifie que l'énergie d'échange et de corrélation 
électronique en tout point de l'espace est fonction uniquement de la densité électronique en ce point 
précis. La fonctionnelle XC est la somme d'une fonctionnelle de corrélation et d'une fonctionnelle 


d'échange : 
Exc — Ex + Ec (L 33) 


Cette approche consiste à appliquer les résultats exacts du gaz uniforme d’électrons aux systèmes 


réels, particulièrement où l'échange doit être évalué : 
= Î dérn (| (3r2n(%) (1. 34) 


n(r) est la densité électronique par unité de volume au point (r), e est la charge d’un électron. 


En cherchant les issues au système d’équation de Schrôdinger : 
1 
E Y(x) = F9 V= Y(x) + V(x) (x) (135) 


On constate que tous les termes sont exprimés en tant que fonctionnels de la densité de charge 
électronique. Ce qui rend ce système plus facile à résoudre ou, plus précisément, nécessite moins 
de calculs que, par exemple, les équations de Hartree-Fock, est que la possibilité efficace est locale. 
Ainsi, la résolution de l’équation de Schrôdinger n’est pas plus complexe que l’approximation de 
Hartree, à condition que l’énergie d’échange-corrélation puisse être représentée en fonction de la 
densité de charge locale. L’énergie d'échange-corrélation d'un système électronique est construite 
en supposant que l'énergie d'échange-corrélation pour chaque électron en un point du gaz 
électronique est égale à celle pour chaque électron dans un gaz électronique identique ayant la 
même densité électronique au même point. L'énergie d’échange-corrélation s’exprime alors selon 


l’équation suivante : 


EXPA[n G] = J n Q ete [n Q)] d5 à (1.36) 
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grom In ©] : L'énergie d’échange et de corrélation d’un gaz homogène d’électrons de densité 


constante est égale à sa valeur en F. 


L.6.2. Approximation du gradient généralisé 


L'étape suivante inévitable dans l'amélioration des formules approchées pour les 
fonctionnelles d'échange-corrélation est l'expansion de Eyç en termes de gradients de la densité. 
La méthode la plus prometteuse pour améliorer la LDA consiste à introduire l’inhomogénéité dans 
l’énergie d’échange-corrélation, en incluant des termes dépendant du gradient de la densité 
électronique [14]. Selon l'approche GGA, cette énergie peut être exprimée de différentes manières 
par la combinaison des termes locaux et des termes dépendant du gradient. Elle se présente sous la 


forme suivante : 
ESA In G)] = f (n@)ei" [n @]Fxe, V [p @), Vn )Jdr (L. 37) 


L'approximation du gradient généralisé (GGA) a attiré beaucoup d'attention pour sa 
simplicité conceptuelle et de la charge de calcul modérée. De nombreux calculs évaluant la 
précision de la GGA ont été rapportés et démontrent communément que la GGA corrige 
considérablement l'erreur LDA dans les énergies de cohésion des molécules et des solides [15]. 
Des études approfondies de référence, menées par Peter BLAHA et al. ont examiné diverses 
fonctionnelles. Les résultats de ces travaux ont montré que parmi les fonctionnelles GGA [16, 17], 
celles avec un petit facteur d'amélioration telles que AMOS [18], PBEsol ou quelques autres sont 
les plus précises pour la constante de réseau et le module de masse, tandis que le PBE standard est 
le meilleur choix pour l'énergie de cohésion [19]. Au niveau MGGA, la fonctionnalité SCAN 
devient de plus en plus populaire et s'est avérée simultanément aussi bonne que les meilleurs GGA 
pour la géométrie (par exemple, PBEsol) et l'énergie de cohésion (PBE) [20]. Cependant, il 
convient de mentionner que SCAN peut poser des problèmes pour les systèmes magnétiques itérant 


ou métaux alcalins [21]. 


L.6.3. Potentiel de Becke Johnson modifié par Tran_Blaha 


Le potentiel de Becke Johnson modifié par Tran_Blaha (TB-mBJ) se compose d'une 
version modifiée du potentiel d'échange BJ, qui reproduit de manière très précise le potentiel KS 


exact des atomes et de la LDA pour la corrélation. Le potentiel mBJ est un outil efficace 
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implémenté dans le code Wien2Kk contribuant à améliorer la structure de bande des matériaux, en 
particulier pour les semi-conducteurs. Son utilisation s'est avérée cruciale pour parvenir à un accord 
plus précis avec les observations expérimentales, élargissant ainsi la portée de la méthode, comme 
l'ont démontré Tran et Blaha dans des études antérieures portant sur un ensemble de semi- 


conducteurs [22]. 


La formule ci-dessous représente ce potentiel : 


Vs) = cvéR(r) + (8c— 2)= : D 


(L 38) 


Où : vËR est le potentiel de Beck Roussel (BR) qui a été proposé pour modéliser le potentiel 
coulombien créé par le trou d'échange et donné par la relation de la référence [23, 24]. C a été 


choisi pour dépendre linéairement de la racine carrée de la moyenne de |V, |/n : 


C= a+ Bg? avec g = (— mo) d°r')1/2 (L. 39) 


Vcell n(r') 
Où a et B sont deux paramètres libres et Ven est le volume unitaire de la cellule. 


Où par défaut : « = —0.012 et f = 1.023 Bohr. 
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Partie B : Méthodologie 


L.1. Introduction 


Cette section commence par une introduction aux concepts de la méthode d'onde plane 
augmentée (APW) et de méthode des ondes planes augmentées plus les orbitales locales 
(LAPW+#lo) utilisées pour résoudre l'équation de Kohn-Sham dans le contexte d’un solide 
périodique [25]. Ensuite, les concepts de la méthode du plein potentiel (FP) et de l'approximation 
du Muffin Tin avec le plein potentiel sont décrits. Enfin, le programme Wien2K et certaines de ses 


applications sont présentés. 


Il existe plusieurs techniques pour résoudre les fonctions à un électron, la base utilisée 
dans ces calculs repose sur les ondes planes (PW), qui correspondent aux fonctions de Bloch. 
Cependant, afin de rendre cette base plus efficace, il est nécessaire de les augmenter car elles se 
révèlent inefficaces pour décrire les fonctions d'onde qui varient rapidement à proximité des 


noyaux. 


L.2. Méthode des ondes planes augmentées 


L'origine de la méthode des ondes planes augmentées sont conçue en 1937 par Slater [26], 
lorsqu'il a introduit l'approximation dite du "muffin-tin" pour construire un ensemble de fonctions 
de base. Par conséquent, cette méthode est alors appelée onde plane augmentée. La méthode des 
ondes planes augmentées réside dans l'utilisation d'ondes planes, dans les régions où le potentiel 
varie lentement, complétées par des fonctions de type orbital atomique dans les régions où les 


fonctions d'onde varient plus rapidement. 


L3. Fonctions de base dans le potentiel Muffin-Tin 


Bien que le potentiel dans le cristal varie de manière continue, dans l'approximation du 
muffin-tin l'espace cristallin est divisé en deux régions distinctes, comme illustré dans la figure 


L2. 


La première région est appelée "muffin-tin" (MT), qui contient des sphères qui ne se 
chevauchent pas, tandis que l'autre région est le reste de l'espace cristallin, appelé région 


interstitielle. Cette division du cristal conduit à une séparation du potentiel en deux parties. Le 
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premier potentiel est un potentiel symétrique presque sphérique, situé dans la région MT ; tandis 


que le deuxième potentiel est un potentiel plat, situé dans la région interstitielle. 


Cette approche repose sur l’hypothèse selon laquelle la fonction d’onde et le potentiel à proximité 
des noyaux sont similaires à ceux d’un atome dans la région MT, mais dans la région interstitielle, 


ils sont plus caractérisés par une plus grande régularité. 


A l'intérieur de la sphère atomique i de rayon Ri (proche des noyaux), les électrons se 
comportent comme s'ils étaient dans un atome libre. Dans cette région, la fonction radiale est 
utilisée. En revanche, dans la région interstitielle restante, (loin des noyaux) où les électrons sont 


libres, des ondes planes sont utilisées. 


Selon cette méthode, la fonction d’onde est donnée par : 


DE TO TORRES 


œ(r) = ne _—— (I. 40) 
a1/2 > Cc exp(i(K + G).f) ,r > Ra 
G 


Où Q est le volume unitaire de la cellule, Tr est la position à l’intérieur de la sphère a, K est le 
vecteur d’onde dans la zone de Brillouin irréductible. Y,,, est la fonction harmonique sphérique du 
moment cinétique des nombres quantiques Let m, et u%, est la solution numérique qui satisfait 


l’équation radiale de Schrôdinger [27] : 


d (+1) 
Re ue À V() — Ej|ruim (7) = 0 (41) 


Où V est la composante sphérique du potentiel et pour une énergie fixe Ey, u, (r) doit être défini 


ES 


à r =0. 


Pour garantir la continuité du potentiel aux frontière des sphères, la contrainte se fait en définissant 


Um en termes de coefficients d’onde plane C4; comme : 


Ami! — 
Um = 7 ——— > Co Je UK + 8D)Yin Ge + 6) (1.42) 


Q72 u (Ra) G 


Où CQ sont les coefficients des ondes planes. 
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Région Interstitielle I 


Sphères (MT) 


Région (MT) 
Il 
Région (MT) 


Région Interstitielle I 


Figure I.2. Division d'une cellule unitaire en régions des sphères (MT) et région interstitielle, 


pour un cas avec deux atomes. 


L.4. Méthode des ondes planes augmentées plus les orbitales locales 


Une approche alternative est proposée par Sjüstedt et al. [28], à savoir la méthode APW+lo 
(orbitales locales). Dans cette approche, l'augmentation est similaire au schéma original APW, mais 
avec la particularité que chaque fonction d'onde radiale est calculée à une énergie de linéarisation 
fixe. Cette modification vise à éviter le problème lié aux valeurs propres non linéaires rencontrées 
dans la méthode APW. Ainsi, pour être valide, la condition de continuité doit être respectée et les 
fonctions de base présentent une coupure à la limite de la sphère. Les orbitales locales peuvent être 


définies comme : 


p®) = ile ro + DE u(F)]Y im) r<Re (1.43) 


O, reEl 
Les orbitales locales sont évaluées à une énergie fixe identique à celle des APW 
correspondants. Les deux coefficients sont déterminés par le processus de normalisation et la 
condition selon laquelle l'orbitale locale a une valeur nulle à la limite de la sphère. Dans cette 
version, à est indépendant des ondes planes, Récemment, il a été démontré que cette nouvelle 
approche peut atteindre la même précision que LAPW mais converge plus rapidement en termes 


du nombre d'ondes planes nécessaires [28]. 
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L'efficacité optimale a été trouvée pour un ensemble de bases mixtes dans lequel les 
nombres quantiques "physiquement importants" / sont traités par la méthode APW+#lo, tandis que 
les nombres quantiques / plus élevés sont traités par la méthode LAPW. Selon la référence [28] il 
a été démontré que des quantités telles que l'énergie totale et les forces convergent de manière 
significative plus rapidement en termes dénombre de fonctions de base par rapport à la procédure 
LAPW pure. Cependant, malgré cette convergence plus rapide, les valeurs atteintes sont 
équivalentes à celles obtenues avec la procédure LAPW classique Dans LAPW, lorsque la base est 
trop petite, la force change de signe et l'atome se déplacerait dans la mauvaise direction. En 
revanche, dans la méthode APW+lo, la force converge en douceur et beaucoup plus rapidement. 
Pour les systèmes de grande envergure, il est possible de réduire la taille de la matrice N d'environ 
la moitié, entraînant une diminution substantielle du coût de calcul, étant donné que la 
diagonalisation évolue avec N°. Ainsi, la méthode APW+lo offre une alternative plus efficace, 


permettant de réduire considérablement les coûts de calcul pour les grands systèmes. 


Le nouveau système combine les meilleures fonctionnalités de toutes les méthodes basées 
sur l’APW disponibles. La LAPW converge plus lentement que la méthode APW, comme l'ont 
déjà souligné Koelling et Arbman [29], Cette différence de convergence est attribuée à la contrainte 
imposée par la LAPW, qui nécessite des fonctions de base différentiables et rend la LAPW moins 


adaptés de manière optimale pour décrire les orbitales à l'intérieur de la sphère. 


Cela justifie le retour à l'APW, mais la base indépendante de l'énergie, introduite dans 
LAPW, demeure cruciale pour éviter le problème des valeurs propres non linéaires et est donc 


également conservée. Les orbitales locales offrent la flexibilité nécessaire [30]. 


L5. Méthode des ondes planes augmentées linéarisées à potentiel total 


L'un des schémas les plus précis pour résoudre les équations de Kohn-Sham est la méthode 
d'onde plane augmentée linéarisée à plein potentiel (FP-LAPW) suggérée par Andersen [31], sur 
laquelle est basé le code Wien. La méthode FP-LAPW combine le choix de l'ensemble de base 
LAPW avec le traitement du plein potentiel et de la densité de charge, sans aucune approximation 


de forme dans la région interstitielle et à l'intérieur des sphères muffins-tin. 


Cette généralisation est obtenue en relâchant le potentiel interstitiel constant Vr et 


l'approximation muffin-tin sphérique Vur(r) grâce à l'inclusion d'un potentiel interstitiel déformé 
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Z VÉe KT, où K sont tous les vecteurs du réseau réciproques jusqu'à la plus grande valeur de Kmax. 


Les termes non sphériques à l’intérieure des sphères Muffin-Tin sont les suivantes : 


7. Vim Yim(®) Muffin — Tin 


V(F) = (1.44) 


Im 
D Veiki Interstitielle 
K 


L'approximation Muffin-Tin, bien qu’assez précise tel que dans les systèmes métalliques, 
a été largement utilisée dans les années 1970. Néanmoins, cette approximation s'avère inadéquate 
et conduit à des incohérences significatives avec les observations expérimentales dans le cas de 


solides covalents ou, les systèmes en couches. 


Dans ces conditions un traitement sans aucune approximation de forme est indispensable. 
Le potentiel et la densité de charge sont développés en harmoniques à l'intérieur de chaque sphère 


atomique et en séries de Fourier dans la région interstitielle. 


L.6. Code Wien2k 


Tous les résultats présentés dans ce travail sont effectués avec le package de programme 
Wien2k, un code amélioré par Blaha, Schwarz et ses collaborateurs. Le programme Wien2k permet 
d'appliquer la théorie de la fonctionnelle de la densité (DFT) pour calculer les structures 
électroniques et magnétiques des solides. Il est basé sur la méthode LAPW+lo, qui est reconnue 


comme l'une des méthodes les plus précises pour les calculs de structure de bande. 


Dans la version la plus récente de Wien2K, un ensemble de base alternatif APW+lo, est 
employé à l'intérieur des sphères atomiques pour les orbitales importantes (ondes partielles) 
présentant des difficultés de convergence, telles que les états p, d ou f de valence la plus externe, 
ou pour les atomes nécessitant l'utilisation de petites sphères atomiques. Pour toutes les autres 
ondes partielles, la méthode LAPW est utilisée. De plus, de nouveaux algorithmes intensifs en 


informatique, pour résoudre le problème général des valeurs propres, sont mis en œuvre. 


Le programme est principalement divisé en deux parties, l'initialisation et le programme 
principal, également appelé cycle auto-cohérent. Chacune de ces parties est son tour subdivisé en 


sous-programmes, reliés par des scripts. 
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Les avancées notables dans les possibilités de simulation de systèmes sur des 
configurations informatique modérées se sont indéniables grâce aux progrès de l’algorithme et de 
la puissance de calcul. Le développement du nouveau package Wien2k [32] est développé en tenant 
compte de multiples facteurs fondamentaux pour la technique informatique moderne. L'objectif est 
d'atteindre une efficacité maximale et des performances exceptionnelles. L'introduction d'une 
nouvelle base mixte APW+ lo/LAPW, avec une taille de matrice réduite, permet de gagner du 


temps de calcul sur l'ordinateur et permet d'étudier des systèmes plus complexes. 


Wien2k est un schéma entièrement électronique, offre une approche avancée en incluant 
des effets relativistes et possède de nombreuses fonctionnalités. La richesse des fonctionnalités de 
Wien2k avec la combinaison d'approximations de la DFT et d'une approche entièrement 
électronique avec prise en compte des effets relativistes renforce la capacité du code à fournir des 
résultats précis et détaillés ainsi que d’explorer divers aspects des propriétés électroniques des 


matériaux, 


L7. Cycle auto-cohérent 


Le package Wien2k se compose de plusieurs programmes indépendants, d'initialisations 
et de programmes de cycle auto-cohérents. Les programmes d'initialisation déterminent le groupe 
spatial de la structure, gênèrent des densités atomiques libres, spécifient le traitement des 
différentes orbitales dans les calculs de structure de bande, élaborent un réseau k-mesh dans la 
partie irréductible de la zone de Brillouin et créent une densité de départ pour le SCF cycle par 
superposition des densités atomiques préalablement générées pour les atomes libres. Ces 
programmes initiaux établissent ainsi les conditions de départ essentielles pour le reste du 
processus de simulation, contribuant à la précision et à la robustesse des calculs réalisés avec 
Wien2k. Dans le cycle SCF, le potentiel est construit à partir de la densité de charge laquelle est 
utilisée pour déterminer l'Hamiltonien servant au calcul des valeurs propres et des vecteurs propres 
des électrons de valence (bandes de valence). Par la suite, une nouvelle densité électronique de 
valence émerge des vecteurs propres calculés. Ensuite, les états centraux et les densités sont 
déterminés par un calcul auto-cohérent entièrement relativiste. Ces deux densités de charge 
(valence et noyau) sont mélangées à la densité d'entrée pour produire une entrée raffinée pour 
l'itération suivante du cycle. Un schéma de mélange le plus simple est représenté par le mélange 


direct, un processus crucial dans le maintien de la cohérence et de la convergence du calcul : 
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nue a ( = Qninput + Qnut (XL. 45) 


Où i fait référence au numéro d'itération et Q est le facteur de mélange. Le processus itératif est 
répété jusqu'à ce que la convergence soit atteinte. Cette convergence est déterminée en comparant 
les anciennes et nouvelles densités électroniques, et le cycle est répété tant que ces densités se 
situent dans les limites de tolérance définies par les critères de convergence du programme. Ce 
processus itératif de mise à jour des densités, de construction du potentiel, et de résolution auto- 
cohérente des équations se répète jusqu'à ce que le système atteigne un état stable où les densités 


électroniques convergent vers des valeurs constantes. 


Lors de la phase d'initialisation, la mise en place de la structure géométrique et 
électronique du système d'intérêt est réalisée. La symétrie du système et les k points qui seront 
utilisés sont déterminés. Lors de l'initialisation, les densités électroniques des atomes libres sont 
calculées et superposées pour obtenir une première estimation de la densité de charge qui servira 


de point de départ dans le cycle auto-cohérent. 


L.8. Applications de Wien2k 
L8.1. Systèmes 


Les problèmes abordés dans les calculs de mécanique quantique en utilisant la méthode 
LAPW telle qu'elle est implémentée dans différentes versions du code Wien2k comprennent une 
gamme variée de matériaux, incluant notamment les isolants, les semi-conducteurs, les métaux de 
transition ainsi que des systèmes à électrons f ou les composés intermétalliques. La structure de 
bandes peut être directement comparée à des expériences dans des cas faiblement corrélés. 
Cependant, la structure électronique de systèmes hautement corrélés tels que les supraconducteurs 
à haute Tc ou les oxydes de métaux de transition, nécessite des traitements allant au-delà de la LD A 
ou la GGA. Dans certains solides, le magnétisme joue un rôle important et tant que les moments 
magnétiques sont ordonnés selon un arrangement colinéaire (par exemple, ferro, ferri ou 
antiferromagnétiques), une description appropriée est relativement facile. Récemment, des 
extensions aux systèmes magnétiques non colinéaires tels que les spins inclinés ou les spirales de 


spin sont proposées [21-33]. 
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L8.2. Volume à l'équilibre et module de compressibilité 


Des séries de calculs d'énergie totale en fonction du volume sont adaptées à une équation 
d'états selon l'équation de Murnaghan : 


: Bo(V) | Vo Vo. 0 


- où FE, est l'énergie minimale, v, est le volume à l'équilibre, B et B' sont le module de 


compressibilité et sa dérivée. 
Avec: B= —-V — (1.47) 


- Le volume à l'équilibre et le module de compressibilité peuvent être déterminés à partir de 


l'équation ci-dessus en traçant l'énergie totale en fonction du volume. 
L8.3. Optimisation des forces et de la structure 


L'optimisation de la structure, souvent nécessaire dans ce contexte, est étroitement liée à 
l'énergie totale du système. Dans les cas où les positions atomiques ne sont pas fortement 
contraintes par la symétrie cristalline. La connaissance des forces agissant sur les atomes permet 
d'optimiser les paramètres de structure. Les forces calculées sont cruciales pour de telles 


optimisations. 


1.8.4. Structure de bandes et densité d'états 


La structure de bande électronique, également connue sous le nom de structure de bande 
d'un solide, en physique du solide, spécifie les plages d'énergie auxquelles un électron est autorisé 
ou interdit d’occuper. Elle est déterminée par les interactions complexes entre les électrons et le 
réseau cristallin. La connaissance de la structure de bande est déterminante pour comprendre 
plusieurs caractéristiques d'un matériau, notamment ses propriétés électroniques et optiques. La 
densité d'états (DOS) est une mesure complémentaire à la structure de bande. Elle indique le 
nombre d'états électroniques disponibles à chaque niveau d'énergie. Une DOS élevée à un niveau 
d'énergie donné indique un grand nombre d'états occupés, tandis qu'une DOS égale à zéro signifie 
qu'aucun état ne peut être occupé à ce niveau d’énergie. La combinaison de la structure de bandes 
d'énergie et de la densité d'états correspondante permet de caractériser de manière exhaustive la 


structure électronique d'un système. 


44 


L'analyse approfondie de la structure de bande et de la DOS fournit des informations 
cruciales sur les propriétés électriques d'un matériau, en identifiant s'il s'agit d'un métal, d'un isolant 
ou d’un semi-conducteur et donne un aperçu détaillé de la liaison chimique et des interactions 
électroniques locales. L’examen les densités partielles d'états décomposées sur site, permet de 
comprendre la force des interactions entre les orbitales des atomes constituants [26]. Une 
cartographie tridimensionnelle de bande par diffraction et photoémission à très faible énergie 
dépendant de l'angle offre une exploration détaillée de la dispersion des bandes électroniques, et 
fournit des informations pour comprendre les propriétés électroniques à un niveau atomique. Ces 
approches jouent un rôle dans la caractérisation et la conception de nouveaux matériaux aux 


propriétés électroniques souhaitées. 


L8.5. Energie totale et transitions de phases 


En calculant l'énergie totale, la stabilité relative des différentes phases peut être évaluée. 
Dans un tel cas, 1l est conseillé de garder autant de paramètres constants que possible afin de 
minimiser les erreurs systématiques. Ces paramètres peuvent inclure la taille des sphères 
atomiques, le maillage K, la fonctionnelle DFT [34], le traitement de la relativité. etc. Étant donné 
que les différences d'énergie sont souvent assez faibles, un traitement cohérent des systèmes à 


comparer contribuera à minimiser les erreurs informatiques [35]. 


L.9. Programmes utilisés dans le code Wien2K 


Comme évoqué précédemment, Wien2k est constitué de nombreux programmes 
individuels, chacun possédant son propre fichier d'entrée. Bien que cela semble très fastidieux au 
début, il existe des entrées par défaut pour tous les programmes et plusieurs outils permettant de 
modifier les paramètres les plus importants de manière dynamique [36]. Dans l'interface w2web, 
le prochain processus consiste à générer les fichiers d'entrée pour le calcul SCF et à vérifier la 
symétrie de la structure nouvellement créée (initialize @ w2web). L'utilisateur peut fournir quelques 
paramètres (nécessaires uniquement si l'on souhaite modifier les valeurs par défaut, voir ci- 


dessous) et exécuter les étapes suivantes en mode batch ou étape par étape : 


NN : trouve les distances entre chaque paire d'atomes jusqu'à deux fois la distance du voisin le 


plus proche. De plus, il vérifie les sphères qui se chevauchent et en cas de chevauchement, émet 
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un message d’erreur. Il vérifie également si des éléments identiques ont le même environnement 
et éventuellement les regroupe en ensembles équivalents. 

SGROUP : identifie le groupe d'espace en examinant la structure. Il regroupe les atomes en 
ensembles d'atomes équivalents en fonction des positions de Wyckoff du groupe spatial 
correspondant. En outre, 1l calcule la plus petite cellule (primitive) possible et si nécessaire 
généré le fichier de structure associé. Par exemple, si on insère la structure NaCI sous forme 
d’une structure cubique primitive avec quatre atomes de Na et quatre atomes de CI, SGROUP 
crée automatiquement une cellule primitive cubique face centrée (FCC) avec un seul atome de 
Na et de CI. 

Symétrie : génère l'expansion LM correspondante pour la densité/potentiel, la symétrie du 
groupe de points de chaque atome et les opérations de symétrie du groupe spatial. 

LStart : produit des densités atomiques en résolvant numériquement l'équation radiale de Dirac 
pour les atomes libres. En utilisant les valeurs propres (ou la localisation dans les sphères 
atomiques) de tous les états atomiques, 1l les regroupe en états de noyau et de valence. Il 
sélectionne automatiquement les LO pour les états semi-cœurs et écrit les paramètres d'énergie 
de démarrage Ex dans case.inl (pendant le cycle SCF, ils sont recherchés et adaptés 
automatiquement pour garantir les meilleurs réglages possibles dans tous les cas). 

KGen : produit, dans la partie irréductible du BZ, un maillage k décalé et équidistant ou non à 
une densité spécifiée par l'utilisateur. 

DStart : produit la densité initiale du cycle SCF en superposant les densités atomiques. Ensuite, 
le cycle SCF implique un ensemble complexe de procédures utilisant plusieurs applications dans 
Wien2k. Les principales étapes sont les suivantes : 

LAPW0 : calcule les potentiels de Coulomb et XC à partir de la densité. 


LAPW1 : calcule les valeurs propres de valence et de semi-noyau ainsi que des vecteurs propres 
à tous les k-points demandés. 


LAPW2 : calcule la densité électronique de valence. 

LCORE : calcule les valeurs propres du noyau et la densité du noyau. 

MIXER : additionne les densités de base et de valence et mélange la densité totale avec les 
densités des itérations précédentes. De plus, 1l peut mettre à jour les positions atomiques en 
fonction des forces calculées, ainsi que les matrices de densité ou les potentiels orbitaux lorsque 


des méthodes DFT+U ou hybrides sur site sont utilisées. 
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CHAPITRE II 


Oxydes magnétiques dilués à large gap et électronique de spin 


L'objectif de ce chapitre est de présenter de manière accessible les phénomènes physiques impliqués dans la spintronique et ses 


applications. Sont également décris l'état de l'art de l'oxyde de magnésium ainsi que les motivations de cette thèse. 


"Tout ce que les êtres vivants font peut-être compris par les tremblements et les mouvements des atomes " 


Richard Feynman 
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Partie A : Spintronic et DMO 


IL.1. Introduction 


La spintronique, un domaine de recherche important, se développe actuellement dans 
diverses nouvelles orientations en physique de la matière condensée. L’objectif du chapitre est de 
présenter de manière accessible les phénomènes physiques impliqués dans la spintronique, en 
évitant les équations complexes au profit de schémas explicatifs. Nous avons également décrit les 
applications de la spintronique, celles d'aujourd'hui et celles qui sont susceptibles d'avoir un impact 
important sur les développements futurs dans les technologies de l'information et de la 


communication. 


La spintronique est généralement définie comme forme innovante de l’électronique qui 
manipule les électrons en agissant non seulement sur leur charge, mais également sur leur spin. 
Son émergence trouve ses origines dans la découverte de la magnétorésistance géante (GMR) dans 
les multicouches magnétiques en 1988 [1]. Durant les premières expériences de la spintronique, les 
courants polarisés en spin sont générés en exploitant l'influence de l'orientation du spin sur les 
propriétés de transport des électrons dans les conducteurs ferromagnétiques. Cette influence, 
suggérée pour la première fois par Mott [2], est démontrée expérimentalement et théoriquement 
décrite dans les premiers travaux [3, 4] bien avant la découverte de la GMR il y a plus de dix ans. 
Ce type de production de courants polarisés en spin par des matériaux magnétiques est utilisé dans 
ce qu'on pourrait appeler la "spintronique classique ” au cours de la première décennie après la 


découverte de la GMR. 


Les événements majeurs de cette époque comprennent les découvertes de la 
magnétorésistance tunnel (TMR) et du couple de transfert de spin. C'est aussi l'époque où des 
concepts importants tels que l'accumulation de spin et le courant de spin pur (courant de spin sans 
courant de charge) sont introduits. Plus récemment, il est devenu possible de générer des courants 
polarisés en spin et des courants de spin purs sans recourir à des matériaux magnétiques en 
exploitant les interactions spin-orbite présentes dans des matériaux non magnétiques. Cela s'inscrit 


dans un domaine émergent appelé spin-orbitronique. 
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IL2. Electronique de spin 
IL.2.1. Pourquoi la spintronique ? 


La miniaturisation des composants microélectroniques a connu une réduction d'environ 
40 fois depuis les débuts des circuits intégrés, vers les années 1970. Au cours de cette période, la 
microélectronique a progressé depuis les premiers circuits intégrés jusqu'aux puces informatiques 
actuelles contenant jusqu'à 100 millions de transistors. Il est désormais reconnu qu’une nouvelle 
réduction de la taille physique des composants électroniques à semi-conducteurs se rapproche 
bientôt d’un obstacle fondamental. Les lois physiques qui régissent le comportement des transistors 
empêcheront leur réduction supplémentaire et leur regroupement en nombre encore plus grand sur 
des puces informatiques. Le rétrécissement continu des transistors entraînera divers problèmes liés 


aux fuites de courant électrique, à la consommation de la chaleur et d'énergie [5]. 


ES 


D'autre part, la miniaturisation des dispositifs électroniques à semi-conducteurs fait 
ressentir aux ingénieurs et aux physiciens la présence imminente de la mécanique quantique, Dans 
ce domaine, des idées contre-intuitives telles que le comportement ondulatoire dominant pour les 
"particules" comme l'électron. Le spin électronique est, après tout, un phénomène quantique. De 
nombreux experts convergent vers l'idée que la spintronique, combinée à la nanotechnologie, 
offrirait la meilleure solution possible aux problèmes liés à la miniaturisation évoqués ci-dessus. 
Les nanosciences et les nanotechnologies englobent l'étude de dispositifs extrêmement petits et des 
phénomènes associés à une échelle spatiale inférieure à un millième du diamètre d'un cheveu 
humain, soit environ la moitié du diamètre d'une molécule d'ADN [6]. Cette échelle nanométrique 
offre de nouvelles perspectives pour la conception de dispositifs électroniques innovants qui 
peuvent surpasser les limites imposées par les contraintes physiques de la miniaturisation classique. 
La combinaison de la spintronique et de la nanotechnologie ouvre ainsi la voie à des avancées 
significatives dans le domaine de l'électronique quantique et de la conception de composants 


électroniques à l'échelle nanométrique. 


IL.2.2. Spintronique des semi-conducteurs 


Malgré les progrès rapides des dispositifs spintronique à base de métaux tels que les 


dispositifs GMR (Giant Magnétorésistance), les chercheurs se sont orientés vers la recherche de 
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nouveaux moyens pour générer et exploiter des courants polarisés en spin dans les semi- 
conducteurs. Ces efforts incluent l'étude du transport de spin dans les semi-conducteurs et 
l'exploration des possibilités d'utiliser les semi-conducteurs en tant que polariseurs de spin et valves 
de spin [7]. Ceci est important car les dispositifs spintronique basés sur les semi-conducteurs 
peuvent être facilement intégrés à la technologie des semi-conducteurs traditionnelle ; ils peuvent 
également servir d'appareils multifonctionnels. De plus, la manipulation et le contrôle des spins 


des semi-conducteurs peuvent être plus facile. 


Les visionnaires envisagent une fusion de l'électronique, de la photonique et du 
magnétique, ouvrant ainsi la voie à de nouveaux dispositifs multifonctionnels basés sur le spin, tels 
que les spin-FET (transistors à effet de champ), les spin-LED (diodes électroluminescentes), les 
spin-RTD (dispositifs tunnel résonants), commutateurs optiques fonctionnant à des fréquences 
Téra Hertz, modulateurs, calcul quantique, etc. [8]. Les progrès réalisés dans ces développements 
dépendent bien entendu de manière cruciale de la compréhension et du contrôle des degrés de 
liberté de spin non seulement dans les semi-conducteurs, mais aussi dans les hétéro structures semi- 


conductrices et les matériaux ferromagnétiques [9]. 


IL3. Magnétorésistance géante (GMR) 


La GMR est le changement brusque de résistance dans les structures magnétiques 
multicouches lorsqu'un champ magnétique externe est appliqué. La découverte fondamentale de 
cet effet spintronique distinctif, est mise au point en 1988 par l'équipe de Peter Grunberg [10] et 
d'Albert Fert à Orsay [11], ce qui a valu à ces chercheurs le prix Nobel de physique en 2007 [12, 
13]. 


Il existe deux manières de mesurer le changement de résistance induit par la GMR. Les 
deux méthodes décrivent ce changement, mais l’une utilise le champ magnétique de saturation 
comme point de référence, tandis que l’autre utilise le point zéro du champ magnétique. La 
première convention s'applique à la configuration antiferromagnétique (antiparallèle) et ne dépasse 
jamais 100 % alors que la seconde convention fait référence à la configuration ferromagnétique 
(parallèle), largement utilisée dans les calculs (figure IL.1). L'équation suivante donne la variation 


de la résistance calculée selon la deuxième convention [14]: 


53 


AR Rar —Rr 


IL 1 
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OÙRF et Ra sont les résistances équivalentes des systèmes dans les configurations 


ferromagnétique et antiferromagnétique, respectivement. 
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Figure IL.1 : Schéma de principe de fonctionnement de la GMR (A) Magnétisation parallèle (B) 


Magnétisation antiparallèle. 


L’important est que pour des changements importants de résistance, les deux méthodes se 
distinguent davantage l’une de l’autre. Les structures elles-mêmes sont constituées de couches 
alternées magnétiques et non magnétiques. Les couches magnétiques sont constituées de matériaux 
ferromagnétiques tels que cobalt, fer, nickel et leurs alliages tandis que les couches non 
magnétiques sont généralement constituées de chrome, de cuivre ou d'argent. Ces couches non 
magnétiques sont connues sous le nom d’espaceurs. L'épaisseur de ces espaceurs est très cruciale 
pour l'effet GMR [15]. L'effet GMR est plus important lorsque l’arrangement antiferromagnétique 
des aimantations dans les multicouches évolue vers un arrangement ferromagnétique. Cette 
transition est significative parce que la GMR dépend de l’orientation relative des magnétisations 
des couches magnétiques successives et donc de la correspondance de la structure électronique 


entre ces couches. 


La GMR à fait l'objet de nombreuses avancées technologiques et trouve une large gamme 
d'applications. Parmi celles-ci, les capteurs GMR (également connus sous le nom de Spin Valves), 


les commutateurs Spin, la détection de mines terrestres, l'électronique non volatile telle que la 
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MRAM, et bien d'autres encore. Notamment, Les têtes GMR utilisées pour la récupération de 
données sont de plus en plus sensibles, ce qui produit des densités de données plus élevées pour les 
disques durs. Le principe sous-jacent à ces applications repose sur le mécanisme de la GMR, qui 


permet une sensibilité plus élevée en raison du changement de résistance [11]. 


IL.4. Magnétorésistance tunnel (MRT) 


Dans une jonction tunnel magnétique (MTJ) où une barrière tunnel ultra-fine sépare deux 
électrodes ferromagnétiques, la MRT fait référence au changement de résistance de la jonction, en 
fonction de la configuration magnétique des électrodes (aimantations parallèles ou antiparallèles). 
Les premières observations de MRT remontent à 1975 avec les travaux de Jullière [16] et en 1982 
avec ceux de Maekawa et Gäfvert [17]. Ces travaux pionniers n’ont pas immédiatement suscité un 
suivi substantiel, probablement en raison de la difficulté à les reproduire. En 1995, des MRT 
importants (environ 20 %) et reproductibles ont été rapportés par les groupes de Moodera [18] et 
Miyasaki [19] sur les MTJ utilisant une barrière tunnel constituée d'une couche d'Al ultra-mince 


oxydée. 


Ces contributions originales ont suscité une activité intense dans le développement des 
MTJ. Une étape importante fut la transition des MTJ utilisant des barrières tunnel amorphes 
(principalement l'aluminium oxydé) vers des MTJ monocristallins, utilisant notamment une 
barrière en M£O. Alors que les premiers résultats obtenus avec la MTJ à base de M£0O n'étaient 
que légèrement supérieurs à ceux des barrières d'alumine [20], une avancée décisive est réalisée en 
2004 avec le premier rapport de taux de MRT élevés (— 200 % à température ambiante) dans les 
MTTJ épitaxiales chez Tsukuba [21] et IBM [22]. Depuis 2004 des progrès continus sont réalisés 
avec les MTJ à base de M£O et les taux de MRT ont augmenté au-dessus de 600 %. 


Le taux de magnétorésistance tunnel est donné par l’expression suivante : 


G,-G Rap — R 2 
MRT= 2 "2 ap Ÿp _ _#PiPz (II. 2) 
Gap Rh 1— P1P2 
Où Gy (Gap) est la conductance de la JMT lorsque les aimantations sont en configuration parallèle 
(antiparallèle), R; (Rap) est la résistance correspondant à un état parallèle (antiparallèle) des 


aimantations des électrodes notées 1 et 2. 
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L'origine des taux élevés de MRT dans MTJ ayant une barrière monocristalline telle que 
M£0 à été attribuée à la sélection de symétrie [23, 24]. Au niveau d'une interface de haute qualité 
entre une électrode magnétique et une barrière tunnel, une fonction d'onde évanescente d'une 
symétrie donnée est connectée aux fonctions de Bloch de même symétrie et de même direction de 
Spin au niveau de Fermi des électrodes. Dans les barrières M£O orientées selon (001), l'une des 
symétries correspondant aux bandes fortement polarisées en spin présente une désintégration 


beaucoup plus lente, ce qui se traduit par une MRT plus élevé. 


D'autres règles de symétrie sont observées dans des barrières constituées de différents 
matériaux, conduisant à des propriétés différentes telles que la polarisation de spin effective 
négative dans les MTJ à base de SrTiO3 [25]. D'autres axes de recherche sont explorés afin 
d'augmenter davantage les taux des MRT, comme le remplacement des électrodes 
ferromagnétiques conventionnelles par des demi-métaux 100% polarisés en spin constitués de 
matériaux oxydes [26] ou d'alliages Heusler [27], ou encore l'utilisation de barrières 


ferromagnétiques en tant que filtre rotatif (figure IL. 2). 


Ferromagnétique 


Isolant 
(a) (b) 


Faible résistance Forte résistance 
électrique électrique 


Figure IL 2 : Etats parallèles et antiparallèles dans les MRT. 


D'un point de vue applicatif, les MRT sont introduits par l'industrie dans les têtes de lecture 
des disques durs (HDD) en 2006, principalement parce qu'ils apportaient des améliorations en 
termes de rapport signal/bruit. Ceci a permis de soutenir l'augmentation continue des densités 
surfaciques des informations stockées et en même temps, la réduction des coûts de ce stockage. Un 
autre domaine d’applications impacté par la MRT est l’émergence des mémoires magnétiques à 
accès aléatoire (MRAM). Alors que le concept de MRAM avait déjà été introduit par Honeywell 
avant la découverte de la GMR [28, 29], son développement n'a été déclenché qu'à la fin des années 


1990 grâce au MRT. 
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En effet, la configuration naturelle des blocs de construction MTJ, avec un courant 
circulant perpendiculairement aux couches, s'est révélée particulièrement bien adaptée au 
développement d'architectures compactes telles que les barres transversales. Cependant, dans la 
première génération de MRAM basée sur les MRT développée au milieu des années 2000, l'écriture 
des informations est réalisée grâce à des champs magnétiques générés par des lignes de "mot" et 
de "bit", ce qui s'est avéré coûteux en termes de consommation d'énergie. Ce problème est résolu 
par l’introduction de l’écriture du couple de transfert de spin [30]. Cela a mené les recherches 
ultérieures vers un contrôle électrique de l’aimantation, en utilisant le couple de transfert de spin 


et les matériaux multiferroïques qui combinant les propriétés ferroélectriques et ferromagnétiques. 
ILS. Applications de la spintronique 


IL.5.1. Mémoires magnétiques à accès aléatoire MRAM 


La mémoire magnétique à accès aléatoire (MRAM) est une classe de circuits de stockage 
à semi-conducteurs qui enregistrent les données sous forme d'états magnétiques stables de 
dispositifs magnétorésistance. Pour lire les données, la mesure de la résistance des dispositifs est 
nécessaire pour déterminer leurs états magnétiques (figure IL. 3 (A)). En pratique, les dispositifs 
magnétorésistance sont intégrés à des circuits métal, oxyde, semi-conducteur complémentaire 
(CMOS) pour créer des puces compatibles avec l'électronique à semi-conducteurs produite en série 
[31]. Bien que divers dispositifs magnétorésistifs sont utilisés dans la conception de ces circuits, 
les produits MRAM disponibles commercialement, ainsi que la grande majorité des technologies 
MRAM en cours de développement pour les futures générations, sont basés sur des dispositifs à 
MTJ. Tous ces circuits fonctionnent comme des mémoires résistives au niveau de l'opération de 
lecture ; c'est la méthode d'écriture de l'état magnétique qui distingue les différents types de 
technologie MRAM. Certaines des méthodes d'écriture les plus étudiées incluent la commutation 
de champ de type Stoner-Wolfarth, la commutation de Savtchenko (également une méthode de 
commutation de champ), la commutation de couple de rotation et la commutation à assistance 


thermique (chaleur avec champ ou couple de rotation) [32]. 


Jusqu'à présent, deux méthodes d'écriture sont commercialisées : la MRAM à bascule, qui 
utilise la commutation de Savtchenko, est en production de masse depuis 2006 et la MRAM à 


commutation de couple de rotation est aux premiers stades de la production commerciale [21]. 


57 


IL.5.2. Diodes électroluminescentes (Spin LED) 


La spin LED est un dispositif optoélectronique de spin qui convertit les informations de 
spin contenues dans une population de porteurs polarisés en spin en lumière polarisée 
circulairement. Les LED sont fortement demandées dans les domaines tels que la signalisation, 
l'affichage avec des écrans couleurs, les sources de lumière blanche et en particulier dans l'éclairage 
public. Dans les sources de lumière polarisées en spin, les électrons polarisés sont injectés depuis 
une couche non magnétique dans une structure semi-conductrice non magnétique par dérive et les 
mécanismes de diffusion sont ajustés [33]. Ils se recombinent radiativement avec des trous non 
polarisés injectés depuis un contact non magnétique dans le milieu actif de la structure. S1 la durée 
de vie du spin des porteurs est supérieure au temps de recombinaison dans le milieu actif, alors 
l'orientation du spin ne se détend pas complètement au moment de la recombinaison [34]. Ainsi, le 
rayonnement issu de la recombinaison des porteurs polarisés en spin sera partiellement polarisé 
circulairement selon les règles de sélection quantique optique. A la différence d’une LED 
conventionnelle, en revanche, les électrons et les trous non polarisés se combinent pour produire 
une lumière non polarisée, la Spin LED exploite la polarisation du spin pour générer une lumière 


partiellement polarisée circulairement (figure IL 3 (B)). 


IL.S.3. Transistors à effet de champ à polarisation de spin (Spin FET) 


Les applications potentielles des semi-conducteurs ferromagnétiques incluent la 
conception de dispositifs nouveaux et améliorés, tels que des transistors de spin. Ceux-ci devraient 
être utilisés comme élément de base de circuits logiques à faible consommation d’énergie, non 
volatiles et reconfigurables. En 1990, Datta et Das ont réintroduit un FET à spin, qui nécessite une 
injection de spin bien organisée dans un semi-conducteur. Premièrement, les spins électroniques 
sont injectés à partir d’une source. Ensuite, les spins sont modulés par la polarisation de grille et 


sont finalement détectés au niveau d'un drain [35]. 


La logique basée sur le spin est actuellement étudiée comme l'une des technologies 
informatiques au-delà du CMOS pour compléter les transistors à effet de champ (FET) à semi- 
conducteurs à oxyde métallique complémentaires dans les microprocesseurs. Des efforts de 
recherches déterminés dans le domaine de la spintronique ont abouti à la commercialisation de la 
MRAM et de la STT-MRAM (mémoires magnétique a accès direct par couples de transfert de spin) 


ayant ainsi une immense influence sur la vie quotidienne [21]. La question qui se pose maintenant 
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est de savoir si les dispositifs spintronique sont suffisamment aptes pour obtenir un succès similaire 


dans le domaine de la logique. 


L'un des principaux concurrents dans ce domaine pourrait être le transistor à effet de 
champ de spin (spin FET) (figure IL. 3. (C)), qui est essentiellement un dispositif à trois bornes 
utilisant des contacts ferromagnétiques dans la source et le drain comme injecteur de spin et 
détecteur. Le concept du Spin FET, visualisé par Datta et Das en 1990, ouvre une passerelle vers 
le traitement de l'information basé sur le spin (Datta 2018 ; Datta et Das 1990) [36]. Bien que les 
recherches sur la spintronique des semi-conducteurs ait rendu possible la manipulation cohérente 
des spins électroniques dans les semi-conducteurs, la réalisation pratique d'un transistor à effet de 
champ de spin fonctionnel pour le traitement de l'information n'a pas encore été réalisée. Les 
principaux défis sont la faible efficacité de l'injection de spin en raison de l'inadéquation des 
résistances, de la relaxation du spin et de la propagation des angles de précession de spin. Supriyo 
Datta a déclaré s'être inspiré des phénomènes optiques pour proposer le concept de transistor de 
spin. Un dispositif électronique similaire aux dispositifs optiques a été proposé en exploitant les 


similitudes entre le spin d'un électron et la polarisation d'un photon. 


IL.5.4 Diodes laser (LD) 


Le laser est un rayonnement électromagnétique présentant des caractéristiques distinctives 
telles qu'une cohérence, une mono chromaticité, une directivité et une intensité toutes élevées. C’est 
l’une des inventions exceptionnelles du XXe siècle et elle a révolutionné des domaines diversifiés 
de la science et de la technologie. On dit souvent que "où que vous alliez et quoi que vous fassiez 
: le Laser vous suivra”. Le laser a également donné naissance à un domaine spécialisé en 
spintronique l'opto-spintronique. Il est important de noter que cela concerne la lumière laser 
conventionnelle qui utilise uniquement les charges des porteurs actuels Il est intéressant de noter 
que les lasers à spin pourraient ouvrir la voie à plusieurs dispositifs spintronique réalistes à 
température ambiante, qui ne se limiteront pas à la seule magnétorésistance. Les diodes laser 
partagent les avantages des LED, mais émettent une lumière laser cohérente et unidirectionnelle. 
Elles sont utilisés dans les pointeurs laser et dans des applications scientifiques et industrielles 
spécialisées (pompage optique d'autres lasers, spectroscopie, trempe superficielle, soudage). Les 


industries des télécommunications et du stockage de données optiques sont les plus gros 
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consommateurs de diodes laser, en partie parce que les diodes laser ont un temps de réponse 


extrêmement rapide, elles peuvent être allumées et éteintes très rapidement [37]. 


IL.5.5. Détecteurs UV 


La couche d'ozone et les gaz de l'atmosphère absorbent une grande partie de l'énergie 
ultraviolette émise par le soleil. Seul un rayonnement de plus de 280 nm atteint la surface de la 
terre. Les détecteurs de rayonnement UV, également appelés "store solaire" produits sur terre, 
doivent détecter les rayonnements dans la bande de 265 nm à 280 nm, caractérisée par un minimum 
de rayonnement parasite. La détection UV a des applications militaires ou civiles telles que les 
dosimètres personnels pour les environnements riches en UV, la détection des incendies et 
l'identification des missiles par leur traînée, ou le guidage des missiles. Les photodétecteurs 
actuellement disponibles sont à base de SiC (Eg=2.9 eV - 430nm) ou de couches minces de diamant 
(Eg=5.4 eV — 230 nm). Cependant, ils ont une mauvaise sélectivité UV/visible et nécessitent 
l'utilisation de filtres. L'alliage AllnGaN dont la bande interdite peut varier de 1.9 eV (650 nm) à 
6.2 eV (200 nm), offre une solution prometteuse pour la détection optique. L'alliage AIGaN peut 
être utilisé pour la détection UV, avec une absorption de 365 nm pour GaN et à 200 nm pour AIN. 
En outre, les nitrures III-V (AlnGaN ) sont stables dans des conditions physiques et chimiques 


extrêmes [38]. 
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Figure IL. 3. (A) Mémoire MRAM (Magnetic Random Access Memory), (B) diode 
électroluminescente (LED) à spin [39], (C) Transistor à effet de champ de spin, (D) Détecteur 
UV. 


IL. 6. Demi-métaux 


Les ferromagnétiques demi-métalliques (HMEF) constituent une voie prometteuse pour le 
développement des matériaux hautement polarisés en spin. Cette classe de matériaux est 
initialement proposée par de Groot et coll. [40] au début des années 1980. Contrairement aux 
ferromagnétiques conventionnels, la polarisation de spin résulte d’un déséquilibre dans la densité 
d’états entre les électrons de spin ascendant (majoritaire) et de spin descendant (minoritaire) [41]. 
Dans les HMF, les propriétés de conduction du canal minoritaire sont complètement différentes de 
celles du canal de spin majoritaire. La bande majoritaire présente une occupation complète des 
états électroniques jusqu'à l'énergie de Fermi, donnant une conduction métallique tandis que les 
états minoritaires ont une bande interdite similaire à celle d'un semi-conducteur. Ceci est représenté 


schématiquement sur la figure IL. 4. 


Depuis les premières recherches sur le ferromagnétisme demi-métallique, notamment sur 
NiMnSb, de nombreux composés sont identifiés comme étant demi-métalliques. Parmi lesquels, 
de nombreux sont des alliages complexes, principalement des alliages Heusler. Bien que le 
ferromagnétisme demi-métallique est observé dans le CrO2 et même prédit dans le graphène, ces 
matériaux présentent des propriétés de conduction inhabituelles. En raison de ces propriétés de 
conduction inhabituelles, les HMEF pourraient potentiellement fournir un matériau avec une 
polarisation de spin de 100 %. Cela fournirait des valeurs extrêmement élevées de TMR et de GMR 
(80 %) [42, 43]) ainsi que des couples de transfert de spin importants pour les futurs dispositifs 
spintronique. 
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Figure IL. 4. Diagrammes schématiques des structures de bandes pour un paramagnétique, 


ferromagnétique et un semi-ferromagnétique métallique [441]. 
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IL. 7. Oxyde magnétiques dilués (DMO) 


Les oxydes magnétiques dilués (DMO) sont des semi-conducteurs non magnétiques dopés 
avec une faible concentration d'éléments magnétiques, généralement des métaux de transition. 
L'ajout intentionnel de ces éléments magnétiques dans la matrice semi-conductrice crée des centres 
magnétiques localisés qui peuvent induire le ferromagnétisme dans le matériau, tout en conservant 
les propriétés semi-conductrices de base. Ces matériaux sont conçus pour être facilement 
intégrables aux semi-conducteurs existants et sont également connus sous le nom de semi- 
conducteurs magnétiques dilués (DMS) [45], en raison de leur forte polarisation en spin. 
Cependant, la découverte et la compréhension de tels matériaux constituent un défi de taille pour 


la science du solide. 


Les premières DMO identifiés qui présentent des propriétés ferromagnétiques à des 
températures supérieures à l’ambiante sont découverts par Dielt et al. dans les années 2000 [46]. 
Cette recherche théorique a impliqué le dopage de l’oxyde ZnO avec du manganèse pour obtenir 
des semi-conducteurs de type p. Suite à cette étude, d’autres chercheurs ont étudié 
expérimentalement les DMO en 2001 en utilisant des couches minces de ZnO et T102 dopé au Co 


[47, 48]. 


Dans l'ensemble de ces matériaux, 1l a été prouvé que le ferromagnétisme est médié par 
les porteurs, offrant ainsi la possibilité de modifier le comportement magnétique par la 
manipulation des charges. Cela a motivé une recherche continue de matériaux présentant une 
température de Curie (Tc) encore plus élevée et un ferromagnétisme médié par les porteurs et a 
conduit à la conjecture selon laquelle les DMS à base d'oxydes seraient des matériaux clés dans le 
développement de dispositifs spintronique [49]. En effet, 1l est important de souligner que la 
capacité de dopage électronique élevé et la masse électronique effective plutôt élevée des semi- 
conducteurs à oxyde pourraient être très efficaces pour atteindre des températures de Curie élevées. 
De plus, la plupart des DMS prévus à base d'oxydes sont des semi-conducteurs à large bande 
interdite (> 3eV) qui peuvent ajouter une dimension optoélectronique à la nouvelle génération de 


dispositifs spintronique. 


Les oxydes magnétiques dilués offrent une combinaison unique de propriétés magnétiques 
et semi-conductrices, ouvrant la voie à de nouvelles possibilités d'applications dans divers 


domaines technologiques. 


62 


IL. 8. Mécanismes sous-jacents aux propriétés magnétiques des DMO 


Dans le cas des oxydes magnétiques, la plupart des interactions magnétiques sont souvent 
attribuées à des mécanismes de double échange ou super échange. Cette caractéristique est liée à 


la nature particulière des matériaux et aux interactions entre les spins des ions magnétiques. 


IL. 8. 1. Mécanisme de double échange 


Le mécanisme de double échange, initialement proposé par Zener constitue une 
interaction magnétique entre deux ions métalliques à travers un ion non magnétique (anion) 
intermédiaire, jouant un rôle déterminant dans les propriétés magnétiques de certains matériaux, 
notamment les oxydes [50]. Ce mécanisme implique le transfert d'électrons d'un ion métallique à 
un autre via l'anion intermédiaire, tout en conservant le spin de l'électron pendant le processus 
figure IL 5. (A) l'interaction d'échange entre les spins des deux cations est dictée par l'orientation 
relative de ces spins, conduisant à deux configurations possibles : ferromagnétique, où les spins 
des cations sont parallèles, et antiferromagnétique, où les spins sont antiparallèles. L'énergie 
d'échange, déterminée par cette orientation, gouverne l'état magnétique du matériau : une énergie 
d'échange positive entraîne un comportement ferromagnétique, tandis qu'une énergie d'échange 


négative induit un comportement antiferromagnétique [51]. 


IL. 8. 2. Super échange 


Le super échange, un type d'interaction magnétique indirecte, intervient entre deux ions 
métalliques (cations) à travers un ion non magnétique (anion) intermédiaire, a été développé par 
John H. Van Vleck dans les années 1930 pour expliquer les propriétés magnétiques de certains 
cristaux se distinguant de l'échange direct par l'implication d'un anion. Ce mécanisme implique le 
transfert d'électrons, où un électron d'un cation peut se déplacer vers l'anion intermédiaire. Le 
couplage de spin survient lorsque le spin de l'électron délocalisé s'aligne avec le spin du cation 
d'origine, influençant ainsi l'alignement des spins des deux cations voisins. Les types de super 
échange incluent le super échange ferromagnétique, où les spins des cations s'alignent 
parallèlement, favorisant une interaction ferromagnétique, et le super échange antiferromagnétique 
[52, 53], où les spins s'alignent antiparallèlement, favorisant une interaction antiferromagnétique 


(figure IL 5. (B)). 


63 


Plusieurs facteurs influencent le super échange, notamment la nature des cations et de 
l'anion, la distance entre les cations, et la structure cristallographique du matériau. Ce type 
d'interaction s'oppose généralement aux propriétés ferromagnétiques des matériaux DMS et 


devient plus prononcé à mesure que la concentration des ions augmente. 


À al L. : V4 d N B d-orbital d-orbital 
_ | \ | \_ p-orbital 
1 eg : \\ : _—_—_—_——— eg | t | 
———|— CD l | 
T0 0 y 0 

Mn?- (d#) Mné- (d°) Ni Ni 


Figure IL. 5. Représentation de mécanisme : (A) du double échange, (B) Super échange direct. 
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Partie B : Etat de l’art et motivation de la thèse 


IL.1 Structures cristallines 


Les réseaux sont des arrangements réguliers de points imaginaires dans l’espace. Un 
cristal réel est constitué d’atomes associés à ces points. Chaque position dans le réseau cristallin 
joue le rôle de base pour le cristal, où chaque point du réseau est lié à des atomes spécifiques de 
l'élément considéré. Cela définit la structure fondamentale du cristal, où chaque point du réseau est 


associé à des atomes particuliers, formant ainsi la base structurale du matériau cristallin. 
Réseau + Base = Structure Cristalline 


Les structures cristallines désignent l'organisation ordonnée et répétitive des atomes, ions 
ou molécules dans un matériau cristallin. Parmi les structures cristallines cubiques il existe trois 
types telle que : la structure cubique simple, la structure cubique centrée et la structure cubique a 
face centrée. Chaque type présente des arrangements distincts des constituants dans l'espace 


tridimensionnel. 


II.1.1 Structure cubique simple (CS) 


La structure cubique simple est l'une des structures cristallines fondamentales, étant la 
plus simple et la moins dense parmi les structures cristallines cubiques (figure IL. 6 (A)). Dans 
cette configuration, chaque atomes est situé à un coin du cube et est en contact avec six atomes 
voisins par unité de cellule. Il est important de noter que dans la structure CS, les atomes n'occupent 


que 52.4 % de l’espace disponible dans la cellule unitaire [54]. 


IL.1.2 Structure cubique centré (CC) 


La structure cubique centrée (CC) est une autre configuration cristalline où, similairement 
à la structure cubique simple (CS) chaque atome est situé à un cois du cube, Cependant, à la 
différence de la structure SC, la structure CC possède un atome supplémentaire placé au centre du 
cube, comme illustré dans la figure IL. 6 (B) Dans cette structure chaque atome est en contact avec 
huit atomes voisins : six à travers les faces du cube et deux à travers les diagonales du cube. Les 
atomes remplissent environs 68 % de l’espace disponible dans la cellule unitaire de la structure 


CC, ce qui en fait une structure plus dense que la structure CS. La présence courante de cette 
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structure dans de nombreux matériaux métalliques, s'explique par sa capacité à conférer une bonne 


résistance et une bonne ductilité des métaux. 


IL.1.3 structure Cubique à faces centrées (CFC) 


Structure cubique a faces centrées (CFC), est une autre structure cristalline dans laquelle 
chaque atome est situé non seulement à chaque coin du cube, mais également au centre de chaque 
face du cube (Figure IL. 6 (C)). Dans cette disposition, chaque atome est en contact avec douze 
atomes voisins : quatre à travers les arêtes du cube et huit à travers les faces du cube et remplissent 
environ 74 % de l’espace disponible dans la cellule unitaire, ce qui en fait une structure très dense. 
Cette configuration crée une alternance de couches atomiques. La structure cubique à face centrées 
est couramment observée dans de nombreux métaux, tels que l’argent, l’or, le cuivre et 


l'aluminium. 


Figure IL 6. (A) La structure cubique simple (CS), (B) La structure cubique centrée (CC), (C) la 


structure cubique à faces centrées (CFC). 


IL 2. Matériaux isolants conducteurs et semi-conducteurs 


Dans le domaine des matériaux, une classification fondamentale repose sur la distinction 
entre conducteurs, semi-conducteurs et isolants, un système qui dépend étroitement de la taille de 
la bande interdite (Figure IT. 7). Cette bande interdite, également connue sous le nom d'espace 
énergétique, joue un rôle déterminant dans le comportement électronique des substances. Les 
conducteurs se caractérisent par une bande de valence partiellement remplie, favorisant une 
mobilité électronique élevée et une circulation fluide du courant électrique. En revanche, les semi- 
conducteurs présentent un faible écart énergétique entre leurs bandes de conduction et de valence, 
généralement inférieur à 4 eV, ce qui leur confère une capacité de régulation du courant électrique. 
Les isolants, dotés d'une bande interdite considérablement large, exigent une quantité significative 
d'énergie pour permettre aux électrons de franchir la barrière entre la bande de valence et la bande 


de conduction. 
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Lorsque des électrons quittent la bande de valence, ils créent des "trous", des 
emplacements vacants pouvant également contribuer à la conduction électrique. La compréhension 
approfondie de ces propriétés électroniques, liée aux principales bandes d'énergie telles que la 
bande de conduction, la bande de valence et l'espace énergétique, est essentielle pour concevoir 


des matériaux adaptés à des applications spécifiques. 


ntercit 
erdite 


Isolant Semi-conducteur Conducteur 


Figure IL. 7. Diagrammes de bandes d’énergie de différents types de matériaux. 
IL. 3. Généralités sur l'oxyde de magnésium (MgO) dans la littérature 


IL. 3. 1. Pourquoi l’oxyde de MgO de la famille II-VI 


Au cours des dernières années, les oxydes alcalino-terreux y compris les semi-conducteurs 
et les oxydes magnétiques dilués (DMS et DMO), ont attiré l'attention dans le domaine industriel 
et sont considérés comme des matériaux prometteurs pour diverses applications, notamment la 
spintronique, la catalyse [55], la microélectronique et l’optoélectronique (laser et diodes 
électroluminescentes (LED) [56-58]. MgO a un large éventail d’applications pratiques. 
Photocatalyseur, réfractaire [59], la catalyse [55], les applications de stockage d'énergie [60, 61], 
les produits supraconducteurs, l'industrie des matériaux et la médecine [62] etc. Ces applications 


diverses soulignent la polyvalence et l'usage intensif de MgO dans divers secteurs. 


Cette thèse, se concentre spécifiquement sur l'oxyde de magnésium (M£gO) en raison de 
ses propriétés physicochimiques critiques : transmission optique élevée, faible densité (3.58 
g/cm°), faible coût, résistance aux températures élevées, nature respectueuse de l'environnement et 


énergies de liaison des excitons élevées (80- 85 meV) [63] (voir Tableau IL. 1.). 
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L'oxyde de magnésium, ou magnésie, est l'un des éléments répandus dans le manteau 
inférieur de la Terre. Il possède une large bande interdite expérimentale égale à 7.8 eV [64]. Dans 
la nature, il se trouve sous différentes structures telles que les structures NaCL, et Zinc Blende. Dans 
des conditions normales, l'oxyde M£O est un semi-conducteur stable dans la structure cubique de 


type NaCI (B1). 


L'oxyde de magnésium est utilisé comme adsorbant pour agents chimiques toxiques pour 
l'élimination des ions de métaux lourds tels que Pb (I), As (D), Ni (I) et Cd (I), car il s'agit d'un 
matériau abondant, non toxique et respectueux de l'environnement [65], Sa capacité à adsorber ces 
substances toxiques est attribuée à sa nature abondante, non toxique et respectueuse de 
l'environnement. De plus, en modifiant sa morphologie pour augmenter sa surface, sa capacité 
d'élimination des métaux lourds peut être améliorée. Il est intéressant de noter que M£O présente 
une capacité d’adsorption exceptionnellement élevée pour l’arsenic, par rapport à celle des autres 


oxydes métalliques, 


Les nanoparticules d'oxyde de magnésium possèdent également des propriétés 
antibactériennes [66] avec une efficacité notable contre les bactéries Gram positives, en particulier 
S. aureus, et une activité comparable contre les bactéries Gram négatives telles qu'E. Coli. Ces 
propriétés antibactériennes ouvrent la voie à diverses applications potentielles dans le domaine 
médical, pharmaceutique ou même dans des applications liées à la prévention de la croissance 


bactérienne dans divers environnements [59, 67]. 


Des études sur les propriétés de l'oxyde de magnésium pure et dopé ont fait l'objet d'une 
attention significative dans la littérature scientifique (expérimentales et théoriques) [68,71]. Il est 
rapporté par Balakrishnan et al. et Taiïb et al. [72, 73] que la bande interdite de M£gO pur, dans les 
phases B1, est une bande interdite directe et est égale à 5.05 eV et 4.756 eV, respectivement, Des 
études ont signalé l'existence du ferromagnétisme dans l'oxyde de magnésium. Plusieurs travaux 
ont confirmé ce phénomène, mettant en évidence la présence de propriétés magnétiques dans M£0O. 
D'autres études ont signalé l'existence du ferromagnétisme dans l'oxyde de magnésium [74-77] et 


ont confirmé ce phénomène, mettant en évidence la présence de propriétés magnétiques. 


Il a été rapporté que le dopage de M£0 avec des éléments du groupe (I, IV, V) du tableau 
périodique et des métaux de transition présente un comportement de ferromagnétisme semi- 


métallique [78]. Expérimentalement, Bilal et al. ont réussi à évaluer les effets du dopage à l'argent 
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sur les propriétés structurales et optiques de MgO en utilisant des techniques de caractérisations 
telles que SEM, EDX, spectroscopie, XRD, FTIR et UV-Vis. L'analyse des données indique que 
le photo catalyseur à l'argent supporté par la magnésie est efficace pour l'élimination du bleu de 
méthylène de l'eau contaminée [79]. Une étude menée par Obeid et al. [80] a porté sur la synthèse 
d'oxyde de magnésium pur et dopé au cobalt en utilisant la méthode de co-précipitation. Les 
résultats obtenus révèlent que les systèmes MgO dopés au Co sont mécaniquement et 
thermodynamiquement stables. De plus, la formation d'un MgO monophasé de taille nanométrique 
pour M£0 pur et dopé a été associée à un redshift lorsque le dopage au cobalt augmente suggérant 


des applications potentielles en photocatalyse. 


Actuellement, une grande partie des travaux expérimentaux se concentre sur MgO dopé 
aux métaux de transition. Par exemple, Narayan et al. ont élaboré des cristaux de M£O dopés au 
Ni en utilisant la technique de fusion à l'arc. Ils ont mis en évidence une propriété paramagnétique 
à des concentrations de dopant de 0.5 % et 4 %, en l'absence de défauts ou porteurs de charge. 
Dutta et al. ont synthétisé du M£O dopé à 10 % en Cu en utilisant la technique d'imprégnation 
primaire a humidité [81]. Ils ont constaté que M£O dopé au cuivre est également paramagnétique. 
Azzaza et al.[82] ont étudié le comportement magnétique et les propriétés optiques de MgO 
substitué par Mn, L'étude a révélé un petit ordre ferromagnétique pour des concentrations étroites 


de Mn de 10 à 20 %. 


D’autres études théoriques ont été menées par Seike et al. [83] sur les propriétés 
magnétiques de composés tels que MgO, Ca0O et SrO dopés à l'azote. Les résultats prédictifs ont 
montré qu’en l’absence d’éléments magnétiques, il induisait le ferromagnétisme appelé 
ferromagnétisme d°. Tandis que Nazir et al. ont prédit le ferromagnétisme semi-métallique en co- 
dopant MgO hôte avec des atomes d'halogène en utilisant des calculs de premiers principes [841]. 
Alors que Klaa et al. ont étudié Mgi:;M;O, M (Fe, Ni) dans la phase de sodium en utilisant 
l’approche GGA- mBJ. Les alliages présentent des caractéristiques de semi-conducteur magnétique 
avec une bande interdite directe pour x = 0.25 et 0.5 et des propriétés semi-métallique pour d'autres 


concentrations (x = 0.75) [85]. Le métal de transition modifie les propriétés de la magnésie. 


Pardo et al. [86]ont approfondi l’étude de N en position substitutionnel dans la série des 
monoxydes alcalino-terreux MgO, CaO, SrO et BaO en utilisant la méthode d'onde plane 


augmentée plus orbitales locales (APW + lo). Leurs observations ont révélé une délocalisation 
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accrue et une diminution de la bande interdite à mesure que le rayon cationique augmentait. Le 
couplage magnétique diminue rapidement à mesure que la valeur de U augmente. En fait, l’effet 
de ce dopage non métallique est similaire à celui des lacunes cationiques dans MgO ou Ca. Pour 
de tels systèmes, le mécanisme diffère du dopage conventionnel des métaux de transition dans les 
oxydes semi-conducteurs, où les états d jouent un rôle important dans la formation du 
ferromagnétisme. Le mécanisme exact d'un tel ferromagnétisme semi-métallique d° continue d'être 


fortement débattu dans la littérature. 


Dans une autre étude, les structures électroniques et le magnétisme, du M£O dopé avec 
des éléments des groupes principaux II, IV et V (éléments métalliques et non métalliques) à une 
concentration de 25 % sont systématiquement étudiés. le premier cas, dans lequel des atomes 
dopants remplacent l'atome O ; ensuite la situation de substitution d'un atome de métal à l'atome 


de Me ; enfin, un système M£gO dopé au F [74]. 


Pour le premier cas, dans lequel les atomes dopants remplacent l'atome O, tous les 
systèmes présentent des propriétés ferromagnétiques semi-métalliques où les électrons d n'ont 
aucun effet. Un dopage de type p, qu’il soit métallique ou non métallique dans M£O, présente des 
propriétés ferromagnétiques semi-métalliques tant que l'énergie orbitale np de l'atome dopant est 
inférieure à celle de l'atome dopé np-O et formant des états localisés après le dopage. Un numéro 
atomique différent des atomes dopants a peu d'effet sur le moment magnétique local pour le même 
groupe principal, alors que différents dopants du groupe principal affectent fortement le moment 
magnétique ; Pour les systèmes dopés, la tendance à la division de spin est la plus simple pour 
M£0 dopé au B, alors que M£O dopé au N est le plus stable; MgO dopé au F ne présente pas de 
propriétés magnétiques semi-métalliques, mais le système est un conducteur électrique. Les 
différents types de dopage entraînent des changements spectaculaires dans les propriétés du MgO 


dopé. 


Dans le second cas, la situation de substitution d'un atome de métal à l'atome de Meg, tous 
les systèmes se trouvent dans l’état fondamental paramagnétique et ont une faible énergie de 
formation. Contrairement au premier cas, la conclusion dérivée est que le ferromagnétisme dans 


M£0 dopé est provoqué par un couplage magnétique médié par des trous dans la bande np. 


Sur la base des détails précédents, il est clair qu’il existe une quantité considérable de 


travaux sur les composés susmentionnés impliquant à la fois des méthodes expérimentales et 
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théoriques. Cependant, aucune étude n’a été publiée sur le calcul des paramètres de l’état 
fondamental, la structure électronique et les propriétés optiques de ces composés en utilisant la 


correction mBJ. Ces lacunes dans la littérature nous ont incités à réaliser ce travail. 


IL. 3. 2. Propriétés et applications de l’oxyde MgO pur 
IL. 3. 2. 1. Propriétés physiques et chimiques 


Avant d'approfondir l’étude sur M£0O, il est indispensable d’aborder ses aspects physique 
et chimique généraux. L'oxyde de magnésium également appelé magnésies présente sous la forme 
d'un solide blanc, souvent sous forme de poudre inodore. À l'échelle microscopique, il se manifeste 
sous la forme de cristaux cubiques translucides [87]. C’est un matériau ionique, il doit ses 
propriétés aux interactions entre les atomes d’oxygène et de magnesium, caractérisées par une 


faible électronégativité pour le Mg et une électronégativité plus élevée pour l’oxygène. 


L'oxyde de magnésium peut être formé en faisant réagir le magnésium métallique avec de 
l'oxygène gazeux. Cette réaction chimique est facilement réalisée en brûlant des rubans métalliques 
de magnésium, où le magnésium s'oxyde et produit de l'oxyde de magnésium. Des précautions 
extrêmes doivent être prises si MgO doit être préparé par cette méthode car la réaction est 
exothermique et lorsque le magnésium s'oxyde, il émet une lumière blanche brillante, riche en 
ultraviolets et difficile à éteindre. Les propriétés physiques de base du M£0O sont détaillées dans le 
tableau IT. 1. Son nom minéral est "périclase" et utilisé comme source de magnésium métallique 
[88]. L'oxyde de magnésium est de nature hygroscopique et il faut veiller à le protéger de 


l'humidité. Il réagit avec l'eau et forme l'hydroxyde : 


Figure IL 8. Structures de l'oxyde de magnésium [89]. 
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Tableau IL. 1. Quelques propriétés chimiques et physiques de MgO [901]. 


Paramètre du réseau 4.231(A) 
Masse molaire 40.3044 g/mol 
Masse moléculaire 40.304 g/mol 
Densité 3.58 g/cm° 
Point de fusion 2825 K 


Indice de réfractions 


1.7255 près de 589 nm 
1.7283 près de 750 nm 


Constante diélectrique 


9:65 


Dureté 5.5 à 6 (Mohs) 
Temperature d’ébullition 3600 °C 
Température spécifique a 27 °C 0.9288 KJ/K£g.K 


Module d’élasticité 


Cu = 294, Ci = 93, Ca= 155 


Capacité thermique 


8.88 cal /Kmol 


Coefficient d'absorption 


0.05 cm'! 


Module d’Young 


249 (GPa) 


IL. 3. 2. 2. Structures cristallographiques de l'oxyde de magnésium 


L'oxyde de magnésium adopte une structure cristalline cubique face centrée, similaire à 
celle du sel de chlorure de sodium (NaCT), dans laquelle chaque atome de magnésium est entouré 
de six atomes d'oxygène il appartient au groupe (II-VD) [91]. Les ions magnésium Mg ** et les ions 
oxygène O7 forment une liaison ionique au sein de la structure cristalline, avec des paramètres 


cristallins égaux (a =b =c=4.2130 À) et des angles droits. (Figure IL. 9 (A)) 


L’oxyde de magnésium M£gO est également présent dans d’autres structures cristallines 


différentes, telles que Zinc Blende et la structure CsCI. 


La structure cristalline de la Zinc Blende (ZB) est également de type cubique et est 
représentée par deux cubes à face centrées décalés l’un par rapport à l’autre, avec un groupe 


d’espace F4 3m (Figure IL. 9 (B)). Dans cette structures, les atomes occupent des positions 
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tétraédriques alternées qui signifie qu’un atome est entouré par quatre (4) autres atomes et formant 


un tétraèdre. 


Pour la structure de chlorure de césium (CsCI), c’est un autre type dans les systèmes 
cubiques centré avec un groupe d’espace de Pm3m, avec 8 atomes de césium (Cs) et 8 atomes de 
chlorure (C1) comme illustré dans la Figure IL. 9 (C). Dans cette structure chaque ion Cs* est situé 


au centre d’un cube et entouré par un arrangement cubique d’ions CI qui placé aux huit coins du 


cube. 
(A) (©) 
19 —@ 
@— e” 
8 
@ Cs: 
ea &-* —". 


Figure IL. 9. Les structure cristalline de l'oxyde de magnésium : (A) (NaCT), (B) ZB, et (C) CsCI. 
IL. 3. 3. 3. Propriétés électroniques et optiques 


Selon les structures de bande électronique du MgO, il se classe en tant qu’un isolant 
électrique avec une largeur de bande interdite élevée d’environs 7.8 eV. Ce qui signifie qu’il 
nécessite une quantité significative d'énergie pour permettre le passage des électrons à travers le 
matériau. Cela lui confère une faible conductivité. D’autre, la capacité d'isolation électrique élevée, 


M£0 est largement utilisé dans la fabrication de tubes électriques et de tubes chauffants [92]. 


L’oxyde de magnésium présente également des propriétés optiques intéressantes. MgO 
pur ne présente aucune absorption notable dans les spectres visible et Infrarouge. Par contre, il 
affiche une réflectance élevée dans le visible et le proche ultraviolet et présente un pic d'absorption 
distinctif à 589 nm, Il est transparent dans l'ultraviolet (UV), le visible et l'infrarouge (IR) lointain, 
avec une transmittance élevée dans le visible (environ 90%). Son indice de réfraction d'environ 


1.74 lui confère une faible réflectance dans le visible (environ 4%), qui augmente dans l'UV et l'IR. 
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IL. 3. 3. 4. Application de l’oxyde de magnésium 


L’oxyde de magnésium se distingue par ses applications diversifiées à travers plusieurs 
domaines en raison de ses propriétés physiques et chimiques exceptionnelles. Son exploitation 
s'étend de l'industrie à la médecine, en passant par la cosmétique et l'agriculture sans oublier 


l’optoélectronique. 


Dans l'industrie, en tant que matériau réfractaire, Mg trouve sa place dans la construction 
de fours et de revêtements nécessitant une résistance aux températures élevées 1l est utilisé dans la 
fabrication de panneaux en céramiques, de sols. Sa contribution à la fabrication du ciment en 
renforce la solidité et la longévité, tandis que son aptitude à neutraliser les acides en fait un allié 
précieux pour la purification de l'eau et des gaz. Dans l'industrie du verre, il agit en tant qu'agent 


stabilisant, améliorant la durabilité et la résistance chimique des produits finis. 


Le domaine médical s'approprie également les propriétés uniques du M£gO. En tant 
qu'antiacide, il combat l'acidité gastrique et soulage les maux d'estomac et son utilisation comme 
complément alimentaire comble les carences en magnésium, un élément vital pour l'organisme, 
soulignant aussi son rôle dans le domaine pharmaceutique. Ses propriétés exfoliantes permettent 
d'éliminer les cellules mortes de la peau, tandis que son pouvoir absorbant régule l'excès de sébum 


dans les produits cosmétiques. 


En agriculture, il enrichit le sol en magnésium, ajustant le pH et fournissant du magnésium 
aux plantes et il permet de précipiter les métaux lourds et de les éliminer. Il agit comme retardateur 


de flamme dans les matériaux de construction (contribution à la protection contre les incendies) 


Ses propriétés électroniques en font un isolant électrique utilisé dans le domaine de 
l'électronique. De plus, ses caractéristiques optiques, telles que sa transparence dans le spectre 


visible, le rend précieux dans la fabrication de composants optiques 


IL. 4. Dopage de M£O par des métaux de transition ou de terre rare 


Dans l'étude menée par Kaltoum Klaa et al. [85] les propriétés structurales, électroniques 
et magnétiques d'un isolant à large bande interdite MgO sont explorées en substituant par Mg par 
des métaux de transition 3d tels que Fe et Ni. L'analyse des structures de bandes électroniques 
révèle un caractère isolant pour M£gO dopé au Ni, tandis que celui dopé au Fe démontre un caractère 


semi-conducteur et semi-métallique. L'analyse des densités d’états électroniques suggère que 
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l'interaction de super-échange p-d confère au matériau un caractère ferromagnétique en raison des 
impuretés de métaux de transition 3d et d'un atome d'oxygène ainsi que l'hybridation p-d observée 


au sommet du bord de la bande de valence. 


Le palladium, un élément aux caractéristiques multiples, est largement utilisé dans divers 
domaines et a de nombreuses applications. Il joue un rôle clé en catalyse. Il est utilisé pour la 
fabrication de prothèses dentaires durables. Dans le domaine de l'électronique, c’est un composant 
essentiel de dispositifs électroniques et de connecteurs et dans les piles à combustible c’est un bon 
catalyseur pour générer de l'énergie électrique de manière propre... . Bien que les publications sur 
le palladium soient relativement limitées, l’étude réalisée par Jun bo Zhong et al. [93] du dopage 
du ZnO par Pd dans le a montré que le dopage augmente l'absorbance du ZnO dans la région visible 


et améliore le taux de séparation des charges photo induites. 


Cependant, J. Beloëevié-Cavor et al. [94] en utilisant des calculs de premiers principes ont 
étudié les propriétés d'intérêt pour l'application photo catalytique, à savoir la bande interdite et le 
coefficient d'absorption, de TiO: rutile dopé de manière substitutionnelle avec divers métaux de 
transition 5d. Dans cette recherche, ils ont testé différentes tailles de super cellules et diverses 
approximations du potentiel d'échange-corrélation. Les résultats calculés indiquent que tous les 
systèmes étudiés présentent un écart d'énergie direct au point let à l'exception de Ta, ils conduisent 
à un rétrécissement de la bande interdite. De plus, ils ont constaté qu’une super cellule (2X2X1) ne 
convient pas pour la description des systèmes étudiés et que pour obtenir des résultats satisfaisants, 
une super cellule (2x2x2) doit être utilisée. La tendance globale des changements de la bande 
interdite reste cohérente avec la DOS obtenu à la fois par les calculs mBJ et hybrides. En tenant 
compte de tous les critères d’une photo catalyseur efficace les chercheurs en concluent que parmi 
les métaux de transition étudiés, le Pt et l'Ir sont les meilleurs candidats pour améliorer les 


propriétés du TiO rutile. 


Une autre étude faite par Yifei Chen et al. [91] Porte sur le mono dopage du Cu ainsi que 
l’effet des lacunes (Cu, Mg ou O) dans M£g0O. Les résultats indiquent que le mono dopage de Cu 
dans M£0 préfère un état fondamental ferromagnétique et un comportement semi-métallique, avec 
un moment magnétique de 1 18 par dopant Cu. Ce moment magnétique provient principalement 
des électrons Cu 3d et O 2p polarisés en spin et qui peut être expliqué en termes de mécanisme 


d'hybridation p — d. Il a également été constaté que la lacune en Mg est propice à la stabilité du 
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ferromagnétisme semi-métallique, tandis que la lacune en O perturbe l’état fondamental 


ferromagnétique dans M£0O dopé au Cu. 


Guoyian Li [95] a également réalisé une étude sur le dopage de ZnO par un métal de 
transition de la colonne 5d (Au). Une augmentation des paramètres de réseau suite au dopage est 
mise en évidence quel que soit la position de l’atome d’or, et que l'énergie de formation de Au en 
position interstitielle est inférieure à celle en position de substitution, ce qui indique une forte 
solubilité de Au en position interstitielle. Les propriétés magnétiques, sont améliorées dans le 
modèle de substitution, tandis qu'en position interstitielle, le matériau demeure paramagnétique. 
La bande interdite diminue et le matériau ZnO dopé à l’or se transforme en semi-conducteur 
dégénéré de type n en position interstitielle, tandis que pour le modèle de substitution, la bande 


interdite diminue aussi mais le matériau ZnO dopé avec Au devient un semi-conducteur dégénéré 
de type p. 


Dans une autre étude menée par Korake et al. [96] des nanotiges de ZnO dopées à l'Eu 
(0.1%, 0.2%, 0.5%) sont préparées par une méthode assistée par micro-ondes. Les résultats ont 
confirmé l'incorporation d'Eu dans la matrice ZnO. Un élargissement et un déplacement des pics 
par DRX, vers des angles inférieur. L’énergie de la bande interdite directe des échantillons, 


diminue avec l’augmentation de la concentration du dopant Eu. 


IL. 5. Objectif générale de cette thèse 


Les semi-conducteurs à large bande interdite ont reçus une attention particulière en raison 
de leurs excellentes applications dans les domaines de l’optoélectronique, tels que les détecteurs 
UV, les diodes laser et les applications électromagnétiques. Parmi ces matériaux, l’oxyde non 
magnétique M£O est largement utilisé dans diverses applications telles que les détecteurs UV- 
Blue, et principalement dans l’électronique de spin. Il est donc impératif de réduire sa bande 
interdite et d’introduire un moment magnétique pour optimiser son efficacité dans les applications 


de l’optoélectronique et le magnéto-électronique. 


L'objectif principal de cette thèse est de se concentrer sur l'utilisation des métaux de 
transition des séries 3d, 4d et 5d en remplacement du cation Mg dans la matrice hôte MgO. Cette 


approche vise à identifier des combinaisons de métaux dopants susceptibles d'améliorer 
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significativement les propriétés électroniques magnétiques et optiques de M£O, ouvrant ainsi la 


voie à de nouvelles applications dans des domaines technologiques avancés. 


Le travail consiste à étudier et à analyser par simulations, les propriétés structurales, 
électroniques, magnétique et optiques de MgO dans sa phase rock salt à la fois pur et dopé avec 
différents métaux de transitions tels que le Nickel (Ni), le palladium (Pd), le platine (Pt), ainsi que 
d’autre métaux de transitions en faisant varier la concentration des dopants tels que le cuivre (Cu), 
l’argent (Ag) et l’or (Au). Par ailleurs, l’effet du co-dopage avec une terre rare telles que l’europium 
(Eu) et les métaux nobles (Cu, Ag, Au) est examiné afin de comprendre et d’exploiter leurs vastes 
applications potentielles. Pour ce faire, la méthode FP-LAPW est employée pour étudier ces 


systèmes. 
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CHAPITRE III 


Simulation des propriétés structurales, électroniques, magnétiques et 
optiques de la matrice hôte MgO et des composés X0.12:Mg0.8750 (X = 
Ni, Pd et Pt) 


La première partie de ce chapitre est une étude comparative des propriétés structurales, élastiques, électroniques, magnétiques et 
optiques et de l'oxyde MgO en utilisant différentes approximations alors que la seconde partie, traite l’ oxyde l'effet du dopage de 
MgO avec des atomes de métaux de transition 3d-Ni, 4d-Pd et 5d-Pt, dans la structure rock salt sur la structure électronique, les 


propriétés magnétiques et optiques en relation avec le numéro atomique des dopants.. 


« Toutes les théories sont légitimes, peu importe. Ce qui compte, c'est ce que vous en faites » 


Jorge Luis Borges 
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III. 1. Etude de la matrice hôte de l’oxyde MgO 


III. 1.1. Introduction 


Les calculs ab initio pour les semi-conducteurs à large bande interdite, tel que MgO, 
revêtent une importance significative dans la recherche en physique des matériaux et l'ingénierie 
électronique. Ces calculs offrent une compréhension approfondie des propriétés électroniques, y 
compris la largeur de la bande interdite et les niveaux d'énergie des électrons De plus, les calculs 
ab initio fournissent des informations sur les interactions atomiques, les forces électroniques et les 
propriétés structurales, contribuant ainsi à la compréhension de la stabilité cristalline de MgO. Ces 
calculs ab initio ouvrent la voie à une utilisation plus efficace de MgO dans diverses applications 


électroniques. 


L’oxyde de magnésium (M£gO) connu sous le nom de magnésie est un matériau céramique 
appartenant à la famille des oxydes métalliques alcalinoterreux [1]. Sa structure cristallographique 
dans les conditions ambiantes est cubique à face centrée (CFC), et sa surface est considérée comme 
basique (TH2014PEST1082) [2]. La maille conventionnelle de l’oxyde MgO contient 8 atomes [3]. 


C’est un semi-conducteur de large gap. 


En raison de ces propriétés remarquables, M£O à fait l'objet de nombreuses recherches 
tant expérimentales que théoriques pour évaluer sa pertinence dans divers domaines, notamment 
les applications électroniques. Sa polyvalence permet son utilisation dans des secteurs variés tels 
que les piles solaires et la géophysique. L'étude du dopage de M£O avec des métaux de transition 
constitue une étape cruciale de ces recherches, qui requiert souvent l'utilisation des techniques 
expérimentales et théoriques. Cette approche est essentielle pour comprendre les changements de 
ses propriétés induites par le dopage et optimiser les conditions de ce processus. L'objectif associé 
à ces travaux est de concevoir des matériaux dotés de caractéristiques particulières afin de répondre 


aux exigences spécifiques des diverses applications technologiques. 


IL. 1. 2. Considérations générales 


Le tableau (LIL. 1) montre les tests effectués pour choisir la valeur optimale de Rur X Kmax. 
Dans ces tests, le même nombre de K points, la même valeur de Rur pour toutes les structures et 
des paramètres du réseau fixes sont utilisés pour chaque calcul. À chaque itération du cycle d'auto- 


cohérence (SCF), l'énergie totale est enregistrée pour chaque exécution et l'objectif cherché est de 
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déterminer l'énergie minimale. À partir des résultats du tableau (TEL. 1) la valeur optimale constatée 


de Rur X Kmax est de 7.5 pour la LDA, et de 7 pour la GGA. 


Pour un système périodique, les intégrales dans l'espace réel sur le système (infiniment 
étendu) sont remplacées par des intégrales sur la première zone de Brillouin (finie) dans l'espace 
réciproque, conformément au théorème de Bloch. Dans le code Wien2K, ces intégrales sont 
calculées en sommant les valeurs de fonction de l'intégrante (par exemple : la densité de charge) 
en un nombre fini de points dans la zone de Brillouin, appelé maillage à K points [4]. Le choix d'un 
maillage suffisamment dense de points d'intégration est crucial pour la convergence des résultats, 


et constitue donc l'un des objectifs majeurs lors de la réalisation de tests de convergence. 


Des tests sont effectués avec les méthodes GGA-WC, GGA-PBEsol et LDA pour 
sélectionner les meilleurs K points, comme indiqué dans le tableau (HI. 2), en cours d'exécution 
(SCF) et en recherchant l’énergie minimale. Il est clair que pour la GGA les meilleurs K points, 
sont au nombre de 1000 avec une matrice de (10 x 10 x 10), tandis que pour la LDA, ils sont au 


nombre de 900 points avec une matrice de (9 x 9 x 9). 


Tableau III. 1. Test de convergence (R mr X Kmax) par les approximations LDA et GGA. 


LDA GGA-WC GGA-PBEsol 
No | RMTXKmax | Etot No | R MTXKmax Etot No | RMTXKmax | Etot 
Il à -548.8353 | 1 |5 -550.8110 |1 à -550.0776 
2 6 -548.9197 12 16 -550.8130 |2 6 -550.0792 
: 6.5 -548.9305 13 16.5 -550.8232 |3 6.5 -550.0808 
+ -548.9330 

-550.8238 , -550.0860 
6 8 -548.9343 16 |8 -550.8236 | 6 8 -550.0832 


Tableau IIT. 2. Test de convergence (K points) par les approximations LDA et GGA. 


LDA GGA-WC GGA-PBEsol 

No | K point Eiot No | K point Eiot No | K point Eiot 

il 600 -548.8215 l 600 -550.7960 | 1 600 -550.0795 
2 700 -548.8220 |2  |700 -550.7980 |2 | 700 -550.0712 
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3 800 -548.8282 


800 -550.8002 | 3 800 -550.0745 


4 900 -550.8019 | 4 900 -550.0839 


-548.8316 
6 1200 -548.8334 


-550.8021 -550.0886 


IIL. 1. 3. Détails de calculs 


Dans ce travail, les calculs sont réalisés en utilisant le code Wien2K. Le potentiel 


d’échange et de corrélation est traité dans le cadre des approximations suivantes : 


Ÿ  L’approximation du gradient généralisé (GGA : Generalized Gradient Approximation) 
paramétrées par Perdew, Burke et Ernzerhof, appelé (GGA-PBE). 

Ÿ”_ L’approximation du gradient généralisé : GGA : paramétrées par Wu-Cohen (GGA-WC) 

Ÿ L’approximation du gradient généralisé : révisée de Perdew, Burke, Ernzerhof du sol solide 
(GGA-PBEsol). 

Ÿ” L’approximation de la densité locale (Local Density Approximation) 

Les calculs des propriétés M£0O sont effectués dans la phase rock salt. Les états de valence 
pour les éléments utilisés sont les suivant : Mg [Ne] (3s?), O [He] (2s? 2p*). Au début de cette étude, 
une série de tests est menée pour déterminer les paramètres adéquats, notamment le choix de Rur 
X Kmax et le nombre de K points. Pour le calcul de ces structures un paramètre de coupure Rmr x 
Kmax = 7 est choisi, des rayons muffin-tin (Run de 1.5 et 1.6 a. u. Bohr. Un ensemble de 1000 K 
points spéciaux est utilisé [5]. Les fonctions de base, la densité de charge et le potentiel sont étendus 
à l'intérieur des sphères en forme de muffin en combinaison avec des fonctions harmoniques 
sphériques avec une coupure [max = 10 et en série de Fourier dans la région interstitielle. Le 
processus d'itération est répété jusqu'à ce que l'énergie totale calculée du cristal converge vers 


moins de 10 Ry assurant ainsi la stabilité et la précision des résultats obtenus 
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III. 1. 4. Résultats et discussions 
III. 1. 4. 1. Propriétés structurales 


L’oxyde de magnésium appartient au groupe spatial Fm3m N° 225 et de constante de 
réseau expérimentale égale à 4.213 À (Fig. IIL. 1) [6]. Les positions des atomes de Mg et O se 
trouvent respectivement sur les sites Wyckoff 4a (0, 0, 0) et 4b (0.5, 0.5, 0.5). 


Figure IIL. 1. Structure de M£O dans la phase rock salt. 


Tout d’abord, les paramètres de réseau de ce composé sont optimisés, en minimisant 
l'énergie totale Erx à l'aide de différentes approximations, en fonction du volume. Cette 
minimisation est réalisée en utilisant l’équation d’état de Murnaghan [7] donnée par la relation 
suivante : 

E(V) = Eo (V) + __— |B (1 _— +) + os _— 1| (HL. 1) 
B'(B' — 1) V V 
Avec : 
E, = énergie totale 
Vo = volume à l’équilibre 
B’ — dérivée première du B par rapport à la pression 
B = coefficient de compressibilité donné par : 


02E 


B=Vy 


(HILL. 2) 


A partir de cette optimisation, il est possible de déterminer différentes caractéristiques de MgO à 
l'équilibre, dont le volume optimal V, les paramètres de maille, le coefficient de compressibilité B, 
B’ la dérivé première du coefficient de compressibilité par rapport à la pression et l’énergie totale 


E. 
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La figure (IL. 2), illustre l'énergie calculée du composé hôte MgO en fonction du volume 
dans la structure rock salt en utilisant les approximations GGA : PBE, WC, PBEsol et LDA, 
respectivement. Les valeurs optimisées obtenues par les différentes approximations sont 
répertoriées dans le tableau (IIL 3), où elles sont comparées aux résultats théoriques et 


expérimentaux. Ceci sert comme test de validation de la méthodologie. 


Dans le cas du composé binaire MgO, la constante de réseau calculée avec 
l’approximation GGA-PBEsol (a = 4.220 À) et GGA-WC (a = 4.222 À) montrent un excellent 
accord avec la valeur expérimentale (a = 4.213 À) [7]. La différence ne dépasse pas 1 %, ce qui est 
notablement meilleure que celle calculée avec l'approximation GGA-PBE utilisée par Chen et al. 
[8]. Pour la LDA, le paramètre de réseau calculé est légèrement inférieur aux valeurs 
expérimentales, mais sera proche d'une valeur de 4.206 À calculée par Taib et al. en adoptant 


l’approche LDA [91]. 


D’après les résultats des GGA, le paramètre de réseau de MgO est surestimé de 4.64 %, 
0.01 % et 0.77 % avec PBE, WC et PBEsol, respectivement par rapport aux autres valeurs 
expérimentales. En revanche, dans le cas de la LDA, le paramètre de réseau sous-estime la valeur 
expérimentale de 4.65 %. Les approximations GGA-PBEsol et GGA-WC sont mieux adaptées pour 


le calcul des paramètres structuraux. 


soi Mgo] 
-551.415 + GGA-PBE GGA-WC] 


-550.806 - 


-551.420 + 
-550.812 - 


Energie (Ry) 
Energie (Ry) 


-551.425 + -550.818 + 


-550.824 - 
-551.430 + 


T L T T T T T T 
110 115 120 125 130 135 140 145 150 155 


T T T T T T T T 
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3 
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Figure IIL. 2. Variation de l’énergie totale du composé hôte MgO en fonction du volume dans la 
structure rock salt par différentes approximations. 

Tableau IIL. 3. Paramètre du réseau ao (À), module de compressibilité Bo (GPa) et sa dérivé B° de 

M£g0. 


Présent travail 


GGA-PBE 42594 148.8113 4.2951 
GGA-WC 42221 157.4312 44202 
GGA-PBEsol 4.2207 156.5702 4.1588 
LDA 4.1665 172.7289 4.3631 
Valeur expérimentale 4.213 [7, 10] 160 

Autres valeurs théoriques 4.168 (LDA) [11] 172.7 [13] 


4.237 (GGA-PW91)[12] 172.72[11] 4.36[13] 
4.166 (LDA) [13] 

4.206 [9] 

4.286 (GGA-PBE) [8] 

4.239 (GGA-PBE) [14] 


IL. 1. 3. 4. Propriétés élastiques 


Compte tenu de la symétrie cubique de l’oxyde MgO, trois constantes élastiques (C11, C12 
et C4) suffisent pour décrire les propriétés élastiques du composé. Les propriétés élastiques sont 
calculées à une pression nulle en utilisant les différentes approximations GGA et LDA. la méthode 
utilisée est celle, récemment développée par M Jamal et intégrée dans le package Wien2k [15]. Les 


résultats sont présentés dans le tableau (HIT. 4). Elles sont comparées à d'autres valeurs théoriques 
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trouvées dans la littérature [16-18]. D'après les résultats de M£gO pur, les valeurs des constantes 
Cu et C4 sont relativement élevées, tandis que la valeur de C12 est plus basse ; avec les deux 
approximations. Cela suggère que le composé est élastiquement stable avec une résistance à la 
compression unidirectionnelle [19]. L'anisotropie élastique A est aussi calculée et les valeurs 
indiquent que M£gO est un matériau anisotropie (A > 1) [20]. Cela signifie que ses propriétés 
élastiques varient selon la direction dans laquelle elles sont mesurées. Selon les calculs du rapport 
B/G l’oxyde MgO est un matériau fragile [21]. De plus, en considérant les modules de cisaillement 


M£gO est un matériau ionique et stable mécaniquement. 


Tableau IIL. 4. Constantes élastiques (C11, C12, C44), module de compressibilité B, module de 


Young E, module de cisaillement G, coefficient de poisson v et paramètre d’anisotropie A. 


GGA-PBE 291.13 87.87 184.14 155.62 332.03 145.06 0.14 1.81 1.07 
GGA-WC 295.50 88.79 147.28 157.70 301.85 127.79 0.18 1.42 1.23 
GGA-PBEsol 295.83 88.99 151.25 157.94 305.81 129.88 0.17 1.46 1.21 
LDA 334.20 96.36 152.07 175.64 327.72 137.81 0.18 127 1.27 
EXP 294.0 [22] 93.0 [22] 155.022] 163.1 121.8 [2] 
[16] 

Autre calculs 297.9118]  95.8[18] 154.4 T18] 

291.0[15] 92.01[15] 139.0 [2] 

338 [16] 91.0 [2] 118 [16] 


IIL. 1. 4. 2. Propriétés électroniques 
IL 1. 4. 2. 1. Structures de bandes (SB) 


En utilisant le code Wien2K, l'exploration de la structure de bandes et le calcul de la bande 
interdite pour M£gO sont effectués de manière facile. Le calcul de la largeur de la bande interdite 
est très important car 1l permet de déterminer si le composé est métallique, semi-conducteur ou 


isolant. En parcourant les "spaghettis", il devient très accessible de calculer la bande interdite. 
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La figure IIL. 3, présente les structures de bandes de MgO dans la phase rock salt 
calculées par les approximations LDA et GGA combinée avec la correction TB-mBJ. Les valeurs 
des bandes interdites (Tableaux IIT. 5) sont de 4.970 eV par LDA, 4.429 eV par GGA-PBE, 4.611 
eV par WC et 4.588 eV par GGA-PBEsol, correspondant à une transition directe au point l'. Quand 
la correction TB-mBJ est appliquée, les valeurs se rapprochent davantage des valeurs 
expérimentales, à savoir 6.876 eV avec GG-PBE, 7.155 eV avec GGA-WC, 7.147 avec GGA- 
PBEsol, et enfin 7.550 avec la LDA. Ces résultats sont en bon accord avec ceux de H. Baltache et 
al. [16]. MgO est un semi-conducteur à grand gap direct avec le minimum de la bande de 
conduction et le maximum de la bande de valence en point [de la zone de Brillouin, selon les 
calculs utilisant les approches de GGA et LDA. Cependant, lorsque la correction TB-mBJ est 
appliquée, MgO se comporte comme un isolant avec une large bande interdite. Dans des conditions 
normales, M£0O possède à la fois les propriétés d'un isolant avec une large bande interdite et celles 


d'un semi-conducteur, avec une grande largeur de bande de valence. 


GGA-PREso| (GGA-WC| 
1 . 10 
al 
6! 
al 
1 


e (eV) 


Energie (eV) 


= & & 2 © e © = æ 
Eu 
Energie (eV) 


Energie (eV) 


& à 2 


W  L T X WKk W OL T X WK 


MBJ-PBEsoI] 


Energie (eV) 

CE 
Es 
© 
œ 
mm 

Y 
El 
S 
= 
5 
2 

a 
5 


SO & à + 
Energie (eV) 


Energie (eV) 


de JE. 16 
Energie (eV) 


SO & & E © 


Figure IL. 3. Structures de bandes de M£O, dans la phase rock salt, calculées par les différentes 


approximations. 
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Tableau IIT. 5. Valeurs des gaps en (eV) pour M£O. 


GGA-PBE 4.429 6.876 
GGA-WC 4.611 7.155 
Direct (T°-T) 
GGA-PBEsol 4.588 7.147 
LDA 4.970 7.550 
Exp 7.8 [23] 
7.1 [24] 
Autre calculs 4.97 [13] 4.06 [23] 
4.83 [25] 7.28 [25] 
4.49 [26] 4.90 [12] 
4.61 [23] 


IIL. 1. 4. 2. 2. Densités d’états électroniques (DOS) 


La densité d'état d'un système décrit le nombre d'états à chaque niveau d'énergie 
susceptibles d'être occupés. La figure III. 4 montre les densités d’états totales et partielles de MgO, 
dans la phase rock salt, calculées par les approximations LDA et GGA, ainsi que la correction TB- 
mBJ. Les courbes de densité d'état sont tracées dans une plage d'énergie allant de -8 eV à 14 eV. 
La densité d'états est caractérisée par deux régions essentielles : la bande de valence et la bande de 
conduction. En examinant la région [-5 eV à O0 eV] au niveau de la bande de valence, celle-ci est 
essentiellement dominée par l'état O-2p. L'autre région de la bande de conduction inférieure, s'étend 
d'environ 4.6 eV à 15 eV, comprend les états O-2p ainsi que les états 3s et 2p de l'atome de Meg. 
Ces observations concordent favorablement avec les résultats des structures de bandes 


électroniques. 
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Figure IIL. 4. Densités d’états totales et partielles de MgO, calculées par les différentes 


approximations. 
III. 1. 4.3. Propriétés optiques 


La réponse optique d'un matériau dans un milieu peut être appréhendée en étudiant ses 
propriétés optiques. Dans le cas des oxydes à large gap, ces matériaux, couvrent une gamme de 
longueurs d'onde allant de l'ultraviolet au visible, ils peuvent être traités comme un milieu continu. 
Lorsqu'ils sont excités optiquement, ces matériaux ont la capacité de générer des électrons et/ou 


des trous. 


La fonction diélectrique complexe e(w) est une propriété optique cruciale qui relie le 
processus physique de transition électronique à la structure de la bande d'énergie. Elle reflète la 
réponse macroscopique linéaire du matériau face au rayonnement électromagnétique. La fonction 
diélectrique est donnée par £(@) — £1(@) + 182(@), où e2() représente la partie imaginaire qui 
caractérise l'absorption optique. D'autre part, la partie réelle £1(), obtenue à partir de la relation 
de Kramers-Kronig [27] en se basant sur la partie imaginaire, offre des informations sur la 
dispersion du rayonnement électromagnétique sur une plage de fréquences w. D'autres paramètres 
optiques tels que la réflectivité, l’indice de réfraction, le coefficient d’extinctions et le coefficient 
d'absorption peuvent être calculés à partir des composantes réelles et imaginaires de la fonction 


diélectrique. 


Pour acquérir une compréhension complète des propriétés optiques de MgO non dopé, 


divers paramètres optiques sont calculés et les courbes correspondantes sont générés (figure III. 5 
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et figure IIL 6), en utilisant l’approximation GGA-WC et WC-mBJ dans la plage d'énergie 
photonique de 0 à 30 eV. 


La figure IIL 5. présente les profils de la fonction diélectrique réelle de l’oxyde pur. Une 
distinction notable, surtout aux basses énergies, est observée entre les profils avec le changement 
de l’approximation. La constante diélectrique statique, £1(0), fait référence à la valeur de la partie 
réelle de la fonction diélectrique à énergie nulle. Des pics nets prononcés, aux valeurs positives, 
indiquant des caractéristiques semi-conductrices, sont observés du côté des basses énergies. Les 
valeurs de £&1 (0) sont 3.17 et 2.40 en utilisant la GGA-WC et WC-mBJ respectivement. Les profils 
1 (©) déduits des deux approximations présentent une tendance globale similaire, se déplaçant vers 
les énergies plus élevées. À une énergie photonique d'environ 4.9, 8.23 et 10.46 eV &1(o) atteint 
ses valeurs maximales avec la GGA-WC, correspondant à des valeurs diélectriques de 4.34, 5.16 
et 4.43. Avec la WC-mBJ, les pics sont situés à une énergie photonique d'environ 7.44, 10.95 et 
13.15 eV avec des valeurs de 3.40, 4.08 et 3.71. Par la suite, £1(©) subit une forte baisse et devient 
négative après 15.36 eV avec la GGA-WC et 19.03 eV avec WC-mBJ. Les valeurs négatives de 
e1(o) indiquent une réflectivité élevée et une excitation plasmique dans la plage d'énergie de 16 à 
24 eV, suggérant un comportement métallique. De plus, à mesure que l’énergie des photons 
augmente, £1(@) montre une tendance à devenir plus stables. Ce comportement suggère une 
limitation ou une diminution de l'interaction du matériau avec les photons incidents à haute énergie. 
Ces résultats sont en accord avec les calculs théoriques et les valeurs expérimentales de £1(@) [25, 


28, 29]. 


La partie imaginaire de la fonction diélectrique fournit des informations spectrales 
précieuses sur la réaction des matériaux au rayonnement électromagnétique incident. De plus, cela 
illustre le phénomène de décalage inter bande, qui résulte de l'interaction entre les électrons dans 
la bande de valence occupée et le rayonnement absorbé dans le matériau en dessous du niveau 
d'énergie de Fermi. Contrairement au comportement de la constante diélectrique réelle, la constante 
diélectrique imaginaire de MgO présente des valeurs positives élevées, suivies d'une forte 
diminution et de plusieurs pics de transition se produisant à des vitesses différentes dans la région 
des basses énergies. Les pics caractéristiques du M£O pur se situent entre 8 eV et 14 eV. L'origine 
de ces pics, si l'on se réfère aux courbes TDOS et PDOS du M£O pur, résulte des transitions 


électroniques impliquant principalement les orbitales 3s et 2p de Mg vers l'orbitale 2p de l’oxygène. 
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Le pic autour de 10 eV provient des transitions entre les états p-O occupés au sommet de la bande 


de valence. 


La figure IIL. 5. représente la partie imaginaire de la fonction diélectrique, calculée à 
l'aide des approximations GGA-WC et WC-mBJ. Les courbes de la partie imaginaire, de M£gO pur, 
montrent des seuils de transition égaux respectivement à 4.66 eV (GGA-WC) et 7.22 eV (WC- 
mBJ). Ces valeurs sont proches des valeurs des bandes interdites. Des pics situés à 8, 6.9, 10.95 et 
15.24 eV (GGA-WC) et 11.30, 13.51 et 17.56 eV (WC-mBJ) sont également observés. Ces pics 
résultent des transitions électroniques entre les états occupés et inoccupés. Les pics autour de 8.0 
et 11.0 eV sont principalement dus à des transitions directes entre le maximum de la bande de 
valence et le minimum de la bande de conduction au-dessus du niveau d'énergie de Fermi. Après 
23 eV, aucune transition n’est observée, similaire à un comportement transparent. Après 20 eV, les 
résultats du mBJ sont comparables à ceux de la GGA, néanmoins avec quelques différences, telles 


que la diminution du nombre de pics. 


Le coefficient d'absorption calculé de M£gO non dopé en fonction de l'énergie est illustré 
à la figure IIL. 6 (a). MgO absorbe la lumière dans la gamme UV selon les approximations GGA- 
WC et WC-mBJ. La majeure partie de la lumière absorbée traverse le matériau ou est utilisée lors 
des transitions électroniques. Ainsi, dans cette gamme d’énergie (10-20 eV), contrairement aux 
basses énergies, MgO se comporte comme un matériau transparent, alors qu’il est opaque aux 
basses énergies, ce qui montre que la variation de l’énergie des photons incidents conduit à des 
comportements optiques différents dans ce composé. Les transitions électroniques de la bande de 
valence O-2p vers la bande de conduction Mg-3s démontrent ce comportement. L'intensité du 
coefficient d'absorption est nulle pour des énergies égales à 4.66 eV (GGA-WC) et 7.22 eV (WC- 
mBJ). Ces données, simulées avec les deux approximations, s'accordent bien avec les données de 
propriétés électroniques ; les seuils d'absorption de M£gO sont situés à proximité de leurs gaps 
énergétiques respectifs. Il est également souligné que la GGA entraîne un décalage vers le rouge, 


tandis que dans la zone UV, le comportement d'absorption des deux approximations est similaire. 


La figure. IIL 6 (b). représente la réflectivité R (©) de MgO en fonction de l'énergie. 
Celle-ci caractérise la puissance réfléchie à l'interface du solide. La réflectivité maximale est 


d'environ 46.7 % à 18.32 eV avec GGA-WC et 37.27 % à 20.93 eV avec WC-mBJ, où la valeur 
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minimale £1(w) a été estimée. L’absorption est moindre et les composés montrent la transition inter 


bande des photons dans cette région. 


L'indice de réfraction n (©) renseigne sur la propagation d'une onde électromagnétique 
dans un solide et sa connaissance est fondamentale pour caractériser les matériaux utilisés dans les 
dispositifs optoélectroniques [52]. La figure IIL. 6 (ec) montre la courbe n (w) calculée entre 0 et 30 
eV. L'indice de réfraction statique est égal à 1.77 avec WC-mBJ et avec GGA-WC sa valeur est de 
1.55, indiquant un comportement semi-conducteur. Dans la région de 6 à 9 eV, où la transition 
d'électrons se produit et E £1 (@) est maximale, il y a une augmentation relative de l'indice de 
réfraction. Mais avec l'augmentation de l'énergie du photon incident, l'amplitude de l'indice de 
réfraction diminue lentement. Après 12 eV, l'indice de réfraction calculé avec les deux 


approximations est inférieur à 1. 


La partie imaginaire du coefficient de réfraction définit le coefficient d'extinction k (œ), 
qui fait référence à la réduction de l'amplitude de la lumière lors de sa propagation à travers le 
matériau. Comme le montrent les fig. IL 5 et 6 (d), la réponse de la matière à la lumière incidente, 
représentée par la courbe la courbe £1 (©) se produit dans la plage de 10 eV pour GGA-WC et 12 
eV pour WC-mBJ. Les maximums de k (w) dans la région ultraviolette, apparaissent dans MgO 
pur autour de 8.72 eV et 11.19 eV avec des amplitudes égales respectivement à 0.87 et 1.54 avec 
GGA-WC. Dans le cas de la correction WC-mBJ, les pics se situent à 13.70 eV et 11.22 eV avec 
des amplitudes égales respectivement à 1.24 et 0.71. En augmentant l'énergie lumineuse, la réponse 
de la matière diminue Cela suggère que bien qu'il y ait une absorption notable de lumière, elle n'est 


pas caractérisée par une perte excessive d'amplitude. 
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Figure IIL. 5. Fonction diélectrique complexe de MgO avec approximations GGA et mBJ. 
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Figure IIL. 6. Paramètres optiques de Mg0O par GGA-WC et WC- mBJ : a) Coefficients 


d’absorption b) réflectivité c) L'indice de réfraction et d) le coefficient d'extinction. 


IIL.2. Simulation des propriétés structurales, électroniques, 
magnétiques et optiques des composés X0.12:Mg0.8750 (X = Ni, Pd et Pt) 


III. 2. 1. Introduction 


Un matériau qui présente à la fois des propriétés ferromagnétiques et à large bande 
interdite offre la perspective passionnante de combiner le stockage magnétique non volatil avec 
l’électronique conventionnelle à semi-conducteurs (isolants) dans un seul dispositif. Cet hybride 
prometteur ouvre de nouvelles voies dans le domaine de la "spintronique", une branche de la 
science et de la technologie qui explore de nouveaux dispositifs électroniques exploitant les 
propriétés de spin des électrons. Dans ce contexte, comment favoriser le magnétisme dans un semi- 
conducteur (isolant) est la clé pour résoudre ce défi. Le dopage avec les métaux de transitions se 
présente comme un outil qui permet de modifier les propriétés physiques des matériaux. Cette étude 
se focalise sur l’oxyde de magnésium M£gO dopé avec trois métaux de transition distincts (Ni, Pd, 


et Pt) chacun caractérisé par un nombre atomique Z différent. 


Les investigations portent sur les propriétés structurales, électroniques, magnétiques et 
optiques du matériau dopé. L'objectif ultime est de comprendre comment le dopage avec ces 
métaux de transition (3d, 4d et 5d) peut influencer ces propriétés, particulièrement en ce qui 
concerne la promotion du magnétisme au sein du semi-conducteur M£0O. Ces avancées pourraient 
contribuer au développement de dispositifs spintronique novateurs, fusionnant les avantages du 


stockage magnétique avec ceux de l'électronique conventionnelle à semi-conducteurs. 


IIL. 2. 2. Méthode et détails de calculs 


Les calculs sont effectués à l'aide de la méthode des ondes planes augmentées linéarisées 
(FP-LAPW), dans le cadre de la théorie fonctionnelle de la densité (DFT), en utilisant le code 
Wien2K. Pour le traitement du potentiel d'échange et de la corrélation, l'approximation du gradient 
généralisé (GGA) révisée par Perdew, Burke, Ernzerhof du solide (GGA-PBE. PBEsol) est 
adoptée. Cependant, le potentiel de Becke-Johnson modifié Tran-Blaha (TB-mBJ) a également été 
appliqué afin d'améliorer la précision des caractéristiques magnétiques et optiques, ainsi que pour 
corriger la bande interdite électronique. Le choix de la TB-mBJ est motivé par sa capacité reconnu 
à s'aligner de manière robuste avec les valeurs expérimentales, renforçant ainsi la fiabilité des 


conclusions obtenues. Dans la région des sphères muffin-tin utilisée pour assurer la convergence 
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des harmoniques sphériques avec Imax = 10, les fonctions d'onde, dans la région interstitielle sont 
étendues en ondes planes avec une énergie de coupure RmTX Kmax = 7 (où Rmr est le petit rayon 
des sphères muffin-tin. L'énergie de coupure définit la séparation entre les états de cœur et de 
valence est de -6 Ry [30]. Les rayons Muffin-Tin (Rmr) sont pris égaux à 2 u.a pour les atomes Ni, 
Pd, Pt, 1.86 et 1.72 pour Mg et O, respectivement. La configuration électronique des atomes est 
Mg [Ne] 3 s?°, O [He] 2 s° 2 p*, Ni [Ar] 3 d 4 s?, Pd [Kr] 4 d'°, Pt [Xe] 4 f*5 d°6s!. Pour simuler 
la structure cristalline des systèmes dopés MgO (Ni, Pd et Pt), une super cellule de 16 atomes (2 x 
2 x 2) est construite. Celle-ci contient 8 atomes de Mg et 8 atomes d’oxygène, sur la base d'une 
cellule unitaire conventionnelle NaCI de M£gO. Pour atteindre une concentration totale d'ions 
d'impuretés de 12.50 %, une substitution est réalisée en remplaçant un atome de (Ni, Pd et Pt) sur 
les sites Mg (c'est-à-dire Xo.125Mgo.8750). Ainsi, la configuration finale de la forme X0.125M£0.875 
O, où X représente l'élément dopant (Ni, Pd ou Pt). Cette substitution vise à reproduire les 
conditions expérimentales où cette concentration spécifique d'impuretés peut être obtenue de 


manière pratique dans les expériences. 


Pour effectuer le calcul de l'optimisation géométrique et des structures électroniques, 
l'intégration de la zone de Brillouin est effectuée en utilisant le schéma Monkhorst - Pack [5] avec 
une grille de 10 x10 x 10 Kk points pour les cellules. L'itération est poursuivie jusqu'à ce que la 
différence d'énergie totale soit inférieure à 10 * Ry, indiquant une convergence satisfaisante du 
calcul. Pour illustrer la structure cristalline, le logiciel Vesta [31]est utilisé. Vesta permet de 
visualiser de manière intuitive les configurations atomiques et les arrangements cristallins, offrant 


ainsi une représentation graphique de la géométrie optimisée résultante. 


IIL. 2. 3. Résultats et discussions 


IL. 2. 3. 1. Propriétés structurales 


Les propriétés structurales notamment les constantes de réseau optimisé (a), le module de 
compressibilité (B), sa dérivée (B°) et l’énergie minimale (E) sont calculés pour les trois composés 
étudier Nio.125Mgo8750, Pdo125Mgos8750 et Pto125Mgos750. Les structures de ces derniers sont 
présentées dans la Figure IL. 7. La variation de l’énergie totale (Ry) en fonction du volume (au) 


est ajustée en utilisant L’équation de Murnaghan (équation II. 1.). 
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Les paramètres structuraux optimisés des systèmes dopés, dans la configuration FM, (la 
phase magnétique stable, Fig. IIL. 8) sont répertoriés dans le tableau IIT. 6. Ceux-ci sont comparés 
à ceux du Mg0O pur (tableau IIL. 3). Pour le nickel, le paramètre de maille est plus petit que celui de 
M£0 pur à cause du plus petit rayon ionique de Ni (0.69 À), par rapport à Meg (0.72 À). En revanche 
pour les deux dopants Pd et Pt les valeurs sont plus grandes (Pd= 0.86 À et Pt= 0.86 À) par rapport 
M£0 pur. 


L'accroissement progressif de la constante de réseau a avec l'augmentation de Z, indiquant 
une expansion du réseau, résultant du rayon ionique plus restreint du Mg par rapport aux dopants. 
Parallèlement, Les modules de compressibilité B montrent une légère décroissance, avec la 
diminution du numéro atomique Z, suggérant une réduction de la dureté pour les systèmes dopés 
[32]. Les résultats montrent que les données simulées sont en bon accord avec les valeurs obtenues 
par d'autres calculs théoriques par rapport aux données expérimentales d'études antérieures sur 


M£0 pur telles que publiées dans les références suivantes [33-35]. 


Figure IIL. 7. Cellules cristallines des composés X0.125Mg0.8750 (X = Ni, Pd et Pt). 
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Figure IIL 8. Variation de l’énergie totale en fonction du volume des composés X0.125Mg0.8750 


(X = Ni, Pd et Pt) dans la phase (RS). 


Tableau IIT 6. Valeurs calculées de la constante de réseau a (À), module de compressibilité B 


(GPa), et sa dérivé B° des composés X0.125M 20.850 (X = Ni, Pd, et Pt). 


Nio.125Mg0.8750 4.2072 165.1479 4.2675 
Pdo.125M g0.8750 4.279 163.9584 4.8058 
Pto.125Mg0.8750 4.306 162.2916 4.2149 


IIL. 2. 3. 2. Propriétés électroniques 
IL 2. 3. 2. 1. Structures de bandes (SB) 


L'examen de la structure des bandes électroniques et de la densité partielle et totale des 
états (PDOS et TDOS) nécessite des recherches sur les propriétés électroniques des matériaux. Le 
gap énergétique est une grandeur physique essentielle utilisée dans les structures électroniques pour 
décrire la différence d'énergie entre le maximum de la bande de valence (VBM) et le minimum de 


la bande de conduction (CBM). 


Des observations antérieures (M£gO pur) indiquent que les bandes interdites des matériaux 


sont généralement sous-estimées lors de l'utilisation de l’approximation WC-GGA en raison de la 
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présence d'erreurs d'auto-interaction [36]. Bien que des résultats plus précis puissent être obtenus 
en employant d’autre méthodes coûteuses en calcul telles que la DFT avec correction d'auto- 
interaction (DET+U) ou les fonctionnelles hybrides (HSE) [37], une approche alternative consiste 
à recourir à la fonctionnelle d'échange mBJ. Cette méthode s’est avérée susceptible de fournir des 
calculs de bande interdite raisonnablement précis sans engendrer les coûts de calculs élevés. Ainsi, 
une étude est menée pour mieux comprendre l'impact du dopage sur les structures électroniques 
des systèmes X0.125M20.8750 (X=Ni, Pd, Pt) en utilisant la correction du potentiel mBJ. La figure 
IL. 9. présentent les structures de bandes polarisées en spin, dans la première Zone de Brillouin 
(BZ) de ces composés. Le niveau de Fermi Er est représenté par une ligne pointillée et est pris égal 


à 0 eV. 


Le niveau de Fermi (Er) reste inoccupé, étant positionné entre la VBM et CBM sans 
chevauchement de ces bandes, 1l ressort que tous les systèmes (Nio.125Mg0.8750 Pdo.125Mg0.8750 et 
Pto.125M£g0.8750) étudiés sont des matériaux semi-conducteurs. Pour les spins majoritaire, le 
maximum de la bande de valence (VBM) et le minimum de la bande de conduction (CBM) sont 
situés au point l', tant pour le composé Nio.123M20.8750 que pour Pdo.125Mg0.8750, conduisant à une 
bande interdite directe large de 5.729 eV et 3.89 eV respectivement. Par contre dans le cas du 
composé Pto.125Mg0.8750 elles sont situées au point X, avec un gap direct de 3.267 eV. Dans le cas 
des spins minoritaire, pour les deux composé N1o.125Mg0.8750 et Pdo.125Mgo.8750 le VBM et le CBM 
est situé au point X, indiquant une bande interdite directe. Pour le composé Pto125Mgo8750 le 
maximum de la bande de valence (VBM) et le minimum de la bande de conduction (CBM) sont 
situés au point l'. L'ampleur de la bande interdite diminue à mesure que le nombre atomique Z 
augmente, avec des valeurs de 4.92 eV, 3.31 eV et 2.02 eV pour le dopage Ni, Pd et Pt, 


respectivement. 


En outre, Il est pertinent de noter que le VBM dans les structures de bandes des trois 
composés sont situées à proximité immédiate de Er, ce qui indique leur caractéristique de semi- 
conducteurs de type p. Ces résultats sont présentés dans le tableau III. 7. Comparé au M£0O pur, 
qui présente le caractère isolant [38], l'introduction de porteurs de type trou par le dopage des 
métaux de transition conduit à une réduction des valeurs de la bande interdite et déplace le niveau 


de fermi EF vers le bas, ce qui entraîne un comportement semi-conducteur. 
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Tableau III. 7. Valeurs des gaps énergétiques pour les canaux à spin majoritaire (up) et minoritaire 


(down), à l'aide des différentes approximations des composés X0.125M2g0.7500 (X : Ni, Pd et Pt). 


Energie (eV) 


nm 


Energie (eV) 


Energie (eV) 


PBEsol-mB1} 
0 : 


Energie (eV) 


Figure IIL. 9. Structures de bandes polarisées en spin des composés X0.125Mg0.8750 (X = Ni, Pd et 


Pd 202750 


Energie (eV) 


Energie (eV) 


mm 


Up Down 
Nio.125Mg0.8750 3.087 1.332 GGA-PBEsol 
5.729 4,929 PBEsol- mBJ 
Pdo.125Mg0.8750 1.309 1.546 GGA-PBEsol 
3.891 4,312 PBEsol- mBJ 
Pto.125Mg0.8750 0.00 00.00 GGA-PBEsol 
3.267 2.002 PBEsol- mBJ 
(PBEsol-mB,} 10 PBEsol-mBJ Pd 22 MO 8750 


Pt) calculées avec l’approximation PBEsol-mBJ. 


E 
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IL. 2. 3. 2. 2. Densité d’états électroniques (DOS) 


Pour mieux élucider la nature de la structure de la bande électronique et pour analyser les 
contributions des différents atomes et orbitales, les densités d'états électroniques totales (TDOS) et 
de partielles (PDOS) par cellule unitaire, pour les systèmes Nio.125Mg0.8750 Pdo125Mgo.8750 et 


Pto.125Mg0.8750 sont calculées et discutées. 


La figure IIL. 10. représentent les densités d’états totale et partiale pour chaque composé 
où la ligne pointillée rouge représente le niveau de Fermi. Il convient de noter que tous les 
composés X0.125M 208750 présentent un caractère semi-conducteur pour les deux spins, en raison 
de l’absence de chevauchement au niveau de Fermi. La valeur de la bande interdite obtenue et la 
nature semi-conductrice sont en bon accord avec les résultats des structures de bandes. 
L'émergence de cette bande interdite dans les TDOS témoigne de la caractéristique semi- 
conductrice des composés X0.125M 208750. Les contributions les plus significatives au TDOS 
proviennent des orbitales Ni-3d (Nio125Mgo.8750), Pd-4d (Pdo125Mg08750), Pt-5d (Pto125Mgo 8750). 
Les spectres DOS sont examinés dans un intervalle d'énergie de — 7.0 à 7.0 eV. Les densités d'états 
partielles, PDOS, sont également calculées pour les atomes X (Ni (3d-eg, 3d-t2g), Pd (3deg, 3dt2g), 
Pt (3deg, 3dt2g)), Mg (35) et O (25, 2p). 


La bande interdite d'énergie dans ces systèmes trouve son origine principalement dans les 
atomes de Ni, Pd et Pt respectivement. Une analyse détaillée des PDOS révèle que les états partiels 
Ni-3d, Pd-4d et Pt-Sd, contribuent majoritairement à la formation des bandes de conduction dans 
leur composé X0.12:Mg0.8750 correspondant. Selon la théorie du champ cristallin, les cinq orbitales 
dégénérées de d-TM subissent une division en deux parties : les niveaux t2g (dx, dxy, dyz) trois fois 
dégénérés de basse énergie et les niveaux eg (d:x2-dy2, et d2) deux fois dégénérés de haute énergie 
symétriques, créant ainsi des états liants et anti-liants, ceci joue un rôle majeur dans l'induction du 


ferromagnétisme dans les composés X0.125Mg0.8750. 


Pour le composé Nio.125Mg0.8750, la bande de valence peut être divisée en trois parties 
distinctes. Dans la première plage, allant du niveau de Fermi jusqu’à -1eV, les états des-Ni dominent 
pour les spins majoritaires, tandis que pour les spins minoritaire d-t2-Ni qui sont dominants. La 
seconde plage, située entre -1eV et -2 eV est principalement occupée par les états d-t2,-Ni. Pour la 


parité entre (-1.5eV et -5 eV), les atomes de p-O dominent avec de faible contributions des autres 


105 


atomes. Dans la bande de conduction, les états d-eg poussent les états d’oxygène vers des énergies 


plus basses, créant ainsi un gap énergétique dans le nouveau composé. 


Pour le composé Pdo.125Mg0.8750 et Pto.125Mg0.8750 la bande de conduction est caractérisée 
par une dominance des états d-eg, entraînant un déplacement des états d’oxygène vers des énergies 
plus basse. En revanche, la bande de valence peut être divisée en trois parties : la première au 
niveau de fermi, est principalement occupée par les états d-eg-(Pd, Pt). Alors que la partie de - 
1.5eV à -3 eV, ce sont les états d-t2,-(Pd, Pt) qui contribuent de manière significative à la bande de 
valence. Enfin, pour la partie de -deV à -6 eV, les orbitales de p-O qui dominent avec une forte 
contribution pour Pd et faible contributions des autres atomes (Mg, et Pt), tandis qu’une 
contribution relativement faible des autres atomes est observée. Dans la bande de conduction, dans 
la région de 2 eV à 3 eV, les états d-eg poussent les états d’oxygène au basse énergie et induit le 


gap énergétique de nouveaux composés. 


Les orbitales O-2p sont repoussées au niveau de Fermi par les orbitales 3d, 4d et Sd des 
métaux de transition. Cela empêche ainsi l'hybridation, conduit à la séparation des orbitales les 
unes des autres et entraîne ainsi la création d’une bande interdite dans les deux états de spin-down 


(minoritaire) et spin-up (majoritaire). 


En résumé, l'origine de la bande interdite d'énergie dans ces systèmes est attribuée aux 
contributions majeures des états partiels des orbitales de Ni-3d, Pd-4d et Pt-5d, conformément à la 
théorie du champ cristallin qui divise les orbitales dégénérées en niveaux t2g et eg, influant ainsi 
sur la formation des bandes de conduction et induisant le ferromagnétisme dans les composés 


X0.125M 20.8750. 
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Figure IIL 10. Densité d’états électroniques total et partielle des composés X0.125M20.8750 (X = 
Ni, Pd et Pt). 


IIL. 2. 3. 3. Propriétés magnétiques 


Les atomes de métaux de transition 3d (Ni), 4d (Pd) et Sd (Pt) substitués dans l’oxyde 
alcalino-terreux M£O, créent un moment magnétique dans les super cellules Xo.125Mg0.850 (X = Ni, 
Pd, Pt). Le moment magnétique de spin total des supers cellules, y compris la contribution de 
chaque atome à la création du moment magnétique dans les composés (super cellules), sont 
regroupés dans le tableau. IL. 8. Le moment magnétique du composé hôte (M£gO non dopé) est 
nul cependant, ce moment magnétique des composée augmente suite au dopage avec des atomes 


de métaux de transition Ni, Pd et Pt. 


La configuration électronique de valence des différents atomes est la suivante, Mg [Ne] 3 
s?, O [He] 2 s° 2 p*, Ni [Ar] 3 dÿ 4 s?, Pd [Kr] 4 d!°, Pt [Xe] 4 fl# 5 d°6s!. Dans les super-cellules, 
la création du moment magnétique est attribuée aux électrons de valence. Cela signifie donc que 
les atomes de métaux de transition jouent un rôle important dans la production de magnétisme dans 
les oxydes alcalino-terreux. Le site cationique de substitution des systèmes Ni (34°), Pd (4d!°), Pt 
(5d°) a une configuration électronique de (tegf, eg?), (teg®, eg“) et (teg°, eg), respectivement, dans 


l'environnement octaédrique. 
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Tableau IIL 8 Moments magnétiques totaux, locaux et interstitiels calculés à l'aide des 


approximations GGA-PBEsol et PBEsol-mBJ, pour les composés X0.125Mg0.8750 (X = Ni, Pd, Pt). 


Nio12:Mg0.8750 PBEsol-mBJ 2.000  0.0087 0.0402  1.7946 0.03792 
GGA-PBEsol 1.995 -0.0154  0.0682 1.5554  0.02792 
Pdo.2Mg0.750 PBEsol-mBJ 2.000  (0.01927  0.09926 1.430 0.02927 
GGA-PBEsol 2.016  -0.02762 0.105908 1.235 0.13129 
Pto2sM 808750 PBEsol-mBJ 1.99  O0.02436  0.02206 1.3990  0.02016 
GGA-PBEsol 045  -0.00572  0.09489 0.264 0.05610 


IIL. 2. 3. 4. Propriétés optiques 


Les propriétés optiques d'un matériau sont étroitement liées à sa structure cristalline et à 
structure de bande électronique. La réponse linéaire des matériaux peut être décrite par la fonction 


diélectrique complexe £ (©) [38]. 


La manifestation de divers phénomènes dans un milieu exposé à la lumière trouve son 
explication à travers plusieurs paramètres optiques qui définissent les propriétés optiques de ce 
milieu, comme indiqué dans la référence [39]. Parmi ces paramètres, la fonction diélectrique, 
représentant la réponse linéaire du milieu à la lumière, se décompose en parties réelle (£1 (œ@)) et 
imaginaire (£2 (@)). Ces deux composantes contribuent à caractériser les propriétés optiques 


globales du milieu. 


La figure (NL. 11) montre les variations des parties réelles £1 (@) et imaginaires £2 (@), de 
la fonction diélectrique, en fonction de l’énergie. Les valeurs de la constante diélectrique statique 
sont 2.41, 2.68, et 2.90 pour Nio.125M20.8750, Pdo.125Mg0.8750, and Pto125Mg0.8750, respectivement. 
Une tendance positive de la variation de £1 (©) est observée jusqu'à 18 eV, indiquant un 
comportement diélectrique, Au-delà de 18 eV, £1(@) devient négative, puis augmente à nouveau à 
partir de 24.5 eV pour devenir positive. £1 (©) augmente avec le numéro atomique des dopants. La 


partie imaginaire de la fonction diélectrique £2 (@) se manifeste à un seuil d'environ 7.47 eV, 6.97 
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eV, et 5.89 eV pour Nio.125Mg0.8750, Pdo.125Mg0.8750, and Pto125M 208750, respectivement. Ceci est 
étroitement lié à la transition directe au point ” résultant de l'excitation d'électrons de la bande de 
valence vers la bande de conduction. Les valeurs maximales £2 (&) sont atteintes aux énergies 
autour 13.94 eV, 13.74 eV, et 14.27 eV, respectivement, pour Nio.125Mg0.8750, Pdo125Mg0.8750, et 


Pto.125Mg0.8750. Au-delà de ces valeurs, à mesure que l’énergie augmente, £2 (@) diminue. 


Le coefficient d'absorption optique a (w), illustré dans la figure (IL. 12), quantifie la 
capacité d'un matériau à absorber l'énergie des photons. En d'autres termes, il renseigne sur la 
fraction de lumière incidente qui est absorbée par le solide. Les composés étudiés sont caractérisés 
par une très faible absorption dans la région visible tandis qu'ils présentent une absorption 
significative dans la région ultraviolette d'énergie de 5 eV à 25 eV. À mesure que le numéro 
atomique augmente de la série 3d à 4d, puis à 5d, la région de forte absorption se décale vers des 
énergies élevées, ce qui entraîne un déplacement du bord d'absorption vers le spectre ultraviolet. 
L'extension du bord d'absorption vers l'UV est due à l'augmentation de la largeur de la bande 
interdite et à la modification de la structure de la bande électronique [40]. Dans la gamme d'énergie 
de 8 à 18 eV, la majeure partie de la lumière incidente traverse le matériau, le rendant transparent. 


A des énergies plus basses, le matériau devient opaque, absorbant la lumière incidente. 


Nio.125M96.8750 
—— Pdo.125M90.875O 
—— Pto125M96 8750 


€, (@) 


€ (@) 


0 2 4 6 8 10 12 14 16 18 20 22 24 
Energie (eV) 


Figure IIL. 11. Partie réelle (a) et partie imaginaire (b) de la fonction diélectrique des oxydes 
ternaires Nio.125Mg0.8750, Pdo.125M 208750, and Pto125Mg0.8750 en utilisant l'approximation 


PBEsol-mBJ. 
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Figure IIL. 12. Coefficient d'absorption des oxydes ternaires Nio.125Mg0.8750, Pdo.125Mg0.875O, et 


Pto.125Mg0.8750 en utilisant l'approximation PBEsol-mBJ. 


III. 2. 4. Conclusion 


La première partie de ce travail est une étude comparative des propriétés structurales, 
électroniques et élastiques de l’oxyde MgO. Les calculs sont basés sur la théorie de la fonctionnelle 
de la densité (DFT) en utilisant la méthode des ondes planes augmentées et linéarisées (FP-LAPW). 
Pour le potentiel d'échange et de corrélation, deux approximations sont utilisées : la densité locale 


et le gradient généralisé. 


Les paramètres structuraux tels que la constante de réseau (a), le module de 
compressibilité (B) et sa dérivée (B') sont calculés dans la structure NaCI (B1). Le composé est 
stable dans la phase non magnétique (paramagnétique). L'analyse et le calcul des paramètres 
élastiques, tels que le module de Young (E), le module de cisaillement (G), le rapport B/G, le 
coefficient de poisson (v), le facteur d'anisotropie (A) montrent que M£O est de nature anisotrope 
et que le composé est mécaniquement stable et fragile. Les propriétés électroniques sont calculées 
en utilisant les approximations mBJ-LDA et mBJ-GGA. Les résultats montrent que l'oxyde MgO 
est un matériau isolant avec une bande interdite directe de 7.55 eV avec mBJ-LDA en bon accord 


avec les résultats expérimentaux. Pour les propriétés optiques, un décalage vers les bases energies 
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dans la zone UV dans l'absorption. Les résultats théoriques sont en bon accord avec les travaux 


expérimentaux et théoriques trouvés dans la littérature. 


Dans la seconde partie, les oxydes alcalino-terreux M£gO, dopés avec des atomes de 
métaux de transition 3d-Ni, 4d-Pd et 5d-Pt, dans la structure rock salt sont étudiés afin d’examiner 
la structure électronique, les propriétés magnétiques et optiques. Nous avons utilisé à la fois 
l'approximation GGA-PBEsol et l’approches TB-mBJ sont utilisées pour le traitement de 
l'échange-corrélation. Les composés dopés Xo125Mg08150 (X = Ni, Pd, Pt) présentent un 
comportement semi-conducteur ferromagnétique, mis en évidence par les structures de bandes et 
les DOS. L'analyse de PDOS révèle une contribution notable des orbitales 3d, 4d, et 5d, entraînant 
la formation de bandes interdites spécifiques dans les canaux de spin majoritaires et minoritaires 
(spin up et down) autour du niveau de Fermi. Cette caractéristique souligne l'importance des 
métaux de transition (Ni, Pd, Pt) dans l'augmentation des moments magnétiques de la matrice hôte 
M£0O. En effet, les valeurs des moments magnétiques totaux pour les super cellules de 
Nio.125Mgo8750, Pdo12:Mg0.8750 et Pto12:Mg08750 sont respectivement de 2, 2, et 1.99 ag. Les 
composantes de la constante diélectrique (réelle et imaginaire), et l’absorption des composés sont 
calculées et analysées. Les composés présentent des propriétés électroniques, magnétiques et 


optiques pouvant être utilisés dans les domaines de la spintronique et l’optoélectronique. 
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CHAPITRE IV 


Etude des propriétés structurales, électroniques, magnétiques et 
optiques des composés TM;M£g1:0 (TM= Cu, Ag, et Au ; x = 
3.125 %, 6.25%, 12.5%, et 25%) 


Dans ce chapitre, les résultats de calculs de la stabilité structurale des propriétés électroniques, optiques et magnétiques des composés TM.Mg,. 
«0 (TM Cu, Ag, et Au ; x = 3125%, 6.25%, 12.5%, et 25%) sont présentés et discutés. L'étude est réalisée dans les phases ferromagnétiques 
(FM) et paramagnétique (PM), en utilisant l’approximation GGA-WEC avec la correction TB-mB] associée à la théorie de la fonctionnelle de la 
densité (DFT). 


1 A LEA x : À * : & +1 LA 4 A I 
Si l'électron nous a servi à comprendre beaucoup de choses, nous n'avons jamais bien compris l'électron lui-même. 


Louis de Broglie 


114 


Table des matières 


CHAPIERB IN 2 am nan en nt este ton en One 114 
Etude des propriétés structurales, électroniques, magnétiques et optiques des composés TM,Mg:.:0 
(TM= Cu, Ag, et Au ; x = 3.125 %, 6.25%, 12.5 %, et 25%)... 114 
IV::L: Introduction.:.::2.,.5118re8 dan ns ne na me ten rene tee free rent ele te tend 116 
IV. 2. Méthode et détails de calculs... 116 
IV. 3. Résultats et discussions... 117 
EVE PEODPIÉTÉS SEUCRUAAIES 2 Li ne ne ne ie Ce 117 
IV.3.2. Propriétés électroniques .…............................................................................... 122 
IV. 3. 2.1. Structures de bandes d’énergie (SB) ..…................................................................... 122 
IV.3.2.2. Densités d’états électroniques (DOS) .…............................................................... 129 
IV. 3.3 Propriétés magnétiques ..…............................................................................. 137 
IV. 3.3. 1. Moment magnétique .…...................................................................... 137 
IV.3.3.2. Constantes de couplage d’échange 138 
IV. 3. 4. Propriétés optiques 140 
IV;4.:. Conclusion: tasntnnininhlaon effet tt ina Rtnsinnlrn 142 
Bibliographie de chapitre LV... 144 


115 


IV. 1. Introduction 


L’oxyde de magnésium, un composé largement étudié, se trouve au cœur de nombreux 
efforts des recherches académiques considérables pour éclairer et contrôler les mécanismes de 
dopage de MgO, avec l'objectif ultime d’améliorer ses propriétés physiques et prédire de nouvelles 
applications prometteuses [1-4]. Dans ce chapitre, les résultats de calculs relatifs à la stabilité 
structurale ainsi qu'aux propriétés électroniques, optiques et magnétiques des composés TM;Mg:- 
xO sont présentés, avec les concentrations qui varient de x = 3.125%, 6.25%, 12.5%, à 25% et aussi 
les effets du dopage avec les métaux de transition des séries 3d a 4d a 5d sont également analysés. 
L’étude prend en compte les phases magnétiques, notamment les états ferromagnétiques (FM), et 
paramagnétique (PM). Les calculs sont menés en se basant sur l’approximation GGA-WC avec la 
correction TB-mBJ associée à la théorie de la fonctionnelle de la densité (DFT). Cette étude fournit 
une analyse approfondie de la structure et des propriétés des composés étudiés, mettant ainsi en 
lumière les influences subtiles du dopage et de la concentration des dopants sur les caractéristiques 


fondamentales de ces matériaux. 


IV. 2. Méthode et détails de calculs 


Dans ce travail, les calculs sont effectués avec le progiciel de simulation Wien2K [5] pour 
étudier les propriétés structurales, électroniques et magnétiques de M£gO dopé par des TM (Cu, Ag, 
et Au). La méthode d'onde plane augmentée linéarisée à plein potentiel plus orbitales locales (FP- 
LAPW + lo) [6] basée sur la théorie de la fonctionnelle de la densité (DET) est utilisée [7]. Pour le 
calcul des propriétés structurales, le potentiel d'échange-corrélation est traité par l'approximation 
de gradient généralisée proposé par Wu-Cohen (GGA-WO) [8]. De plus le potentiel d'échange de 
Becke-Johnson modifié par Tran-Blaha (TB-mBJ) est employé pour calculer les propriétés 
électroniques, offrant une estimation précise de la bande interdite des semi-conducteurs [9-11]. 
Les calculs sont basés sur trois types de super cellules (2 x 2 x 2), réseau cubique à faces centrées 
(FCC), réseau cubique centré (BCC) et réseau cubique simple (SC), appartenant aux groupes 
d'espace N° (225), N° (229), et N° (221), respectivement. La concentration obtenue est de x = 0.25, 
0.012, 0.06 et 0.03 avec le remplacement des atomes de M£ par des atomes de cuivre, d'argent et 


l’or. 
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Les rayons Muffin-Tin (MT) sont égaux à 2, 1.86 et 1.72 a. u. pour TM, Mg et O, 
respectivement. L'équation d'états de Murnaghan est utilisée pour ajuster l'énergie totale en 
fonction du volume et pour déterminer les paramètres structuraux. Les configurations électroniques 
utilisées dans ce travail pour caractériser les différents atomes sont les suivantes : Mg [Ne] 3 s?,O 
[He] 2 s° 2 p*, Cu [Ar] 3 d!°4s!, Ag [Kr] 4 d!°5s!, Au [Xe] 5 d'°6s!. L'intégration de la première 
zone de Brillouin est réalisée avec un maillage Monkhorst Pack de dimensions (12 x 12 x 12) [12], 
et un maillage de points k de dimensions (17 x 17 x 17) est utilisé pour échantillonner et pour 
obtenir les propriétés électroniques. Les paramètres d'entrée, à savoir Lmax, RKmax €t Gmax SOnt fixés 
respectivement à 7, 10 et 12 pour toutes les configurations étudiées. La séparation entre l'état de 
valence et de cœur est établie à de -6 Ry. Les cycles auto-cohérents sont interrompus lorsque 
l'énergie totale converge vers 10* Rydberg. Pour illustrer la structure cristalline, le logiciel 


Xcrysden est utilisé [13]. 


IV. 3. Résultats et discussions 
IV.3.1. Propriétés structurales 


L'oxyde M£O cristallise dans une structure type NaCI (B1) dans des conditions normales 
[14]. Pour optimiser les paramètres structuraux et étudier la stabilité des phases structurales, les 
courbes des énergies totales en fonction du volume sont ajustées en utilisant l'équation d'état 
empirique de Murnaghan [15]. Les caractéristiques structurales des systèmes TM;Mg: ;0 (x = 
3.125 %, 6.25 %, 12.5 % et 25 %) sont obtenues par substitution d'un atome de TM dans une super 
cellule MgO. Un atome de Mg est remplacé par un atome de TM dans les trois types de super 
cellules cubiques (respectivement SC, BCC et FCC) pour les concentrations (x = 3.125 %, 6.25 %, 
12.5 %) et deux atomes de M£ pour la concentration 25 %. Les structures géométriques sont 


illustrées par la figure IV. 1. 


Les figures (IV. 2-5) représentent, l'énergie calculée, dans les phases non magnétiques 
NM (état non spin polarisé) et ferromagnétiques FM (état spin polarisé) des composés dopés 
TMo25Mg0.750, TMo.125M 20.8750, TMo.062Mg0.9380, TMo.031Mgo.60, en fonction du volume. Tous 
les systèmes susmentionnés sont énergétiquement stables dans l’état de spin polarisé, avec une 


énergie plus faible à un volume donné. 
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Figure IV. 1. Cellules cristallines TM;Mg:1.:0 (A) 25% (FCC), (B) 12.5% (FCC), (C) 6.25% 
(BCC), et (D) 3.125 % (SC). 
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Figure IV. 2. Variation de l’énergie totale en fonction du volume pour TMo25Mg0750 (TM = Cu, 
Ag, et Au) dans les états NM et FM. 
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Figure IV. 3. Variation de l’énergie totale en fonction du volume pour TMo.125Mg0.8750 (TM = 
Cu, Ag, et Au) dans les états NM et FM. 
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Figure IV. 4. Variation de l’énergie totale en fonction du volume pour TMo.062Mg0.9380 (TM = 
Cu, Ag, et Au) dans les états NM et FM. 
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Figure IV. 5. Variation de l’énergie totale en fonction du volume pour TMo.031Mg0.9690 (TM = 
Cu, Ag, et Au) dans les états NM et FM. 
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Les valeurs de la constante de réseau (a), du module de compressibilité (B) et de sa première 
dérivée (B') pour les composés ternaire TM;M2£: ;O, calculées à l'aide de l'approximation GGA- 


WC, sont présentées dans le tableau IV. 1. Ces résultats sont comparés aux données théoriques de 


la littérature. 


Tableau IV.1. Valeurs calculées de la constante de réseau a (À), module de compressibilité B 


(GPa), et sa dérivé B’° des composés TM;Mg:1.x O (TM = Cu, Ag, et Au ; x = 3.125 %, 6.25 %, 12.5 


% et 25 %). 

CuxMg1:x O a (A) B (GPa) B’ 
Cu/MgO (CS) 3.125 % 4.2196  158.8178 4.2063 
Cu/Mg0O (CC) 6.25 % 4.2194  159.8691 4.2018 
Cu/Mg0O (FCC) 125% 4.2192  161.8911 4.3052 
Cu/Mg0O (FCC) 25 % 4.2073  169.4236 4.3283 

AgMg1x O a(À)  B(GPa) B’ 
Ag/MgO (CS) 3.125 % 4.242 156.024 4.124 
Ag/MgO (CC) 6.25 % 4.264 156.353 4.124 
Ag/MgO (FCC) 12.5 % 4.303 157.142 4413 
Ag/MgO (FCC) 25 % 4.354 159.359 4.191 
Résultas Théoriques 12.5 % 4.37 [16] 

AuxMg1+ O a(À)  B(GPa) B’ 
Au/MgO (CS) 3.125 % 4.2541  156.3004 4.0802 
Au/Mg0O (CC) 6.25 % 4.2841 155.7414 4.1654 
Au/Mg0O (FCC) 12.5 % 4.3395  154.9697 4.1293 
Au/Mg0O (FCC) 25 % 44100  154.5794 4.1889 
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L'analyse des données de ce tableau montre que la constante de réseau calculée, des 
systèmes dopés, augmente avec le numéro atomique du dopant pour toutes les concentrations, 
puisque plus le rayon ioniques des atomes dopants est grand Cu (0.54 À), Ag (1.26 À), et Au (1.37 


À), plus la constante de réseau est grande. 


La Figure IV. 6 montre la variation du paramètre de réseau en fonction de la concentration 
des dopants. Les valeurs des paramètres de réseau des composés ternaires Ag;Me1 ,;O et Aux;Meg:- 
xO sont supérieures à celles de la structure binaire en raison de l'effet de taille des rayons ioniques 
entre Mg (0.72), Ag (1.26 À) et Au (1.37 À). Cette augmentation du paramètre de réseau est 
d'autant plus marquée que la concentration en éléments dopants Ag ou Au s'accroît, reflétant 
l'influence directe de la substitution de Mg par des ions plus larges sur la structure cristalline. En 
revanche, hormis pour l'ion Cu, une tendance générale à la diminution du paramètre de maille est 
observée avec l'augmentation de la concentration des dopants comparativement à celle de la 
matrice hôte MgO, phénomène qui peut être expliqué par la différence de rayon ionique entre les 
1ons concernés. Cette variation, à la fois d'augmentation pour certains et de réduction pour d'autres, 
des paramètres structuraux est cohérente avec des observations similaires rapportées dans la 


littérature pour d'autres systèmes dopés tels que Ag-ZrO2, Ca-CeO2 Cu-MgO et Ag-CeO [17-201]. 


Par ailleurs, il convient de noter l'absence d’études théoriques permettant une comparaison 
avec les valeurs prédites des propriétés structurales pour les concentrations de x = 3.125%, 6.25%, 
et 25% pour les trois métaux de transitions. Uniquement, pour le composé TMo.125Mg0.8750, les 
résultats sont comparables à ceux publiés par Teli et al. [16] ayant utilisé la symétrie Fm3m. Ainsi, 


ces résultats constituent une référence théorique pour orienter les futures recherches. 


En conséquence, les propriétés électroniques et magnétiques des composés étudiés sont 


évaluées en utilisant les constantes de réseau calculées dans l’état ferromagnétique spin polarisé. 
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Figure IV. 6. Variation du paramètre de réseau a (À) en fonction de la concentration des dopants 
en utilisant l'approximation GGA-WC. 


IV.3.2. Propriétés électroniques 


Dans cette section, les propriétés électroniques des composés étudiés, telles que les 
structures de bandes, les densités d’états et les densités de charges, sont examinées par les deux 
approches GGA-WC et GGA-WC-mBJ. Les résultats obtenus sont comparés les uns aux autres et 
avec d’autres résultats disponibles. Désormais, toutes les propriétés mentionnées sont simulées en 


utilisant la constante de réseau d'équilibre précédemment calculée avec l'approximation GGA-WC. 


IV. 3. 2.1. Structures de bandes d’énergie (SB) 


Les structures de bandes de M£O pur dans la phase non magnétique (NM) sont calculées 
et présenté dans le chapitre (SIL. 1. 4. 2. 1). Le composé hôte MgO (x = 0%) est un semi- 
conductrice, avec une bande interdite directe aux points de symétrie I” égale à 4.49 eV en utilisant 
GGA-WC, tandis qu’en utilisant l'approximation WC-mB)J, la valeur du gap est égale à 7.15 eV, 
cette valeur est proche des données expérimentales. Les structures de bandes électroniques 
polarisées en spin des composés TM;Me:.,(TM = Cu, Ag, et Au ; x = 3.125 %, 6.25 %, 12.5 %et 
25 %) (Fig. VI 8-10) montrent un caractère métallique dans le spin minoritaire (spin down) car la 
bande d'énergie traverse le niveau de Fermi au point de I” de haute symétrie. Pour le spin majoritaire 
(spin up), seuls les composés TMo25Mg0750 (Cu, Ag, et Au) présentent un comportement 
métallique. En revanche, TMo:125Mg0.8750, TMo.062M 2809380, TMo031Mg0960O préservent le 


comportement semi-conducteur. Les structures de bande de TMo.125Mg0.8750, TMo.062Mg0.9380 


122 


indiquent une bande interdite directe, alors que TMo031Mg0960 présentent une bande interdite 
indirecte qui passe du point L'au point X pour le composé Cuo.031Mg096°0 et de T° au point M pour 
le composé Ago.031Mg0.690 en utilisant les deux approximations GGA-WC et WC-mBJ. Pour le 
composé Auo.031Mg0.9690, la bande interdite indirecte passe du point R au point I’. En conséquence, 
pour des concentrations de x = 3.125 %, 6.25 % et 12.5 %, les nouveaux composés ternaires sont 
ferromagnétiques semi-métalliques avec une polarisation de spin de 100 %. Le même phénomène, 
le caractère métallique, est constaté lorsqu'un essai est effectué avec une autre super cellule (1x1x1) 


a une concentration de 25 %. 


Les bandes de spin minoritaires du MgO dopé au TM (Cu, Ag, et Au) présentent deux 
caractéristiques fondamentales qui revêtent une grande importance dans le domaine des dispositifs 
spintronique. La première caractéristique est le gap ferromagnétique (Eg), qui correspond à la 
différence d'énergie entre le VBM et le CBM. Ce gap ferromagnétique définit la bande interdite 
fondamentale. La seconde caractéristique, le gap demi-métallique (Gxm), représente l'écart 
énergétique entre EF et le niveau d'énergie le plus proche entre le bas de la bande de conduction et 
le haut de la bande de valence [21]. Le tableau IV. 2 présentes les résultats calculés des gaps 


énergétiques Eg (eV) et des gaps demi-métalliques (Gxm) des composés dopés. 
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Figure IV. 7. Variations de gap E, et gap Gam des composés TM;M£:.,0 en fonction de la 
concentration avec les deux approximations GGA-WC et TB-mBJ. 


Les valeurs des gaps énergétiques calculées pour les composés TM;M£:.0 dans le canal 
de spin majoritaire sont de 2.59 eV, 3.43 eV et 3.62 eV pour les composés Cu;M£::0. Pour les 
composés Ag,M21:0, les valeurs sont de 1.93 eV, 2.41 eV et 2.55 eV et elles sont de 0.51 eV, 
1.19 eV et 1.38 eV, pour le composé Au;M£1,0, correspondant aux concentrations respectives de 
x = 12.5 %, 6.25 %, et 3.125 %. Ces écarts énergétiques sont évalués par rapport au niveau de 


Fermi en utilisant la GGA-WC. 


Les valeurs sont de 3.11, 5.93, et 6.09 eV pour Cu;M£g1,0. Pour Ag;Mo1.,0 les valeurs 
équivalentes sont de 4.32 eV, 4.76 eV, 4.89 eV. Pour Au;M21.,0, les écarts énergétiques sont de 
2.52 eV, 3.21 eV, 3.43 eV pour les concentrations respectives de x = 12.5 %, 6.25 %, et 3.125 % 
en utilisant l'approche WC-mBJ. Ces résultats mettent en lumière les propriétés énergétiques 
distinctes de chaque composé dopé, fournissant des informations cruciales pour la compréhension 


et la conception de dispositifs basés sur ces matériaux dans le domaine de la spintronique. 


Les données du tableau démontrent une tendance significative : l'énergie de la bande 
interdite diminue à mesure que la concentration augmente pour les trois métaux de transition, à 
savoir le cuivre (Cu), l'argent (Ag) et l'or (Au). Les résultats d’après le tableau et la figure IV. 7. 
montrent que l'énergie de la bande interdite diminue avec l'augmentation de la concentration pour 


les trois métaux de transition (Cu, Ag et Au), jusqu’à x — 0.125. 


Plus la teneur en TM est élevée, plus les énergies de la bande interdite sont faibles. Le 


rétrécissement observé de la bande interdite dans les matériaux M£O dopés au TM est dû à la 
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présence des niveaux d'énergie 3d, 4d et 5d supplémentaires dans la structure de la bande 
électronique et la perturbation du réseau cristallin du l'incorporation de ces ions dopant dans le 


réseau cristallin qui modifie les interactions entre les ions Mg et O. 


Une tendance monotone est observée, en augmentant le numéro atomique du dopant, 
l'énergie de gap diminue pour toutes les concentrations. Cette diminution peut être attribuée aux 
électrons des couches internes des dopants 4d et 5d qui ont un effet d'écrantage plus important que 
ceux des dopants 3d et un élargissement des orbitales des dopants 4d et 5d qui sont plus diffuses 


que celles des dopants 3d, ce qui favorise la délocalisation des électrons et la diminution de Eg. 
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Figure IV. 8. Structures de bandes polarisées en spin des composés Cu;M£1,0 calculées avec 


l’approximation GGA-WC et TB-mBJ. 


Les résultats pour les composés TMo.125Mg0.8750, TMo.062Mg0.9380, TMo.031Mg0.9690 pour 
tous les dopants révèlent une diminution du Gam avec l'augmentation de la concentration de dopage 
de 3.125 % à 6.25 %, jusqu’à 12.5 % tels que illustrés par les courbes de la figure IV. 7. Sur une 
perspective du numéro atomique, l'analyse des données montre que les valeurs calculées, du gap 
Guam des systèmes dopés, diminuent avec l’augmentation de numéro atomique du dopant pour 


toutes les concentrations. Ces valeurs étant assez élevées suggèrent la possibilité d’atteindre des 
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températures de Curie élevées. Les valeurs de Eg obtenues via WC-mBJ dans ce travail sont 
comparées à celles obtenues par H. Absike et coll. (l’approche TB-mBJ). Les résultats montrent 
une concordance, suggérant que la méthode mBJ améliore considérablement la précision de la 
bande interdite [22]. Ces observations renforcent la robustesse des résultats obtenus et mettent en 
évidence la pertinence de l'approche mBJ dans l'évaluation des propriétés électroniques des 


composés dopés [201]. 
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Figure IV. 9. Structures de bandes polarisées en spin des composés Ag;M£1,0 calculées avec 


l’approximation GGA-WC et TB-mBJ. 
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Figure IV. 10. Structures de bandes polarisées en spin, des composés Au;M£g:1.,0, calculées avec 


l’approximation GGA-WC et TB-mBJ. 


Tableau IV. 2. Gaps énergétiques Eg (eV) et gap demi-métallique Gam des composé TM;Mg1.x0 
calculés avec les deux approximations GGA-WC et l’approche TB-mBJ 


Energie (eV) 


Energie (eV) 


Energie (eV) 


Energie (eV) 


CuxMg1x O 

Cu 0.031 Mg 0.268 O (T-M) 6.09 2.20 TB-mBJ 
3.62 0.77 GGA-WC 

Cu 0.062 Mg 0.938 O (T'-T) 5.93 2.14 TB-mBJ 
3.43 0.66 GGA-WC 

Cu 0125 Mg 0.875 O (T-F) LL 2.00 TB-mBJ 
2.59 0.53 GGA-WC 

Cu 0.25 Mg 075 O Metallic Metallic TB-mBJ 
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Metallic Metallic GGA-WC 


Autres calculs 


MgO 4.49 [16] 


AgxMg1-x O 
Ag 0.031 M£g 0.969 O (T-M) 4.89 1.084 WC-mBJ 
(SC) 295 0.42 GGA-WC 
Ag 0.062 M£g 0.938 O (T°-T°) 4.76 1.03 WC-mBJ 
(BCC) 2.41 0.34 GGA-WC 
Ag 0.125 M£g 0.875 O (T°-T) 4.32 0.85 WC-mBJ 
(FCC) 1.93 0.16 GGA-WC 
Ag 0.25 Mg 0.75 O (FCC) Metallic Metallic WC-mBJ 
Metallic Metallic GGA-WC 
Autre calculs 
Ag 0.125 Mg 0.875 O (FCC) 2.20 [16] / GGA 


AuxMg1-x O 

Au 0.031 Mg 0.969 O (R:- T°) 3.43 1.02 WC-mBJ 
(SC) 1.38 0.36 GGA-WC 
Au 0.062 Mg 0.938 O (T°-F) 3.21 1.01 WC-mBJ 
(BCC) 1.19 0.35 GGA-WC 
Au 0.125 Mg 0.875 O (T°-T) 2:52 0.83 WC-mBJ 
(ŒCC) 0.51 0.14 GGA-WC 
Au 0.25 Mg 0.75 O (FCC) Metallic Metallic WC-mBJ 


Metallic Metallic GGA-WC 


IV.3. 2.2. Densités d’états électroniques (DOS) 


Pour mieux comprendre et expliquer les propriétés électroniques ainsi que l'origine du 
ferromagnétisme au sein des composés TM;M£g::0, les courbes de densités d'états totales (TDOS) 
et de densités d'états partielles (PDOS) dans la première zone de Brillouin (BZ) pour les deux 
orientations de spin up et down sont examinés. Les densités d'états totales et partielles constituent 


une description fondamentale de la nature du ferromagnétisme, en fournissant des informations à 
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la fois qualitatives et quantitatives sur les contributions atomiques à divers états électroniques 
présents dans les bandes d'énergie des composés TM;M2g1:0. Les courbes TDOS et PDOS sont 
présentées dans les Fig. 11 jusqu’à 13, où l'énergie de Fermi est prise égale à zéro. Les courbes 
sont tracées sur un intervalle d'énergie de — 10 à 10 eV. La densité d'états est caractérisée par deux 
régions essentielles : la bande de valence et la bande de conduction. En utilisant la densité d’état 


partielle, il est possible de définir le caractère de chaque région pour MgO et TM;M£::x0. 


Pour la densité totale et partielle des états de MgO non dopé (x = 0), la région de —-5 à 0 
eV, au niveau de la bande de valence, les orbitales O-2p dominent. À l'inverse, dans la région de 
la bande de conduction inférieure, qui s’étend d’environ 4.6 à 15 eV, une contribution notable 
provient des états O-2p ainsi que les états 3 s et 2p de l’atome de Mg (($SIIL. 1. 4. 2. 2). Comme 
mentionné précédemment, dans la partie de la structure de bande, les composés TMo.125Mg0.8750, 
TMo.062M 80.380, TMo.031Mg0.690, se comportent en tant que semi-conducteurs dans le canal de 
spin up et comme des conducteurs métalliques dans le canal de spin down. En revanche, le caractère 
métallique du composé TMo25M 0750 est justifié par la présence d'états électroniques au niveau 
de Fermi dans les deux orientations de spin. Ces observations concordent avec les résultats de 
densités d’états électroniques. En examinant les courbes PDOS (Fig. IV. 11-13) pour le cas x = 
25%, la bande de valence est dominée par l'état O-p dans la plage [-— 5, 3.5] et l'état d-TM dans la 
plage [— 3.5, 2.5]. Au niveau de Fermi, les orbitales d-eg-TM se dispersent dans les deux canaux 
de spin et s'élargissent considérablement dans la plage d'énergie compris entre — 2 et 1 eV, ce qui 
confirme le caractère métallique du composé TMo25Mg0.750. Dans le cas de x = 3.125 %, 6.25 % 
et 12.5 %, les courbes PDOS montrent qu'au maximum de la bande de valence et près du niveau 
de Fermi, 1l y a une asymétrie entre le spin-up et le spin-down. Les états TM-3d, 4d et 5d se 
chevauchent avec ceux de O-2p à x = 12.5%, au niveau de Fermi, en raison de la forte hybridation 
subsistant entre les états TM-d et O-p et de la division importante par échange de spin des états 
TM-d. Selon la théorie du champ cristallin, les cinq orbitales dégénérées des orbitales d-TM sont 
divisées en deux parties : les trois orbitales à basse énergie (t29) (dxz, dxy, dyz) et les deux orbitales 


à haute énergie générées (eg) (dx2-dy2, et d2), de manières symétriques. 


Par conséquent, cela donne lieu à la création d'états liants et anti-liants, jouant un rôle 
majeur dans l'induction du ferromagnétisme HM dans TM;M£g1:0 pour les concentrations x = 0.6, 
et 0.3. Les orbitales 4d (t24,) sont localisées au niveau de Fermi. Les états TM- 3d, 4d et 5d sont 


repoussés du niveau de Fermi, ce qui génère un écart énergétique qui donne le comportement semi- 
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conducteur aux composés résultants. Dans la partie inférieure de la BV, de — 7 à — 2 eV, les orbitales 
TM-d t>, sont dominantes avec une contribution significative des états O-p. Dans la bande de 
conduction, de 4.5 à 11 eV, la contribution principale provient des atomes O-p, avec une petite 
contribution des états s-Mg. La contribution de s-Mg et de s-O, en général est négligeable. À partir 
des PDOS, les valeurs de Eg diminuent lorsque la concentration de dopage augmente pour tous les 
dopants (Cu, Ag, Au). Alors qu'à x = 3.125 %, la valeur est de 3.62, 2.55, et 1.38 eV, pour Cu, Ag 
et Au respectivement avec GGA-WC. Des valeurs de 6.09, 4.89, et 3.43 eV pour Cu, Ag, et Au 
avec WC-mBJ., Les orbitales TM- 3d, 4d, 5d-eg sont repoussées vers EF et Eg diminue avec 
l’augmentation des concentrations et du numéro atomique des dopants. À forte concentration, x = 
25 %, les composés deviennent métalliques en utilisant les deux approximations. Ce résultat est 


similaire à celui observé dans ZnO dopé avec Ag [23]. 
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Figure IV. 11. TDOS et PDOS calculées des composés Cu;Mg10, par GGA-WC et WC-mBJ. 
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Figure IV. 12. TDOS et PDOS calculées des composés Ag,;Mg:1.0, par GGA-WC et WC-mBJ. 
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Figure IV. 13. TDOS et PDOS calculées des composés Au;M£1,0, par GGA-WC et WC-mBJ. 


IV.3.2.3 Densité de charge électronique 


Pour illustrer la nature des liaisons chimiques entre les atomes et décrire leur transfert de 
charge, la densité de charge de MgO non dopé et de MgO dopé à l’Ag pour des concentrations de 
= 3.125 %, 6.25 %, 12.5 % et 25 % est calculée en utilisant la méthode WC-mBJ, comme indiqué 
sur les figures. IV. 14 et 15. L'atome de Mg est moins électronégatif que les atomes d'Ag et d'O. 


Leurs valeurs d'électronégativité sont respectivement de 1.31, 1.93 et 3.44. 


La densité de charge de M£O projetée sur le plan (100) montre une nature ionique, car les 
contours électroniques de charge sont isolés. Pour les composés dopés, la Fig. IV. 15 indique 
clairement qu'ils présentent une nature de liaison mixte. La liaison M£g-O est ionique, tandis que 
pour Mg-Ag et Ag-O, il y a un chevauchement entre le cation (Ag*) et l'anion (07°), ainsi qu’entre 
(Ag*) et (Mg*?), ceci peut être considérer comme une liaison covalente, Cela implique l'existence 
d'une charge partagée (électron) entre les atomes liés pour les deux canaux de spin. Les mêmes 
liaisons covalentes se forment lorsque MgO, SrO et CaO sont dopés avec Ag et Cu [24], la règle 
de Pauling permet d'identifier le type de liaison entre les atomes. C'était une méthode très efficace, 


qui estime le facteur d'iconicité Fi par l'expression suivante [25]: 


(Xa — XB)? 


F, = 1— exp(— TD (IV. 1.) 


Où XA et X8B représentent respectivement les électronégativités des atomes A et B. Une valeur "F;" 
proche de zéro indique la prédominance du caractère covalent, tandis qu'une valeur proche de 1 
représente une liaison ionique. Les résultats du facteur d'iconicité sont : Fi : Mg-O = 0.98, F; : Mg- 


Ag = 0.10, Fi: Ag — O = 0.44, qui confirment les résultats des densités de charges. 
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Figure IV. 14. Densité de charge de MgO non dopé dans le plan (100) utilisant la fonctionnalité 
WC-mBJ. 
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Figure IV. 15. Courbe des densités de charge des composés Ag,Mg1.,0 dans le plan (100) en 
utilisant WC-mBJ. 
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IV. 3.3 Propriétés magnétiques 
IV. 3.3. 1. Moment magnétique 


Le tableau. IV. 3. résume les moments magnétiques totaux et locaux calculés par atomes 
TM = Cu, Ag, Au, Meg, O ainsi que pour les sites interstitiels en utilisant les approximations GGA- 
WC et l’approche WC-mBJ. Le composé hôte MgO est non magnétique. Après dopage avec des 
atomes de métaux de transition (Cu, Ag et Au) non magnétiques à différentes concentrations (x = 
3.125 %, 6.25 %, 12.5 % et 25 %), les composés dopés sont magnétiques. Ces résultats concordent 
avec ceux du M£O dopé au Cu [26], du ZnO dopé avec Ag [23] SnO: dopé à Ag [27] et TiO2 dopé 
(Cu, Ag et Au) [28]. 


En raison du transfert d'électrons de type p-d, de petits moments magnétiques locaux 
permanents sont produits sur les sites non magnétiques des atomes de Mg et O, qui présentent un 
signe positif indiquant des arrangements de spins parallèles. Les atomes de Cu, Ag et Au contribue 
considérablement au moment magnétique total des composés, et la valeur de moment magnétique 
diminue à mesure que la concentration augmente. Le site interstitiel contribue également au 


moment magnétique total de la super cellule avec des valeurs positives. 


Les composés ternaires TM;M£g1:0 sont obtenus en incorporant de Cu, Ag et Au avec la 
configuration électronique 4s! 3d!° (t26 — eg), 5s! 4d!0 (té — eg), et 6s! 5d!° (t$ — eg‘), 
respectivement. Le moment magnétique provient de la division des états d des atomes de métaux 
de transition Cu, Ag, et Au (t2° — eg*). Les atomes de transition (TM) contribuent avec un seul 
électron aux états de liaison des porteurs des bandes de valence et deviennent Cu *! [4s° 3d!° (t2,f 
— eg*)], Ag*! [5s° 4d!° (125 — es*)], et Au *! [65° 5d!° (1,2 — eg*)]. Alors que l'atome Mg possède 
deux électrons de valence (Mg*?), ce qui respecte la loi de Hund et Pauli [29]. Selon les mécanismes 
du M£gO dopé au TM, il peut être décrit par la création d'états d'impuretés attribués aux atomes de 
TM et par l'hybridation entre les orbitales TM-3d, 4d et Sd, et les états O-2p proches pour 12.5% 


et TM-d pour les autres concentrations. 


Les résultats calculés montrent que TMo.:125Mg0.8750, TMo.062M£0.9380, TMo.031Mg£0.9690 
offrent à l'expérimentateur une nouvelle voie pour des applications potentielles dans les dispositifs 


de spintronique. 
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Tableau VI. 3. Moment magnétique total et local en 18 pour les composés TM;M£:,0 en utilisant 


l’approximation GGA-WC et l’approche TB-mBJ. 


Cux Mgix O 
Cu 0.031 Mg 0.268 O GGA-WC 0.9 0.62479 0.00108 0.05734 0.01451 
TB-mBJ 1 0.79578 0.00167 0.00085 0.05081 
Cu 0.062 Mg 0.938 O GGA-WC 0.9 0.62377 0.00182 0.05602 0.00262 
TB-mBJ 1 0.790468 0.00052 0.04351 0.05324 
Cu 0.125 Mg 0.875 O0 GGA-WC 0.9 0.61310 0.00642 0.06286 0.00144 
TB-mBJ 1 0.770945 0.00090 0.04802 0.06371 
Cu 025 Mg 0750 GGA-WC 0.0 0.00410 0.0002 0.00046 0.00038 
TB-mBJ 0.1 0.10137 0.0000 0.00593 0.00130 
«Composé Approximations Total(w)  AgQu)  Mgu) Ou)  Interstitiel(u) 
Agx Mgix O 
Ag 0.031 Mg 0.268 O GGA-WC 1 0.47884 0.00589 0.07544 0.05722 
TB-mBJ 1 0.55622 0.00347 0.08396 0.06449 
Ag 0.062 Mg 0.938 O0 GGA-WC 0.9 0.47369 0.00415 0.07493 0.04735 
TB-mBJ 1 0.55209 0.00303 0.08779 0.06779 
Ag 0.125 Mg 0.875 O0 GGA-WC 0.9 0.46208 0.00754 0.08727 0.03439 
TB-mBJ 1 0.53697 0.00349 0.09764 0.10750 
Ag 0.25 Mg 0750 GGA-WC 0.0 0.00571 0.00028 0.00266 0.00912 
TB-mBJ 0.0 0.00187 0.00001 0.00838 0.00064 


AUx Mgix O 

AU 0.031 Mg 0.968 O GGA-WC 0.9 0.51987 0.005412 0.05734 0.07370 
TB-mBJ Il 0.59052 0.00405 0.0664 0.039839 

Au 0.062 Mg 0.938 O GGA-WC 1 0.51813 0.00746 0.06877 0.07254 
TB-mBJ 1 0.59005 0.00368 0.07607 0.05347 

Au 0.125 Mg 0.875 0 GGA-WC 0.9 0.51322 0.00891 0.07725 0.09072 
TB-mBJ 1 0.58076 0.00442 0.08594 0.038316 

Au 0.25 Mg 075 O GGA-WC 0.004 0.004 0.00005 0.00015 0.00489 
TB-mBJ 0.000 0.0000 0.00000 0.00001 0.00000 


IV.3.3.2. Constantes de couplage d’échange 


Les constantes d'échanges s-d (Noa) de la bande de conduction et p-d (No) de la bande de 


valence sont calculées pour déterminer la nature de l'attraction dans TM;Mg:,0 selon les 


expressions de la théorie des champs suivante : 


(IV. 2) 


(IV. 3) 
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Où AE. et AE, représentent respectivement l'énergie de séparation des bords de la bande de 


conduction et de valence. Leurs équations sont : 
AEc = EcB Min — ECB Min (V.4) 
AE, = Ep max — ECB Max AV.5) 
x = 0.250, 0.125, 0.062, 0.031, sont les concentrations d'impuretés (TM) 
< s > est le moitié de moment magnétique total calculé par ion de métaux de transition. 


Le tableau IV.4 résume tous les résultats calculés à savoir AËk, AË,, Noa et Nob. Les valeurs de 
Noa et Nof sont positives pour tous les composés avec différentes concentrations, indiquant un 
couplage ferromagnétique entre les états d (TM) et la bande de conduction et de valence. 
Cependant, pour les deux composé Cuo12:Mg08750 et Auo.031Mgo680, il semble que le couplage 
d'échange sd entre la bande de conduction des niveaux de M£O et de d-Cu soit antiferromagnétique, 
expliquant le signe négatif de la constante Noa. De plus, la constante No est négative pour 
AUo.031Mg0.9680, suggérant un échange p-d entre la bande de valence et les niveaux de d-Au de type 


antiferromagnétique. 


Noa et NB augmentent lorsque la concentration des atomes d’impurités diminue pour le 
cas de Ag,Me1.,0 et Au;M£1,0. Par contre, dans le cas de Cu;M£:1.,0 les valeurs de Noa augmente 
puis diminue avec la diminution de la concentration. La constante d'échange conserve la valeur 
positive, indiquant que le couplage est ferromagnétique et confère le caractère FM. Ces résultats 


de prédiction sont en accord avec les résultats théoriques rapportés pour d’autres DMS [30-31]. 


Dans la plupart des travaux existants sur les semi-conducteurs HME, la valeur de Noa est 
généralement estimée comme positive, alors que la valeur de Nof est toujours négative [32-33]. 
Cette tendance est complètement inversée dans la situation actuelle. Néanmoins, la diminution de 
la symétrie et l’éventuel effet de confinement quantique pourraient être responsables du signe 
inversé de Noa et Nof. De plus, Mac et coll. sont les premiers à prédire le signe positif de Nof dans 
le semi-conducteur magnétique dilué Zn1,Cr,Se, attribuant le changement de signe notable à la 
position de l'état d dans la bande de valence [34]. Les résultats confirment l'effet magnétique 


retrouvé dans ce matériau. 
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Tableau IV. 4. Séparation des bords de bande de conduction AE (eV), séparation des bords de 
bande de valence AE, (eV) et constantes d'échange Noa et NOB en utilisant l'approche WC-mBJ. 


Cu 0.125 Mg 0.875 0 0.125 -0.024 2.1257 -0.384 34.011 
(ŒÆCC) 

Cu 0.062 Mg 0.938 O0 0.063 0.0486 2.1588 1.542 68.53 
(BCC) 

Cu 0.031 Mg 0.969 O 0.031 0.0045 2.1841 0.29 140.64 
(SC) 

Ag 0.125 Mg 0.875 O0 0.125 0.04114 0.44829 0.65824 7.17 
(FCO) 

Ag 0.062 Mg 0.938 O 0.063 0.02331 0.41884 0.7519 13.5109 
(BCC) 

Ag 0.031 Mg 0.969 O 0.031 0.038907 0.4303 2.514 27.7612 
(SC) 

Au 0.125 Mg 0.875 O0 0.125 0.00213 1.0834 0.0196 17.33 
(ŒÆCC) 

Au 0.062 Mg 0.938 O 0.063 0.0265 1.0699 0.0841 33.96 
(BCC) 

AU 0.031 Mg 0.969 O 0.031 0.0209 -0.9685 1.34 -6.248 
(SC) 


IV. 3. 4. Propriétés optiques 
Pour étudier les propriétés optiques des oxydes ternaires, la fonction diélectrique 
complexe en fonction de l'énergie des photons, £ (©) = £1 (@) + 1 2 (w), est calculée. &1 (@) et 2 


(wo) sont les parties réelles et imaginaires de la fonction diélectrique [35]. 


Les résultats des composantes réelles et imaginaires calculées de la fonction diélectrique 
sont présentés dans la Figure IV. 16 Les constantes diélectriques statiques obtenus à l'aide de la 
fonctionnelle TB-mBJ sont 1.70, 1.72, et 1.82 pour Cuo1:25Mgo8750, Ago125M 208750, et 
Auo.125Mg0.8750, respectivement, ce qui révèle une variation positive de £1 (@), suggérant ainsi un 
comportement diélectrique. La partie imaginaire £2 (©) émerge à un seuil d'environ 6.81 eV, 4.93 
eV, et 2.92 eV pour Cuo.:125Mg0.8750, Ago.125Mg0.8750, et Auo.125Mg0.8750, respectivement. Ce seuil 
est associé à la transition directe au point [résultant de la migration des électrons de la bande de 
valence à la bande de conduction. Les valeurs maximales de £2 (@) se situent autour de 13.86 eV, 


13.70 eV, et 14.29 eV, respectivement, pour Cuo.125M 208750, Ago.125M 20.8750, et Auo.125Mg0.8750. 
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L'absorption caractérise la profondeur à laquelle une longueur d'onde spécifique de la 
lumière peut pénétrer dans un matériau avant qu'elle ne soit absorbée. Cette propriété peut être 
calculée à partir de la partie imaginaire de la fonction diélectrique [36]. La Figure IV. 17 révèle 
une forte absorption pour tous les oxydes ternaires dans la région de 4 eV à 16 eV. Cette région de 
forte absorption présente un déplacement vers des énergies plus élevées à mesure que le numéro 
atomique passe de 3d à 4d puis à 5d. En raison des gaps énergétiques élevés, tous les spectres 
d’absorption des composés dopés sont essentiellement confinés dans la plage du spectre 


d’ultraviolet. 


L'indice de réfraction n(o) est un paramètre optique essentiel, mesurant la déviation d'un 
rayon lumineux lors de son passage d'un milieu à un autre. Les indices de réfraction calculés des 
composés M£O dopés sont présentés dans la Figure IV. 18 sur une plage d'énergie de 0 à 20 eV. 
Les indices de réfraction statiques sont 1.30, 1.31, et 1.34 pour les oxydes ternaires dopés 
Cuo:125Mg0.8750, Ago.:125Mg0.8750, et Auo125Mg0.8750 respectivement. Les valeurs calculées de 
l'indice de réfraction augmentent proportionnellement à l'énergie et atteignent des valeurs 
maximales autour de 13.61, 13.59, et 10.21 eV respectivement pour les composés Cuo.125Mg0.8750, 
A g0.125Mg0.8750, et Auo.125Mg0.8750, L'évolution de l'indice de réfraction n (w) suit la tendance de 
la partie réelle £1 (w). Ces valeurs se manifestent principalement dans la région ultraviolette 
suggérant des applications potentielles dans le domaine des LED et de détecteurs 


électromagnétiques UV. 


=— Cuo.125M9o.875O 

=— A9o.125M9o.875O 
J D 2 ———… AL, M le] 

2.0 = 0.125 M9o.875 


2.5 + (a) 


ë (@) 


ë, (@) 


o 2 2 6 8 10 12 14 16 18 20 22 24 
Energie (eV) 
Figure IV. 16. Partie réelle (a) et partie imaginaire (b) de la fonction diélectrique des oxydes 
ternaires Cuo.125Mgo.8750, Ago.125Mg0.8750, and Auo 125Mg0.8750 en utilisant l'approximation WC- 


mBJ. 
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Figure IV. 17. Coefficient d'absorption des oxydes ternaires Cuo125Mgo.8750, Ago.125Mg0.8750, et 
Auo.125Mg0.8750 en utilisant l'approximation WC-mBJ. 
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Figure IV. 18. Indice de réfraction pour oxydes ternaires Cuo.125Mg0.8750, Ago.125Mg0.8750, et 


Auo.125M£0.8750 en utilisant l'approximation WC-mBJ. 


IV. 4. Conclusion 


Cette étude a utilisé la méthode de calcul DET pour examiner le dopage de l'oxyde de 
magnésium M£gO par des métaux de transition 3d, 4d et 5d (Cu, Ag et Au) pour de faibles 
concentrations [x = 3.125 % (SC), 6.25 % (BCC), 12.5 % (FCC) et 25 % (FCC)]. La structure, 
TDOS, PDOS, la densité de charge, les propriétés optiques et les moments magnétiques sont 


rigoureusement analysés en utilisant les approximations GGA-WC et WC-mBJ. 
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Les paramètres du réseau augmentent en fonction de numéro atomique ainsi qu’avec la 
concentration de dopant pour Ag et Au. Par contre il démunie avec la concentration de dopant Cu. 
L'analyse des structures de bandes, des TDOS et PDOS Pour les composés ternaires, révèle que 
TMo:125Mg0.8750, TMo.o62M 20.380, TMo.031Mg09690 ont un caractère ferromagnétique semi- 
métallique, avec une transition d’un gap direct à un gap indirect et une polarisation de spin à 100 
%. Pour x = 0.25, le composé présente un comportement métallique. L'hybridation p-d entre les 
états 3d -4d et 5d des TM et les états 2p-O pour la concentration de 12.5 % (FCC) est responsable 
de la demi-métallicité et du ferromagnétisme induits dans les super cellules. Pour les autres 
concentrations, 3.125 % (SC) et 6.25 % (BCC), les états 3d, 4d et Sd génèrent le comportement 


semi-métallique. 


De plus, lorsque les concentrations de dopage des TM augmentent, une diminution du gap 
semi-métallique, passant de TMo.03:Mgo.9600 à TMoo62Mg09380 puis à TMo:125Mg08750 et un 
moment total par atome TM égal à 1 LB est observée, en utilisant l'approximation WC-mBJ. 
L'analyse de la densité de charge a identifié deux types de liaisons entre les atomes : une liaison 
covalente pour Ag-O, Mg-A£ et une liaison ionique associée M£g-O. Grâce au processus de dopage, 
les matériaux dopés présentent des changements notables dans les propriétés optiques, ouvrant 


ainsi la voie à un large éventail d’applications technologiques. 


Comparé au DMS conventionnel, MgO se distingue par l'absence d’ions magnétiques 
intrinsèques. Le dopage avec des ions de TM (Cu, Ag et Au) non magnétiques lui confère des 
propriétés magnétiques prometteuses, le rendant attractif pour la spintronique. Cette observation 
offre fournie une base théorique solide pour la synthèse et le développement de ces matériaux en 


vue de l’application future. 
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Chapitre V 


Propriétés structurales, magnétiques et optoélectroniques de MgO 
de type rock salt co-dopé avec Eu-TM (TM = Cu, Ag, et Au) en vue 
d’applications en spintronique et détecteurs UV 


Au cours de ce chapitre, les résultats de calculs des propriétés structurales, électroniques, optiques et magnétiques des composés 
co-dopés Euo.125TMo:125Mgo.7500, à une concentration x = 12.5 %, sont présentés et interprétés. Les calculs sont effectués en se 


basant sur l'approximation du gradient généralisé (GGA) pour décrire le potentiel d'échange et de corrélation. 


«Une théorie qui a donné de tels résultats doit contenir une grande part de vérité » 


(HA. Lorentz) 
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V.1. Introduction 


Le domaine du co-dopage de MgO demeure largement inexploré d'où l'importance de 
l’étude du co-dopage de MgO spécifiquement par les métaux de transition (3d (Cu), 4d (Ag) et 5d 
(Au)) combinés avec un élément des terres rares 4f (Eu). Ceci est effectué en remplaçant deux 
atomes de Mg par un atome de Eu et un atome parmi les métaux de transition (3d/4d et 5d). 
L'objectif de cette étude, basée sur des calculs de premiers principes, est de prédire les propriétés 


optiques, électroniques et magnétiques des composés co-dopés. 


V.2. Méthodes de calcul 


Cette étude vise à examiner en profondeur les propriétés essentielles de MgO co-dopé 
avec l'europium, ainsi que les trois métaux de transition (Cu, Ag et Au) à une concentration de x = 
12.5 % pour chaque élément. La méthode des ondes planes augmentées linéarisée à potentiel total 
plus orbitales locale (FP/LAPW+lo), implémentée dans le code Wien2Kk [1], est utilisée pour 
résoudre l'équation de Kohn-Sham dans le contexte de la théorie de la fonctionnelle de la densité 
polarisée en spin (SP-DFT) [2]. Les approximations du gradient généralisée (GGA-WC) et la 
correction de Becke-Johnson modifiée par Tran- Blaha (WC-TB-mBJ) sont sélectionnées pour 


traiter le potentiel d'échange-corrélation [3]. 


Une super cellule (2 x 2 x 2) est utilisée, où deux atomes de M£ sont remplacés par un 
atome de TM et un Eu, situés respectivement aux coins (0, 0, 0) et aux bords (0.5, 0, 0). Le 
paramètre de coupure Rur x Kmax est fixé à 9, avec une expansion angulaire maximale des 
fonctions d'onde de valence étant fixée à lnax = 12 et la densité de charge expansée de Fourier est 
estimée à Gmax = 14 (a.u) !. Les électrons Mg [3 s°], O [2 s? 2 p*], Cu [4 s! 3 d'], Ag [5 s! 4 dl", 
Au [6 s! 5 d!° 4 fM] et Eu [6 s? 4 f”] sont traités comme des électrons de valence. Un maillage de 
points k (10 x 10 x 10) est utilisé pour échantillonner la Zone de Brillouin (B2Z), composé d'un total 
de 1000 points [4]. De plus, les rayons des sphères muffin-Tin (Rur) sont égaux à 1.98, 2.00, 2.02, 
1.84, 1.86 et 2.02 (a, u) pour Cu, Ag, Au, Me, O et Eu, respectivement. 
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V.3. Résultats et discussion 
V.3.1. Propriétés structurales 


Pour calculer les paramètres structuraux optimisés de la matrice hôte MgO et des 
composés co-dopés Euo.125TMo.125Mg07500 (TM = Cu, Ag et Au), comme précédemment 
mentionné ($ IL. 1. 4. 1), l'approximation du gradient généralisé de Wu et Cohen (WC-GGA) est 
utilisée en raison à sa précision supérieure par rapport à la LDA et la GGA-PBE. L'énergie totale, 
la constante de réseau d'équilibre (a), le module de compressibilité (B), sa dérivée de pression (B') 
et les écarts énergétiques (AE = Exm-Erm) sont ajustées à l'équation d'états de Birch-Murnaghan 
donnée par l'expression suivante [5] : 


= BV _ prf1 + o8 — 
Eror(V) = Eo (V) + BBD) [86 ( nu) ne 1| (V.1) 


La figure V.1 montre les structures cristallines des composés Euo.125TMo.125Mg0 7500 (TM 
= Cu, Ag et Au) dessinées par le logiciel Vesta [6-7] ceux-ci ont des structures cristallines de type 


NaCI et appartiennent au groupe spatial Fm3m. 


Figure V. 1. Structure cristalline simulée des composés Euo.125TMo.125Mg07500. (TM = Cu, Ag 
et Au). 


La figure V. 2 représente l'optimisation énergétique totale, en fonction du volume, des 
composés co-dopés Euo.125TMo.125M 807500 dans les états ferromagnétique (FM) et non magnétique 
(NM). Il est clair que pour tous les composés, l’état FM est la phase la plus stable présentant une 


énergie inférieure à celle de l’état NM. 
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Figure V. 2. Variation de l’énergie totale en fonction du volume, des composés 
Euo.125TMo.125M£0.7500 (TM = Cu, Ag et Au), dans les deux états, non magnétique et 


ferromagnétique. 


Le Tableau V.1 présente les paramètres de réseau optimisés et le module de compressibilité 
issus de cette étude. Étant donné l'absence de donnée théorique ou expérimentale préalables dans 
la littérature pour les composés co-dopés Euo.1251Mo.125Mg07500, la présente étude est considérée 


comme une analyse prédictive détaillée. 


Tableau V.1 : Constante de réseau a (À), module de compressibilité B (GPa), sa dérivée B', écart 


énergétiques (AE) et l'énergie de formation de Euo.1251Mo.125Mg07500 (TM : Cu, Ag et Au). 


Eu 0.125 Cu 0.125 Mg 07500 FM  4.3911 FM  148.6474 FM 4.2690 0.4494  -6.4800 
NM 4.3903 PM  144.7959 NM 1.4238 

Eu 0125 Ag 0125 Mg 0.750 O0 FM  4.4620 FM  145.6023 FM 42617 04552  -6.1100 
NM 44526 PM  146.7815 NM  4.2924 

Eu 0125 Au 0.125 Mg 07500 FM  4.4906 FM  144.7312 FM 42722 0.4578  -6.0300 
NM 44824 PM  147.8287 NM  4.1528 
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Le composé Euo.125Auo.125Mg07500 à une constante de réseau plus grande que les autres 
composés, De plus, le module de compressibilité de Euo.125Cuo.123Mg0.7500 est supérieur à celui des 
deux autres composés. Le composé Euo.125Auo.125Mg0.7500 est plus résistant à la compression et a 
tendance à être rigide et moins déformable sous pression par rapport aux autres composés. 
L'énergie de formation AE diminue uniformément avec la diminution du numéro atomique et de la 
constante de réseau. En progression de Cu à Ag puis à Au la constante de réseau augmente. Celle- 
c1 est liée à l'augmentation des rayons ioniques respectifs de Cu (0.54 À), Ag (1.26 À) et Au (1.37 


À) ainsi qu'au numéro atomique des TM. 


Pour évaluer la possibilité de synthétiser expérimentalement ces composés, une analyse 
de leur stabilité thermodynamique est réalisée en calculant l'énergie de formation Er donnée par 


l’équation [8] : 
A Hs = Erota (codoped MgO) — A Erm — 1 Epu — M Emg — n Eo (V.2) 


Où Eroui représente les énergies des systèmes co-dopés, et À Erm, lEu, M Emgin Eo représentent 
respectivement les énergies global calculées par DFT de Cu, Ag, Au, Eu, Mg et O dans la phase 
cubique. Les énergies de formation calculées sont également répertoriées dans le tableau V.1. Le 
signe négatif de ces énergies indique que ces systèmes sont relativement stables. En raison de sa 
plus faible énergie de formation, Euo.125Auo.125Mg07500 est la configuration la plus stable parmi les 
différentes structures. Les paramètres optimisés répertoriés dans le tableau V.1 sont utilisés pour 
calculer les autres propriétés électroniques, optiques et magnétiques en utilisant l'approximation 


WC-GGA et la fonctionnelle WC-mBJ. 


V. 3.2. Propriétés électroniques 


La plupart des propriétés physiques des solides sont régies par leurs structures de 
bandes électroniques et leurs densités d'états. Par conséquent, la connaissance de ces 
caractéristiques est essentielle pour leur utilisation efficace dans les dispositifs optoélectroniques, 
magnéto-optiques et électromagnétiques. Cette section étudie les propriétés électroniques du MgO 
co-dopé, dans la phase ferromagnétique, avec l’europium (Eu) et les métaux de transition (TM — 


Cu, Ag et Au). 


La structure de bandes du composé hôte MgO dans la phase non magnétique (NM) est 


simulée et présentée dans la figure ($SIIL. 1. 4. 2. 1) La GGA-WC sous-estime l'écart énergétique, 
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ce qui entraîne un gap direct égal à 4.61 eV. Les résultats pour la fonctionnelle TB-mBJ se situent 
dans la région des résultats expérimentaux avec une valeur de 7, 15 eV. Les figures V. 3. et V. 4. 
représentent les structures de bandes électroniques polarisées en spin des composés 
Euo.125TMo.125M 207500 (TM = Ag, Cu et Au) pour les spins majoritaires et minoritaires en utilisant 
les approximations GGA-WC et GGA mBJ, le long de la direction de plus haute symétrie dans la 


zone de Brillouin. 


L'approximation GGA-WC révèle un comportement demi-métallique pour les composés 
M£g0 co-dopés avec Eu-Ag et Eu-Cu, mais elle présente un caractère métallique pour le M£O co- 
dopé avec Eu-Au. Ces résultats sont conformes aux approximations limitées de la GGA pour 


calculer le gap énergétique. 


Eu Cup aMlgo0 
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Figure V. 3. Structures de bandes polarisées en spin des composés Euo.125TMo.125M 807500 (TM : 
Cu, Ag et Au) pour les états de spin up majoritaires (T) en couleur noire et les états de spin down 


minoritaires (|) en couleur bleue en utilisant l'approximation GGA-WC. 


Les valeurs de la bande interdite (gap énergétique) calculées pour Eg et HM sont 
répertoriées dans le tableau V. 2. Le calcul de la bande interdite est basé sur la différence entre le 
minimum de bande de conduction (CBM) et le maximum de bande de valence (VBM), où l'écart 
HM est défini comme la différence entre les valeurs de l'énergie supérieure des bandes de valence 
de spin minoritaires (ou majoritaires) et l'énergie inférieure des bandes de conduction de spin 
minoritaires (ou majoritaires) par rapport au niveau de Fermi. De plus, le CBM et le VBM se 
trouvent à des points symétriques différents, ce qui entraîne une bande interdite indirecte au point 
L-T, le VBM est déplacé vers le bas et s'éloigne du niveau de Fermi. À l’inverse, le CBM se déplace 
vers le bas et se rapproche de Er. En raison du co-dopage, la largeur de la bande interdite se rétrécit 
et de nouveaux niveaux d'énergie émergent autour du niveau de Fermi et dans la région de la bande 
de valence. Les électrons peuvent passer plus facilement de la bande de valence à la bande de 
conduction, cela se produit au sein des états majoritaires de Eu-4f, qui sont partiellement occupés, 
alors que tous les états minoritaires sont vides, ce qui indique la nature demi-métallique du système. 
Les valeurs de la bande interdite calculées à l'aide de l'approximation GGA-WC sont 


respectivement de 0.87, 0.60 et 0 eV pour TM = Cu, Ag et Au. 
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Tableau V.2 : Valeurs calculées, des gaps ferromagnétiques (Eg) et des gaps demi-métalliques 
(HM), des canaux à spin minoritaire, à l'aide des approximations GGA-WC et WC-mBJ pour les 


composés Euo.125TMo.125Mg0.7500 (TM : Cu, Ag et Au). 


Eu 025 Cu 0.25 Mg 0.750 O 0.87 0.19 GGA-WC 
2.53 1.18 WC- mBJ 
Eu 0125 Ag 0.125 Mg 0.750 O 0.60 0.06 GGA-WC 
2.49 0.80 WC- mBJ 
Eu 6.125 AU 0.125 Mg 0.750 O 0.00 / GGA-WC 
1.39 0.67 WC- mBJ 


Le potentiel d'échange Becke-Johnson modifié (mBJ) est employé pour améliorer les 
résultats des propriétés électroniques, à savoir la structure de bande et des densités d'états pour les 
composés Euo.125TMo.125M 207500. Cette méthode offre des améliorations notables de la largeur de 
bande interdite (Eg) des semi-conducteurs, des isolants, les positions de l'état d, et les propriétés 
optiques par rapport aux résultats GGA. Ces avantages découlent de l'utilisation du potentiel 
d'échange BJ qui démontre son efficacité pour divers types de matériaux. Ces observations peuvent 
être clairement déduites en examinant les résultats du tableau V. 2, où les valeurs de Eg et AM 
affichent un comportement linéaire en relation avec l'augmentation du numéro atomique des 
métaux de transition, passant de 3d à 4d à Sd. L'explication réside dans le processus de co-dopage 
de Eu : TM qui induit des états d'impureté dans la bande interdite de M£O, ce qui provoque un 
rétrécissement substantiel et abrupt de la bande interdite. À mesure que le numéro atomique 
augmente, la teneur en impuretés fait décaler simultanément le VBM et le CBM, avec une réduction 
nette de la bande interdite de 0.00 eV en utilisant GGA et de 1.39 eV en utilisant la méthode TB- 
mBJ-GGA par rapport à l'oxyde de M£gO pur. Ce décalage vers le bas est due à la diminution de 
l'interaction du cation M£ et de l'anion O avec le dopage du TM en fonction de la répulsion orbitale, 


ce qui réduit considérablement la bande interdite. 
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Figure V. 4. Structures de bandes polarisées en spin des composés Euo.125TMo.125Mg0.7500 (TM : 
Cu, Ag et Au) pour les états de spin up majoritaires (T) en couleur noire et les états de spin down 


minoritaires (|) en couleur bleue en utilisant l'approximation WC-mBJ. 


Les densités d'états fournissent une explication qualitative des états énergétiques de la CB 
et de la VB, ainsi que la distribution de chaque orbitale atomique. C’est pourquoi les densités d'états 
sont calculées, en utilisant la TB-mBJ, dans le domaine d'énergie de -10 à 10 eV pour les trois 


composés et sont présentés dans la figure V. 5. 


Le spectre TDOS peut être divisé en différentes parties. La partie inférieure de la bande 


de valence de -7 eV à -4 eV, est principalement constituée des états p de l'oxygène et un petit 


155 


nombre d'états d-eg des métaux de transition (Cu, Ag et Au). D'autre part, les états d-t24 des métaux 
de transition dominent la région comprise entre -3 eV et -2 eV avec une certaine contribution des 


états p-O. 


Au niveau de l'énergie de Fermi (Er), il existe une hybridation entre les états f-Eu, les 
états p-O et les états d-eg des métaux de transition (Cu, Ag et Au). Cette hybridation contribue à la 
formation du gap dans le canal spin-down et est à l'origine du comportement du HM. En outre, la 
bande de conduction montre que les états f et d-t2, de l’Eu dominent de 3 à 6 eV, tandis que p-O 
n'apporte qu'une petite contribution. Les valeurs des bandes interdites calculées à partir des 
structures de bandes sont similaires à celles calculées à partir des densités d'états. Cependant, les 
données théoriques et mesures expérimentales n'ont pas encore été rapportées dans la littérature, 


les résultats de ces travaux sont donc prédictifs pour les trois composés. 
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Figure V. 5. Densités d'états totales et partielles (spin polarisées) pour les composés 
Euo.125TMo.125Mg0.7500 (TM : Cu, Ag et Au) pour les états majoritaires de spin up (T) et 


minoritaires (|), en utilisant la fonctionnalité WC-mBJ. 


V. 3.3. Propriétés magnétique 


Le tableau V. 3 présente les moments magnétiques totaux, locaux et interstitiels pour Eu, 
Cu, Ag, Au, Mg et O, obtenus à partir de calculs polarisés en spin en utilisant les approximations 
WC-GGA et TB-mBJ. Il convient de noter que M£O non dopé présente un comportement non 
magnétique, car il ne possède aucun moment magnétique. Le co-dopage avec des atomes de métaux 
de transition des groupes 3d, 4d et 5d, et un élément de terres rares 4f, a abouti à la création de 


moments magnétiques dans le composé co-dopé d'une valeur de 6 18. Malgré divers dopants, la 
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valeur du moment magnétique total (Mror) est restée constante au sein de la super cellule et est 
principalement générée à partir de l'élément Eu, avec des contributions mineures des atomes de 


métaux de transition (Cu, Ag et Au) et de la région interstitielle. 


Un électron de TM ; Cu (3d!° 4s!), Ag (4d!° 4s!), Au (5d!° 6s!) et un électron de Eu (4f? 
6s?) forment des liaisons avec deux atomes O voisins (2 p* 2 s?). En revanche, les six électrons 
restants de Eu contribuent aux moments magnétiques du composé. Cela démontre que le TM ne 
contribue pas directement au moment magnétique, mais influe plutôt sur d'autres ions, qui 
influencent le moment magnétique résultant dans le système M£gO co-dopé. En conséquence, 
l’impureté Eu est cruciale pour améliorer les propriétés magnétiques de MgO. Comme l'illustre la 
figure V. 5, les composés semi-métalliques présentent une polarisation de spin à 100 %. Pour le 
composé Euo.125Auo.125Mg0.7500, la P % est égal à zéro (0 %) en utilisant GGA-WC car les métaux 


ont une bande interdite nulle. 


Tableau V. 3 : Valeurs des moments magnétiques totaux et locaux à l'aide des approximations 
GGA-WC et WC-mBJ, pour les composés FM Euo.125Au0.125Mg0.7500, dans les sphères à Muffin- 


Tin et dans les sites interstitiels. 


Euc.125Cu0.125Mg0.7500 WC-mBJ 6 0.0006 -0.0298 -0.031 6.180 0.186 


GGA-WC 6 0.0013 -0.2948  -0.2399 6.1623 0.2454 


Eur.5Ag 0125Mg 0.7500 WC-mBJ 6 0.0004 -0.019 -0.149 6.068 0.197 


GGA-WC 6 0.0011 -0.039 -0.1602 6.1263 0.268388 


Euo.25AU0.12:M go.7500 WC-mBJ 6 0.0010  -0.0080 -0.233 6.215 0.195 


GGA-WC 5.94 0.00183 -0.0470  -0.2569 6.2487 0.2378 


V.3. 4. Propriétés optiques 


L'analyse de la réponse des composés au rayonnement électromagnétique permet 
d'évaluer leurs propriétés optiques pour le développement de dispositifs optoélectroniques. Comme 
on le sait, les propriétés optiques de la matière sont déterminées par le tenseur de la fonction 
diélectrique £& (@), qui décrit la structure de bande et les caractéristiques optiques des matériaux et 
constitue l'interface entre la transition électronique et la structure électronique des solides. Il existe 
deux contributions à £ (@), à savoir les transitions intrabande et interbande, responsables de la 


réponse optique dans les métaux et les semi-conducteurs, respectivement. Dans la zone Brillouin, 
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les transitions interbandes peuvent être divisées en transitions directes et indirectes d'un état occupé 
à un état inoccupé, seuls les métaux sont significatifs pour la contribution des transitions 


intrabandes. 


La fonction diélectrique £ (©) = £1 (@) + 1 2 (w), qui reflète la réponse du matériau à un 
champ électrique, se compose d'une composante réfléchissante, £1 (w), la partie réelle, et d'une 
composante absorbante, £2 (w), la partie imaginaire de la fonction diélectrique. £1 (@) et £2 (©) sont 
fournis respectivement par le calcul théorique des éléments de la matrice de moment des états 


occupés et inoccupés et par les relations de dispersion de Kramers-Kronig [9-10]. 


Les paramètres optiques, telles que l'indice de réfraction, le coefficient d'absorption et le 
coefficient de réflexion, sont dérivées de la fonction diélectrique déjà calculée. Les propriétés 
optiques de tous les composés sont étudiées à l'aide des approximations GGA-WC et TB-mBJ dans 
le domaine d'énergie des photons allant de 0 à 20 eV au-dessus du niveau de Fermi (Figures V. 6 


a V. 8). 


Après co-dopage et en utilisation l'approximation WC-mB)J, les valeurs de la fonction 
diélectrique statique sont de 2.02, 1.90 et 1.85 pour Euo:125Cuo.125M 807500, Euo.125A go.125Mg0.7500, 
Euo.:125Au0.125Mg0.7500 respectivement. Ils sont inférieurs à ceux du MgO non dopé. Ceci est en 
accord avec le modèle de Penn [11], où la constante diélectrique statique £1 (0) est fortement liée à 
la bande interdite, qui diminue de la même manière. Les maximums sont décalés vers des énergies 
plus basses, de Cu à Ag puis à Au et l'énergie de coupure est absente. Ici aussi, on assiste à une 
évolution vers les hautes énergies grâce au GGA-WC. Des tendances similaires ont été observées 


dans d’autres oxydes co-dopés précédemment étudiés [12-14]. 


La figure V.7 représente la partie imaginaire de la fonction diélectrique, calculée à l'aide 
des approximations GGA-WC et WC-mBJ pour tous les composés étudiés. La partie imaginaire 
des courbes après co-dopage de M£g0O, 2 (@) des composés Eu : TM (Cu, Ag et Au)/M£gO démontre 
un déplacement du seuil de transition, en utilisant l'approximation GGA-WC à 1.18 et 1.15 pour 
Cu et Ag, respectivement. En revanche, il révèle un maximum vers les basses énergies (4.63, 5.67 
pour Cu, et 1.62, 4.91 pour Ag). Avec l'approximation WC-mBJ, £2 (@) montre un déplacement 
des seuils de transition (2.82, 3.08 et 1.91 pour Cu, Ag et Au) et des maximas vers les basses 
énergies à 5.10, 6.19 pour Cu, 4.93, 6.51 pour Ag, et 3.95, 4.77 pour Au. Les courbes semblent être 


similaires. 
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Figure V. 6. Partie réelle de la fonction diélectrique pour Euo.125Cuo.125Mg07500 (TM : Cu, Ag et 
Au) en utilisant les approximations GGA-WC et WC-mBJ. 
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Figure V. 7. Partie imaginaire de la fonction diélectrique pour Euo.1251Mo.125M£g0.7500 (TM : Cu, 


Ag et Au) en utilisant les approximations GGA-WC et WC-mBJ. 


La valeur du coefficient d'absorption reflète la capacité du matériau à générer des paires 
électron-trou et décrit le transfert d'énergie de l'onde électromagnétique vers le matériau. 
L'absorption définit jusqu'où la lumière d'un matériau spécifique d'une longueur d'onde particulière 
peut pénétrer dans un matériau avant d'être absorbée. L'absorption optique joue un rôle important 
en influençant l'activité photo catalytique. Cette propriété peut être calculée à partir de la partie 


imaginaire de la fonction diélectrique [15, 16]. 


Le coefficient d'absorption calculé des composés co-dopés, en fonction de l'énergie, est 


illustré à la figure V. 8. MgO absorbe la lumière dans la gamme UV selon les approximations 
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GGA-WC et WC-mBJ. Les transitions électroniques de la bande de valence O-2p vers la bande de 
conduction M£-3s le démontrent. Pour les composés co-dopés Eu:TM (Cu, Ag et Au)M£0O et en 
utilisant l'approximation GGA-WC, les pics d'absorption d'intensité 11.33, 24.54 et 78.49 sont 
situés à 4.74, 5.70 et 8.93 eV, respectivement, dans le cas de Euo125Cuo125Mg0 7500, 
Euo.:125A go.125Mg0.7500 a des pics situés à 4.85, 6.02 et 8.72 avec des intensités de 14.70, 17.17 et 


69.60, respectivement. Le caractère métallique de Euo125Auo125Mg0 7500 n'a pas été traité. 


Le même schéma avec un déplacement des pics vers les hautes énergies est observé suite 
à l'approximation WC-mBJ pour Cu et Ag co-dopés. Les résultats de Euo.125Auo.125Mg07500 sont 
aux positions 4.77, 6.70 et 7.93, avec des pics d'intensité d'absorption autour de 9.20, 17.92 et 
32.97. Après co-dopage, Eu : TM (Au, Ag et Cu)/ M£gO présentent un décalage vers la lumière 
visible. Le décalage du bord d'absorption se produit en raison de changements significatifs dans les 
structures électroniques dus à l'introduction de niveaux d'énergie d'impuretés dans la bande de 
valence. Les résultats suggèrent que les composés co-dopés pourraient constituer des matériaux 
potentiels pour des applications optoélectroniques et photo catalytiques, car ils présentent une 


absorption optique de la lumière visible et ultraviolette. 
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Figure V. 8. Coefficient d’absorption pour Euo.1251Mo.125Mg0.7500 (TM : Cu, Ag et Au) en 
utilisant les approximations GGA-WC et WC-mBJ. 


V. 4. Conclusion 


En employant la méthode ab initio, les propriétés optoélectroniques et magnétiques de la 


structure électronique de MgO co-dopé avec Eu-TM, adoptant la structure de type NaCI, sont 
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étudiées. Les calculs sont effectués en prenant en compte les effets de spin polarisé dans le cadre 
de l'approximation du gradient généralisé adoptée pour le potentiel d'échange-corrélation. Le terme 
TB-mBJ est intégré au terme d'échange-corrélation pour une exploration plus précise des propriétés 


optiques et magnétoélectriques. 


Les résultats de l'optimisation structurale du M£0O co-dopé révèlent une augmentation des 
constantes de réseau et une diminution du module de compressibilité avec l'augmentation du 
numéro atomique des métaux de transition, favorisant ainsi la stabilité dans la phase magnétique. 
De plus, les valeurs d'énergie de formation confirment la stabilité de la structure cubique Fm3m, 


la plus basse énergie de formation étant observée pour Euo.125Cuo.125Mg07500. 


Pour les structures de bandes électroniques, la valeur du gap au niveau de Fermi suggère 
la demi-métallicité ferromagnétique de Euo.125TMo.125Mg07500 avec une large bande interdite dans 
le canal à spin minoritaire, en utilisant l'approche TB-mBJ. Les diagrammes de densités d'états 
partielles montrent que Eu-4f, TM-d et O-2p présentent une hybridation évidente dans la région de 


Fermi. 


De plus, une étude du magnétisme est menée, révélant une origine principale de 
l'interaction p-d-f entre les orbitales des métaux de transition (TM), des terres rares (Eu), et les 
orbitales p de l'oxygène (O0). Un moment magnétique entier de 6 u8 confirmant la caractéristique 


des matériaux demi-métalliques. 


Enfin, les propriétés optiques sont été étudiées, révélant un potentiel évident dans les 
domaines des LED et des détecteurs UV. Ces résultats constituent un point de départ prometteur 


pour explorer davantage leur potentiel, tant sur le plan expérimental que théorique. 
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Conclusion générale 
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L‘oxyde de magnésium, en tant que matériau transparent dans le visible, suscite un intérêt 
croissant pour ses applications potentielles en optoélectronique. Sa bande interdite étendue, dans 
son état pur, en fait un candidat attrayant pour des technologies telles que les écrans plats et les 
LED. Lorsqu'il est dopé avec des métaux de transition, M£gO acquiert des propriétés magnétiques, 
optiques et de transport uniques, le transformant ainsi en un matériau multifonctionnel adapté à 
une gamme variée d'appareils innovants. Les recherches continues et les avancées dans le domaine 
du dopage et de la synthèse, ouvrent la voie, dans le futur, à des développements révolutionnaires 


notamment dans le domaine de l'optoélectronique et d'autres domaines connexes. 


Cette thèse vise à améliorer les propriétés physiques de MgO, un oxyde à large bande 
interdite, en le dopant et en le co-dopant avec des métaux de transition (Ni, Pd, Pt, Cu, Ag, et Au) 
et un élément des terres rares (Eu). Pour ce faire, une étude théorique approfondie est effectuée en 
utilisant la méthode d’ondes planes augmentées linéarisées à plein potentiel (FP-LAPW). Cette 
méthode est basée sur la théorie de la fonctionnelle de la densité polarisée en spin (SP-DET) et 
l'approximation du gradient généralisé (GGA), pour décrire le potentiel d'échange et de corrélation. 
De plus et pour affiner les résultats des propriétés optoélectroniques, le potentiel d'échange modifié 


de Becke-Johnson (TB-mBJ) est également utilisé. 


La première étape de cette recherche est une étude comparative des propriétés structurales, 
élastiques, électroniques, magnétiques et optiques de l’oxyde MgO dans la structure de type NaCI 
(B1). Cette comparaison vise à identifier la méthode appropriée pour l'étude de ces composés. Les 
deux approximations : la densité locale (LDA) et le gradient généralisé (GGA) sont utilisées pour 
le traitement du potentiel d'échange et de corrélation. Les résultats obtenus avec la GGA, pour les 
paramètres structuraux (constantes de réseau et modules de compressibilité), sont cohérents avec 
les données expérimentales disponibles de la littérature. L'analyse et le calcul des paramètres 
élastiques, tels que le module de Young (E), le module de cisaillement (G), le rapport B/G, le 
coefficient de poisson (v) et le facteur d'anisotropie (A) révèlent que MgO est de nature anisotrope 
et est mécaniquement stable et fragile. L'analyse des structures de bandes et de densité d'états a 
permis de confirmer que l'oxyde de magnésium (M£gO) est un matériau isolant avec une bande 
interdite directe. L'utilisation de l’approche TB-mBJ a permis d'obtenir des valeurs des gaps 


énergétiques en bon accord avec les résultats expérimentaux. Pour les propriétés optiques, la GGA 
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a induit un décalage vers les bases energies dans la zone UV dans l'absorption. Cette concordance 


valide l'utilisation de ces potentiels pour l'étude des propriétés électroniques du MgO 


La deuxième étape de la recherche consiste à étudier l'effet du dopage de M£O avec des 
métaux de transition : le nickel (Ni), le palladium (Pd) et le platine (Pt). En faisant varier le numéro 
atomique (2), tandis que la concentration est fixée à 12.5 %. L'objectif est d'évaluer l'influence de 
ces dopants sur les propriétés du M£gO. Les composés dopés X0.125Mg0.8750 (X = Ni, Pd, Pt) 
présentent un comportement semi-conducteur ferromagnétique. L’analyse des PDOS, révèle une 
contribution significative des orbitales 3d, 4d et 5d des dopants à la structure électronique du 
matériau. Cette contribution est à l'origine de la modification de bandes interdites dans les canaux 
de spin up et down autour du niveau de Fermi. Le dopage a pour effet d'augmenter 
considérablement les moments magnétiques par rapport à la matrice hôte MgO. Les moments 
magnétiques total des super cellules Nio.125Mgo.8750, Pdo.125Mg0.8750, Pto125Mgo.8750 sont égaux 
respectivement de 2, 2, et 1.99 18. Cette augmentation est due à la contribution des moments 
magnétiques localisés des dopants dus à leurs électrons 3d, 4d ou 5d non appariés et à leur 
interaction d'échange avec les électrons du M£O. Les paramètres optiques tels que la constante 
diélectriques (réelles et imaginaires), et l’absorption des composés sont calculés. Les résultats 
montrent une large absorption dans le domaine UV. Cette absorption est attribuée aux transitions 


électroniques entre les bandes nouvellement créées par le dopage. 


La deuxième partie de la recherche se focalise sur une analyse approfondie des effets 
spécifiques du dopage du M£gO avec trois métaux de transition : le cuivre (Cu), l'argent (Ag) et l'or 
(Au). L'influence de l'incorporation de ces dopants et les variations de leurs numéros atomiques et 
de leurs concentrations sont étudiées en détail. Pour ce faire, les propriétés structurales, 
électroniques et magnétiques d’une série de composés TM;M21:0 à différentes concentrations [x 


= 0%, 3.125% (SC), 6.25 % (BCC), 12.5 % (FCC) et 25 % (FCC)] sont étudiés 


L'étude des paramètres de réseau du M£gO dopé avec des métaux de transition (TM) révèle 
des tendances distinctes en fonction du numéro atomique et de la concentration du dopant. Pour les 
dopants Ag et Au, une expansion du paramètre de réseau est observée, en raison du plus grand 
rayon ionique de l'Ag et de l'Au par rapport au Mg, avec une augmentation tant du numéro atomique 
que de la concentration. En revanche, le dopage avec le Cu présente un comportement opposé, avec 


une contraction du paramètre de réseau en fonction de la concentration, attribuée à la différence de 
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taille entre les ions Cu et Mg ainsi que l'effet induite par le dopage dans la structure cristalline du 
M£0. L'étude révèle que les composés nouvellement formés après le dopage du MgO avec des 


métaux de transition sont structuralement stables en phase ferromagnétique. 


L'étude des structures électroniques des composés TM;M£g:1,0, (TM=Cu, Ag et Au, x = 
0.125, 0.062 et 0.031) révèle qu'ils sont semi-métalliques avec un comportement ferromagnétique 
et un moment magnétique total de 1 8 ainsi qu’ils sont 100 % polarisés en spin. Ces observations 
permettent de conclure que les atomes d'impuretés TM (Cu, Ag et Au) sont les principaux 
contributeurs au moment magnétique total des composés. En effet, MgO pur est intrinsèquement 
non magnétique, et l'introduction des dopants TM induit le ferromagnétisme. Avec l'augmentation 
de la concentration d'impuretés et à une concentration de 25 %, ils se comportent comme des 
métaux et l'espace ferromagnétique semi-métallique se rétrécit. Il convient de noter que les 
composés présentent une transition en fonction de la concentration de dopants ; il y a passage d'un 
semi-conducteur direct (x = 0%) à un semi-métal direct (6.25% (BCC) et 12.5% (FCC)) puis à un 
semi-métal indirect (3.125% (SC)). 


L'hybridation p-d entre les états 3d -4d et 5d TM et 2p-O à une concentration de 12,5 % 
(FCC) est responsable de la demi-métallicité et du ferromagnétisme induit dans les super cellules. 
À des concentrations plus faibles, soient 3.125 % (SC) et 6.25 % (BCC), les états 3d, 4d et 5d 
génèrent le comportement semi-métallique. Les propriétés magnétiques notamment les constantes 
d'échange No et Nog sont également étudiées et confirment le comportement ferromagnétique. 
L'analyse de la densité de charge pour le composé Ag;M£:1,0 indique l'existence de deux types de 
liaisons entre les atomes : une liaison covalente avec Ag-O, Mg-Ag ainsi qu’une liaison ionique 


avec M£g-O. 


La dernière partie de cette étude explore l'impact du co-dopage du MgO avec de 
l'europium (Eu) associé à l’un des métaux de transition (Cu, Ag et Au) à une concentration de 
12.5% (pour chacun des dopants) sur ses propriétés physiques. Cette approche vise à améliorer les 
propriétés du MgO de manière synergique en exploitant les caractéristiques individuelles de chaque 
dopant. Les résultats de l'optimisation structurale montrent que les constantes de réseau 
augmentent, tandis que le module de masse diminue. Le module compressibilité des DMO diminue 


avec l'augmentation du numéro atomique des TM. Les oxydes magnétiques dilués résultants sont 
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stables dans la phase ferromagnétique polarisée en spin, et les valeurs d'énergie de formation 


confirment leur stabilité dans la structure cubique Fm3m, 


De même, les calculs des propriétés électroniques démontrent que les trois composés 
présentent un comportement semi-métallique avec une bande interdite indirecte à L-T' égale à 2.53, 
2.49 et 1.39 eV pour Euo125Cuo125Mg07500, Euo.125Ago.125Mg07500 et Euo:125Au0.125Mg07500, 
respectivement. Ils possèdent un moment magnétique total entier de 6 18, confirmant leur nature 
ferromagnétique. L'analyse de la densité d'états révèle que les états 4f-Eu, d-TM et p-O contribuent 
à la création de nouveaux niveaux d'énergie au niveau de Fermi. Cette hybridation f-p-d est à 
l'origine des propriétés ferromagnétiques observées dans ces composés. En outre, la fonction 
diélectrique et les spectres d'absorption des composés co-dopés sont calculés et analysés. Les 
résultats prédisent un comportement optique intéressant avec une absorbance optique plus élevée 


que M£0O pur et dopé. 


L'exploration du dopage et du co-dopage de MgO avec des métaux de transition et 
l’europium ouvre des perspectives prometteuses pour la conception et la réalisation de matériaux 
(DMO) aux fonctionnalités avancées. La recherche dans ce domaine est susceptible de 
révolutionner le domaine de la spintronique, de l'optique et d'autres technologies basées sur les 


propriétés magnétiques et optiques des matériaux. 


En perspective, poursuivre l’étude en explorant d’autres propriétés telles que les propriétés 
thermodynamiques, thermoélectriques et les phonons. Il est également envisagé de modifier 
l'élément des terres rares en variant les concentrations. Une compréhension globale de cet oxyde, 
à la recherche de nouveaux matériaux solides multifonctionnels, serait ainsi obtenue. Ces 
investigations plus approfondies devraient permettre une meilleure compréhension des matériaux 
étudiés et ouvrir de nouvelles perspectives pour leur utilisation dans des applications 


technologiques innovantes. 
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